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The physics of quantum degenerate atomic Fermi gases in uniform as well as in harmonically 
trapped configurations is reviewed from a theoretical perspective. Emphasis is given to the effect of 
interactions which play a crucial role, bringing the gas into a superfluid phase at low temperature. 
In these dilute systems interactions are characterized by a single parameter, the s-wave scattering 
length, whose value can be tuned using an external magnetic field near a broad Feshbach resonance. 
The BCS limit of ordinary Fermi superfluidity, the Bose- Einstein condensation (BEC) of dimers and 
the unitary limit of large scattering length are important regimes exhibited by interacting Fermi 
gases. In particular the BEC and the unitary regimes are characterized by a high value of the 
superfluid critical temperature, of the order of the Fermi temperature. Different physical properties 
are discussed, including the density profiles and the energy of the ground-state configurations, the 
momentum distribution, the fraction of condensed pairs, collective oscillations and pair breaking 
effects, the expansion of the gas, the main thermodynamic properties, the behavior in the presence 
of optical lattices and the signatures of superfluidity, such as the existence of quantized vortices, the 
quenching of the moment of inertia and the consequences of spin polarization. Various theoretical 
approaches are considered, ranging from the mean-field description of the BCS-BEC crossover to 
non-perturbative methods based on quantum Monte Carlo techniques. A major goal of the review 
is to compare the theoretical predictions with the available experimental results. 
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I. INTRODUCTION 

An impressive amount of experimental and theoreti- 
cal activity has characterized the last ten years of ultra- 
cold atom physics. The first realization of Bose-Einstein 
condensation (BEC) in dilute vapours of alkali atoms 
(Anderson et ai, 1995; Bradley et ai, 1995; Davis et 
al, 1995) has in fact opened new stimulating perspec- 
tives in this area of research. Most of the studies in 
the first years have been devoted to quantum gases of 
bosonic nature and were aimed to investigate the impor- 
tant consequences of Bose-Einstein condensation which, 
before 1995, had remained an elusive and inaccessible 
phenomenon. 

Major achievements of these studies have been, among 
others, the investigation of superfluid features, includ- 
ing the hydrodynamic nature of the collective oscillations 
(Jin et ai, 1996; Mewes et ai, 1996), Josephson-like ef- 
fects (Cataliotti et al, 2001; Albiez et al, 2005) and the 
realization of quantized vortices (Matthews et ai, 1999; 
Madison et al, 2000; Abo-Shaeer et al, 2001); the ob- 
servation of interference of matter waves (Andrews et at, 
1997); the study of coherence phenomena in atom laser 



configurations (Mewes et al., 1997; Anderson and Kase- 
vich, 1998; Bloch, Hansen and Esslinger, 1999; Hagley 
et ai, 1999), the observation of four- wave mixing (Deng 
et ai, 1999) and of the Hanbury-Brown— Twiss effect 
(Schcllckcns et ai, 2005): the realization of spinor con- 
densates (Stenger et ai, 1998); the propagation of soli- 
tons (Burger et ai, 1999; Denschlag et ai, 2000; Strecker 
et ai, 2002; Khaykovich et ai, 2002) and the obser- 
vation of dispersive schock waves (Dutton et ai, 2001; 
Hoefer et ai, 2006); the transition to the Mott-insulator 
phase (Greiner et ai, 2002), the observation of interac- 
tion effects in the Bloch oscillations (Morsch et ai, 2001) 
and of dynamic instabilities in the presence of moving 
optical lattices (Fallani et ai, 2004); the realization of 
low dimensional configurations, including the ID Tonks- 
Girardeau gas (Kinoshita, Wenger and Weiss, 2004; Pare- 
des et ai, 2004) and the Berezinskii-Kosterlitz-Thouless 
phase transition in 2D configurations (Hadzibabic et al, 
2006; Schweikhard, Tung and Cornell, 2007). 

On the theoretical side the first efforts were devoted 
to implement the Gross-Pitaevskii theory of weakly in- 
teracting Bose gases in the presence of the trapping con- 
ditions of experimental interest. This non-linear, mean- 
field theory has proven capable to account for most of 
the relevant experimentally measured quantities in Bose- 
Einstein condensed gases like density profiles, collective 
oscillations, structure of vortices etc. The attention of 
theorists was later focused also on phenomena which can 
not be accounted for by the mean-field description like, 
for example, the role of correlations in low dimensional 
and in fast rotating configurations as well as in deep op- 
tical lattices (for general reviews on Bose-Einstein con- 
densed gases see Inguscio, Stringari and Wieman, 1999; 
Dalfovo et ai, 1999; Leggett, 2001; Pethick and Smith, 
2002; Pitaevskii and Stringari, 2003). 

Very soon the attention of experimentalists and theo- 
rists was also oriented towards the study of Fermi gases. 
The main motivations for studying fermionic systems are 
in many respects complementary to the bosonic case. 
Quantum statistics plays a major role at low tempera- 
ture. Although the relevant temperature scale providing 
the onset of quantum degeneracy is the same in both 
cases, being of the order of ksTdeg ~ h 2 n 2 ^ 3 /m where n 
is the gas density and m is the mass of the atoms, the 
physical consequences of quantum degeneracy are differ- 
ent. In the Bose case quantum statistical effects are as- 
sociated with the occurrence of a phase transition to the 
Bose-Einstein condensed phase. Conversely, in a non- 
interacting Fermi gas the quantum degeneracy tempera- 
ture only corresponds to a smooth crossover between a 
classical and a quantum behavior. Differently from the 
Bose case the occurrence of a superfluid phase in a Fermi 
gas can be only due to the presence of interactions. From 
the many-body point of view the study of Fermi superflu- 
idity opens a different and richer class of questions which 
will be discussed in this review paper. Another impor- 
tant difference between Bose and Fermi gases concerns 
the collisional processes. In particular, in a single com- 
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FIG. 1: Gallery of molecular BEC experiments. Bimodal spatial distributions were observed for expanding gases at JILA 
(Greiner, Regal and Jin, 2003) with 40 K, at MIT (Zwierlein et al, 2003) and at ENS (Bourdel et al., 2004) with 6 Li. They 
were instead measured in situ at Innsbruck (Bartenstein et al., 2004a) and at Rice University (Partridge et al., 2005) with 6 Li. 



ponent Fermi gas s-wave scattering is inhibited due to 
the Pauli exclusion principle. This effect has dramatic 
consequences on the cooling mechanisms based on evap- 
oration, where thermalization plays a crucial role. This 
has made the achievement of low temperatures in Fermi 
gases a difficult goal that was eventually realized with the 
use of sympathetic cooling techniques either employing 
two different spin components of the same Fermi gas or 
adding a Bose gas component as a refrigerant. 

First important achievements of quantum degeneracy 
in trapped Fermi gases were obtained by the group at 
JILA (De Marco and Jin, 1999). In these experiments 
temperatures of the order of fractions of the Fermi tem- 
peratures were reached by working with two spin com- 
ponents of 40 K atoms interacting with negative scatter- 
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FIG. 2: Transition temperature in units of the Fermi energy 
Ef as a function of the interaction strength along the BCS- 
BEC crossover, calculated using BCS mean-field theory (from 
Sa de Melo, Randeria and Engelbrecht, 1993). The diamond 
corresponds to the theoretical prediction by Burovski et al. 
(2006a) based on a Quantum Monte Carlo simulation at uni- 
tarity. 



ing length. According to the Bardeen-Cooper-Schrieffcr 
(BCS) theory this gas should exhibit superfluidity at suf- 
ficiently low temperature. However, due to the extreme 
diluteness of the gas, the critical temperature required 
to enter the superfiuid phase was too small in these ex- 
periments. Quantum degeneracy effects were later ob- 
served in 6 Li Fermi gases (Truscott et al., 2001; Schreck 
et al., 2001) using sympathetic cooling between 6 Li and 
the bosonic 7 Li isotope. Fermion cooling using different 
bosonic species has also proved very efficient for instance 
in the case of 40 K- 87 Rb (Roati et al, 2002) as well as 
6 Li- 23 Na (Hadzibabic et al, 2003). 

It was soon realized that a crucial tool to achieve su- 
perfluidity is provided by the availability of Feshbach res- 
onances. These resonances characterize the two body in- 
teraction and permit one to change the value and even the 
sign of the scattering length by simply tuning an exter- 
nal magnetic field. Feshbach resonances were first inves- 
tigated in bosonic sytems (Courteille et al., 1998; Inouye 
et al, 1998). However, inelastic processes severely limit 
the possibility of tuning the interaction in Bose conden- 
sates (Stenger et al., 1999). Strongly interacting regimes 
of fermionic atoms were achieved by O'Hara et al. (2002) 
and Bourdel et al. (2003) working at the resonance where 
the scattering length takes a divergent value. In this case 
three-body losses are inhibited by the Pauli exclusion 
principle, leading to a greater stability of the gas (Petrov, 
Salomon and Shlyapnikov, 2004). The resonant regime, 
also called unitary regime, is very peculiar since the gas 
is at the same time dilute (in the sense that the range of 
the interatomic potential is much smaller than the inter- 
particle distance) and strongly interacting (in the sense 
that the scattering length is much larger than the inter- 
particle distance). All the length scales associated with 
interactions disappear from the problem and the system 
is expected to exhibit a universal behavior, independent 
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of the details of the interatomic potential. Bertsch (1999) 
and Baker (1999) first discussed the unitary regime as 
a model for neutron matter based on resonance effects 
in the neutron-neutron scattering amplitude (for a re- 
cent comparison between cold atoms and neutron matter 
see Gezerlis and Carlson, 2007). The critical tempera- 
ture of the new system is much higher than in the BCS 
regime, being of the order of the quantum degeneracy 
temperature, which makes the realization of the super- 
fluid phase much easier. Thanks to Feshbach resonances 
one can also tune the scattering length to positive and 
small values. Here, bound dimcrs composed of atoms of 
different spin are formed and consequently the system, 
that was originally a Fermi gas, is transformed into a 
bosonic gas of molecules. The possibility of tuning the 
scattering length across the resonance from negative to 
positive values, and vice- versa, provides a continous con- 
nection between the physics of Fermi superfluidity and 
Bose-Einstein condensation, including the unitary gas as 
an intermediate regime. 

On the BEC side of the Feshbach resonance molecules 
can be created either by directly cooling the gas at pos- 
itive values of the scattering length a, or by first cool- 
ing the gas on the BCS side and then tuning the value 
of a through the resonance. At low enough tempera- 
tures, Bose-Einstein condensation of pairs of atoms was 
observed through the typical bimodal distribution of the 
molecular profiles (see Fig. [1]) (Greiner, Regal and Jin, 
2003; Jochim et al. 2003; Zwierlein et al, 2003; Bourdel 
et al., 2004; Partridge et al., 2005). Condensation of pairs 
was later measured also on the fermionic side of the res- 
onance (Regal, Greiner and Jin, 2004b; Zwierlein et al., 
2004). Other important experiments have investigated 
the surprisingly long lifetime of these interacting Fermi 
gases (Strecker, Partridge and Hulet, 2003; Cubizolles et 
al., 2003), the release energy (Bourdel et al., 2004) and 
the density profiles (Bartenstein et al., 2004a) along the 
crossover. 

Many relevant experiments have also focused on the 
dynamic behavior of these interacting systems, with the 
main motivation of exploiting their superfluid nature. 
The first observation of anisotropic expansion (O'Hara 
et al., 2002) and the measurements of the collective os- 
cillations (Bartenstein et al., 2004b; Kinast et al., 2004), 
although nicely confirming at low temperatures the pre- 
dictions of the hydrodynamic theory of superfluids, can 
not be however considered a proof of superfluidity since a 
similar behavior is prediced also in the collisional regime 
of a normal gas above the critical temperature. The mea- 
surement of the pairing gap observed in radio-frequency 
excitation spectra (Chin et al., 2004) was an important 
step toward the experimental evidence of superfluidity, 
even though it was not conclusive since pairing correla- 
tions are present also in the normal phase. A convincing 
proof of superfluidity was actually provided by the obser- 
vation of quantized vortices which were realized on both 
sides of the Feshbach resonance (Zwierlein et al., 2005b). 

More recent experimental work (Zwierlein et al. 2006a; 



TABLE I: Ratio T c /Tf of the transition temperature to the 
Fermi temperature in various Fermi superfluids. 





T C /T F 


Conventional superconductors 
Superfluid 3 He 

High-temperature superconductors 
Fermi gases with resonant interactions 


10" 5 -10" 4 

icr 3 

10" 2 
~0.2 



Partridge et al. 2006a and 2006b) has concerned the 
study of spin polarized configurations with an unequal 
number of atoms occupying two different spin states. 
In particular the Clogston-Chandrasekar limit, where 
the system looses superfluidity, has been experimentally 
identified at unitarity (Shin et al., 2006). These config- 
urations provide the unique possibility of observing the 
consequences of superfluidity through sudden changes in 
the shape of the cloud as one lowers the temperature, in 
analogy to the case of Bose-Einstein condensates (Zwier- 
lein et al., 2006b). Another rapidily growing direction of 
research is the study of Fermi gases in periodic poten- 
tials (Modugno et al, 2003). First experimental results 
concern the effect of the periodic lattice on the bind- 
ing energy of molecules across the Feshbach resonance 
(Stoferle et al, 2006) and some aspects of the superfluid 
to Mott insulator transition (Chin et al., 2006). A moti- 
vation of the investigations in this field is the possibility 
of implementing an important model of condensed mat- 
ter physics, the Hubbard hamiltonian, in analogy to the 
Bose-Hubbard counterpart already realized in bosonic 
systems (Greiner et al., 2002). Fermi gases in periodic 
potentials are of high interest also in the absence of inter- 
actions. For example they give rise to long-living Bloch 
oscillations that were observed in spin polarized Fermi 
gases (Roati et al., 2004). 

On the theoretical side the availability of interacting 
Fermi gases with tunable scattering length has stimulated 
an impressive amount of work. Differently from the case 
of dilute Bose gases, where the Gross-Pitaevskii equa- 
tion provides an accurate description of the many-body 
physics at low temperature and small densities, an analog 
theory for the Fermi gas along the BCS-BEC crossover is 
not available. The theoretical efforts started in the con- 
text of superconductors with the work by Eagles (1969), 
where it was pointed out that for large attraction between 
electrons the equations of BCS theory describe pairs of 
small size with a binding energy independent of density. 
A thorough discussion of the generalization of the BCS 
approach to describe the crossover in terms of the scatter- 
ing length was presented in the seminal paper by Leggett 
(1980) (see also Leggett, 2006). This work concerned 
ground-state properties and was later extended to finite 
temperatures by Nozieres and Schmitt-Rink (1985) and 
by Sa de Melo, Randeria and Engelbrecht (1993) to calcu- 
late the critical temperature for the onset of superfluidity. 
These theories describe the properties of the many-body 
configurations along the BCS-BEC crossover in terms of 
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a single parameter related to interactions, the dimension- 
less combination kpa where kp is the Fermi wavevector. 
In Fig. [2] we report theoretical predictions for the criti- 
cal temperature, showing that T c is of the order of the 
Fermi temperature in a wide interval of values of /cp|a|. 
For this reason one often speaks of high-T c Fermi super- 
fluidity (see Table [T| . Furthermore, the results shown in 
Fig. [U suggest that the transition between BCS and BEC 
is indeed a continous crossover. 

The first application of the concept of bound pairs to 
the case of Fermi gases with resonant interactions was 
proposed by Timmermans et al. (2001) and Holland 
et al. (2001). These extensions of the BCS mean-field 
theory are however approximate and, even at zero tem- 
perature, the solution of the many-body problem along 
the crossover is still an open issue. Different methods 
have been developed to improve the description of the 
BCS-BEC crossover in uniform gases as well as in the 
presence of harmonic traps. These methods include the 
solution of the four-body problem to describe the inter- 
action between molecules on the BEC side of the reso- 
nance (Petrov, Salomon and Shlyapnikov, 2004), appli- 
cations of the BCS mean-field theory to trapped config- 
urations with the help of the local density approxima- 
tion, extensions of the mean-field approach using dia- 
grammatic techniques (Pieri, Pisani and Strinati, 2004; 
Chen et al., 2005; Haussmann et al., 2006) and the de- 
velopment of theories based on the explicit inclusion of 
bosonic degreees of freedom in the Hamiltonian (Ohashi 
and Griffin, 2002; Bruun and Pethick, 2004; Romans and 
Stoof, 2006). At the same time more microscopic calcula- 
tions based on quantum Monte Carlo (QMC) techniques 
have become available providing results on the equa- 
tion of state at zero temperature (Carlson et al., 2003; 
Astrakharchick et al., 2004a; Juillet, 2007) and on the 
critical temperature for the superfluid transition (Bul- 
gac, Drut and Magierski, 2006; Burovski et al., 2006a; 
Akkineni, Ceperley and Trivedi, 2006). In addition to 
the above approaches, aimed to investigate the equilib- 
rium properties of the system, a successful direction of 
research was devoted to the study of dynamic proper- 
ties, like the expansion and the collective oscillations, 
by applying the hydrodynamic theory of superfluids to 
harmonically trapped Fermi gases (Menotti, Pedri and 
Stringari, 2002; Stringari, 2004), the pair-breaking exci- 
tations produced in resonant light scattering (Torma and 
Zoller, 2000) and the dynamic structure factor (Minguzzi, 
Ferrari and Castin, 2001). A large number of theoreti- 
cal papers has been recently devoted also to the study 
of spin polarization effects, with the aim of revealing the 
consequences of superfluidity on the density profiles and 
on the emergence of new superfluid phases. 

Since the number of papers published on the subject 
of ultracold Fermi gases is very large we decided to limit 
the presentation only to some aspects of the problem 
which naturally reflect the main interests and motiva- 
tions of the authors. In particular, we have tried to 
give special emphasis on the physical properties where 



an explicit comparison between theory and experiment 
is available, focusing on the effects of the interaction and 
on the manifestations of superfluidity exhibited by these 
novel trapped quantum systems. Most of the results pre- 
sented in this review are relative to systems at zero tem- 
perature where the theoretical predictions are more sys- 
tematic and the comparison with experiments is more 
reliable. A more complete review, covering all the in- 
teresting directions of theoretical research would require 
a much bigger effort, beyond the scope of the present 
paper. Some advanced topics related to the physics of 
ultracold Fermi gases are discussed, for example, in the 
Proceedings of the 2006 'Enrico Fermi' Varenna School 
(Ketterle, Inguscio and Salomon, 2007) and in the review 
article by Bloch, Dalibard and Zwerger (2007). 

II. IDEAL FERMI GAS IN HARMONIC TRAP 
A. Fermi energy and thermodynamic functions 

The ideal Fermi gas represents a natural starting point 
for discussing the physics of dilute Fermi gases. In many 
cases the role of interactions can in fact be neglected, 
as in the case of spin polarized gases where interactions 
are strongly suppressed at low temperature by the Pauli 
exclusion principle, or treated as a small perturbation. 

The ideal Fermi gas in the harmonic potential 

Vho = \^nu? x x 2 + \^mu? y y 2 + ^mu 2 z z 2 (1) 

is a model system with many applications in different 
fields of physics, ranging from nuclear physics to the 
more recent studies of quantum dots. For this reason 
we will mainly focus on the most relevant features of the 
model, emphasizing the large N behavior where many 
single-particle states are occupied and the semiclassical 
approach can be safely used. The simplest way to intro- 
duce the semiclassical description is to use a local den- 
sity approximation for the Fermi distribution function of 
a given spin species: 

/(r ' P) = exp[p(p 2 /2m + V ho (r)- /i)] + 1 ' (2) 

where /3 — 1/ksT and fi is the chemical potential fixed 
by the normalization condition 

(2irh) 3 J Jo exp[/3(e - /z)] + 1 

(3) 

N a being the number of atoms of the given spin species 
which is supposed to be sufficiently large. The semiclas- 
sical approach accounts for the Fermi pressure at low 
temperatures. In Eq. ([3]) we have introduced the single- 
particle density of states g(e) whose energy dependence is 
given by g(e) = e 2 /2(ftw ho ) 3 where u ho = (u> x u y u z ) 1/3 is 
the geometrical average of the three trapping frequencies. 
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In terms of the density of states one can easily calculate 
the relevant thermodynamic functions. For example, the 
energy of the gas is given by the expression 



£7(T) 



dc- 



(4) 



At T = the chemical potential fj, coincides with the 
Fermi energy 



E 



HO = k B T*° = (6AT CT )V3^ o ^ 



(5) 



and the energy takes the value £7(0) = 3/AE^°N a . 
Eq. ([I} fixes an important energy (and temperature) 
scale in the problem, analog to the expression Ep = 
(fi 2 /2m)(67r 2 n cr ) 2 / 3 of the uniform gas, where n a is the 
density of a single spin component. 

It is worth noticing that the Fermi energy ([5]) has the 
same dependence on the number of trapped atoms and on 
the oscillator frequency uj) lo as the critical temperature 
for Bose-Einstein condensation given by the well known 
formula k B T BE c ^ 0.94/iaWV 1 / 3 . 

An important quantity to investigate is also the release 
energy £7 re i defined as the energy of the gas after a sud- 
den switching off of the confining potential. The release 
energy is directly accessible in time-of-flight experiments 
and, as a consequence of the equipartition theorem ap- 
plied to the ideal gas with harmonic confinement, is al- 
ways equal to £7 xe ] = £7/2, where £7(T) is the total energy 
(111) of the gas. At low T the energy per particle deviates 
from the classical value ZkgT due to quantum statistical 
effects as clearly demonstrated in the JILA experiment 
(De Marco and Jin, 1999; De Marco, Papp and Jin, 2001) 
reported in Fig. [3] 



In Eq. {7} aho = \J fr/rriLUho denotes the harmonic oscil- 
lator length. The Fermi wavevector 



fc ° = J_(48)V«j\rV« 
n aho 



(9) 



fixes instead the width of the momentum distribution 

21 3/2 



»v(p) 



N„ 



^ (p° )3 



(10) 



This result is obtained by integrating the T = dis- 
tribution function in coordinate space. Equations ([5]) 
and (|10p . which are normalized to the total number of 
particles N a , hold for positive values of their arguments 
and are often referred to as Thomas-Fermi distributions. 
Equation (fT0|) is the analogue of the most familiar mo- 
mentum distribution 3iV (T /(47rp'p)0(l — p 2 /pp) charac- 
terizing a uniform Fermi gas in trems of the Fermi mo- 
mentum pp. The broadening of the Fermi surface with 
respect to the uniform case is the consequence of the 
harmonic trapping. Notice that the value of k F de- 
fined above coincides with the Fermi wavevector k F = 
[67r 2 n cr (0)] 1 / 3 of a uniform gas with density ?v(0) calcu- 
lated in the center of the trap. It is worth comparing 
Eqs. (J5J) and (|10[) with the analogous results holding for 
a trapped Bose-Einstein condensed gas in the Thomas- 
Fermi limit (Dalfovo et al., 1999). The shapes of Fermi 
and Bose profiles do not look very different in coordi- 
nate space. In both cases the radius of the atomic cloud 
increases with N although the explicit dependence is 
slightly different (TV 1 / 5 for bosons and iV 1 / 6 for fermions). 
Notice however that the form of the density profiles has 



B. Density and momentum distributions 

The Fermi energy (O can be used to define typical 
length and momentum scales characterizing the Fermi 
distribution in coordinate and momentum space respec- 
tively: 

R° = ^E^/mcof ; p° F = ^2mEf° . (6) 
The Thomas-Fermi radius (i — x, y, z) 

i? s ° =a h0 (48N a ) 1/6 uho/ui , (7) 

gives the width of the density distribution at T = 0, 
which can be calculated by integrating the distribution 
function in momentum space: 

f \ 8 N ° 

« >-(i) ! -G|) ! -Gi) 2 " 




0.0 0.5 1.0 1.5 

T/T F 

FIG. 3: Evidence for quantum degeneracy effects in trapped 
Fermi gases. The average energy per particle, extracted from 
absorption images, is shown for two-spin mixtures. In the 
quantum degenerate regime the data agree well with the ideal 
Fermi gas prediction (solid line). The horizontal dashed line 
corresponds to the result of a classical gas. From De Marco, 
Papp and Jin (2001). 
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a deeply different physical origin in the two cases. For 
bosons it is fixed by the repulsive two-body interactions, 
while in the Fermi case is determined by the quantum 
pressure. 

In momentum space the Bose and Fermi distributions 
differ instead in a profound way. First, as a consequence 
of the semiclassical picture, the momentum distribution 
of the Fermi gas is isotropic even if the trapping poten- 
tial is deformed [see Eq. (|T0|) ] . This behavior differs from 
what happens in the BEC case where the momentum dis- 
tribution is given by the square of the Fourier transform 
of the condensate wavefunction and is hence sensitive to 
the anisotropy of the confinement. Second, the width of 
the momentum distribution of a trapped Bose-Einstein 
condensed gas decreases by increasing N while, accord- 
ing to Eqs. © and (TIT))) , the one of a trapped Fermi gas 
increases with the number of particles. 

It is finally useful to calculate the time evolution of the 
density profile after turning off the trapping potential. 
For a non-interacting gas the distribution function fol- 
lows the ballistic law /(r, p, t) = /o(r — pt/m, p), where 
fo is the distribution function at t = given by @. By 
integrating over p one can easily calculate the time evo- 
lution of the density and one finds the following result 
for the mean square radii 



(rh = 



E(T) 1 



(11) 



The asymptotic isotropy predicted by Eq. (|11|) is the con- 
sequence of the absence of collisions during the expansion 
and reflects the isotropy of the momentum distribution 

ma. 



If one neglects small relativistic spin interactions, the 
problem of describing the collision process between two 
atoms reduces to the solution of the Schrodinger equation 
for the relative motion. For positive energy e, the s-wave 
wavefunction in the asymptotic region r Rq can be 
written as ipa(r) cx sin[fer + 5o(k)]/r, where r = |ri — r^l 
is the relative coordinate of the two atoms, 5o(k) is the 
s-wave phase shift and k — ^/2m r e/h is the wavevector of 
the scattering wave with m r the reduced mass of the pair 
of atoms (m r = m/2 for identical atoms). The s-wave 
scattering amplitude, fo(k) = [— k cot S (k) + ifc]" 1 , docs 
not depend on the scattering angle and when k — ► 
it tends to a constant value: fo(k — > 0) = —a. The 
quantity a is the s-wave scattering length, which plays a 
crucial role in the scattering processes at low energy. By 
including terms to order k 2 in the expansion of the phase 
shift So(k) at low momenta one obtains the result 



/o(fc) 



1 - k 2 R*/2 + ik ' 



(13) 



defining the effective range R* of interactions. This 
length scale is usually of the same order of the range Rq, 
however in some cases, e.g. close to a narrow Feshbach 
resonance (see Sec. MI Bj) . it can become much larger than 
Rq providing a new relevant scale. In the limit a — » oo, 
referred to as "unitary limit", the scattering amplitude 
([13]) at wavevectors k <C 1/|-R*| obeys to the universal 
law fo(k) = i/k, independent of the interaction. 

For positive scattering lengths close to the resonance 
(a 3> Rq) shallow s-wave dimers of size a exist whose 
binding energy does not depend on the short-range 
details of the potential and is given by 



III. TWO-BODY COLLISIONS 

A. Scattering properties and binding energy 

Interaction effects in quantum degenerate, dilute Fermi 
gases can be accurately modeled by a small number of 
parameters characterizing the physics of two-body col- 
lisions. In the relevant regime of low temperature and 
large mean interparticle distance, the spatial range i?o of 
the interatomic potential is much smaller than both the 
thermal wavelength Xp = ^J2nh 2 jmk^T and the inverse 
Fermi wavevector kp : 



Rq <C At 



Rq < kp 1 



(12) 



Under the above conditions the main contribution to 
scattering comes from states with I = component of 
angular momentum, i.e. s-wave states. Another con- 
straint comes from the antisymmetry of the wavefunc- 
tion of identical fermions which excludes s-wave scatter- 
ing between spin-polarized particles. As a consequence, 
only particles with different spin can interact. 

In this Section we briefly recall some results of the 
theory of elastic scattering in the s-wave channel (see, 
e.g., Landau and Lifshitz, 1987). 



2m r a 2 



(14) 



The binding energy of molecules of two fermionic 40 K 
atoms formed near a Feshbach resonance was first mea- 
sured by Regal et al. (2003a) using radio frequency spec- 
troscopy. 

In the many-body treatment of interactions it is con- 
venient to use an effective potential V e f / instead of the 
microscopic potential V . Different model potentials can 
be considered as the description of low-energy processes 
is independent of the details of V(r). In many applica- 
tions one introduces the regularized zero-range pseudo- 
potential defined as (Huang and Yang, 1957) 



d 

Veff(r)=gS{r)—r 



(15) 



where the coupling constant g is related to the scattering 
length by the relationship g = 2ith 2 a/m r . This potential 
has a range Rq = and results in the scattering ampli- 
tude f(k) = —l/(a~ 1 + ik). For a > 0, a bound state 
exists having the binding energy (|14p and corresponding 
to the normalized wavefunction 



(16) 
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Notice that the differential operator (d/dr)r in (fT5| elim- 
inates the singular 1 jr short-range behavior of the wave- 
function. The use of the pseudopotential (fT5| in the 
Schrodinger equation is equivalent to the following con- 
tact boundary condition imposed on the wavefunction 
ip(r) (Bethe and Peierls, 1935) 



d(rtp)/dr 



rip 



1 



(17) 



r=0 



Another model potential which will be considered in 
this review in connection with quantum Monte-Carlo 
simulations is the attractive square- well potential defined 
by 



Ve ff (r) 



-Vo (r<R n ) 
(r > Ro) • 



(18) 



The s-wave scattering parameters can be readily calcu- 
lated in terms of the range Ro and of the wavevector 
Kq = ^2m r Vo/h 2 . The scattering length is given by 
a = Rq[1 — tan(KoRo)/ (KoRo)] and the effective range 
by R* =Ro- Rl/Za 2 - 1/K 2 a. 



two particles (referred to as open channel) becomes close 
to the energy of a bound state of the pair in a different 
spin state (closed channel). 

If the magnetic moments of the pairs of atoms in the 
two channels are different, one can go from a situation 
where the bound state in the closed channel is just below 
the threshold of the continuum spectrum in the open 
channel to a situation where the same bound state is 
just above threshold. 

The transition between the two situations takes place 
at some value (denoted by Bg) of the magnetic field. In 
the absence of coupling, the existence of the bound state 
in the closed channel has no effect on the scattering in 
the open channel. However, in the presence of small cou- 
pling induced for example by exchange interactions, the 
scattering length will be large and positive if the state is 
below threshold and large and negative in the opposite 
case. As a function of the magnetic field B the scattering 
length can be parametrized in the following form 



a bg 1 



B-B 



(19) 



B. Fano-Feshbach resonance 

The recent experimental achievements in the field of 
ultra-cold Fermi gases are based mainly on the possibil- 
ity of tuning the scattering length a, in particular to val- 
ues much larger than the mean interatomic distance, by 
changing an external magnetic field. This situation ex- 
ists near the so called Fano-Feshbach resonances (Fano, 
1961; Feshbach, 1962). These resonances take place when 
the energy associated with the scattering process between 




0.0 0.5 1.0 1.5 2.0 
Magnetic field [kG] 



FIG. 4: Magnetic field dependence of the scattering length 
in 6 Li, showing a broad Feshbach resonance at Bo ~ 834 G 
and a narrow Feshbach resonance at Bo — 543 G (can not be 
resolved on this scale). From Bourdel et al. (2003). 



where is the width of the resonance and ab g is the 
background scattering length away from the resonance. 

An important distinction concerns broad and narrow 
resonances, which in a Fermi gas involves the comparison 
of kp and the effective range of interactions \R*\. Broad 
resonances correspond to fcp|ii*| <C 1. In this case, the 
effective range is irrelevant at the many-body level and 
the properties of the gas near the resonance can be de- 
scribed only in terms of £;f|«| (Partridge et al., 2005). 
On the contrary, for narrow resonances corresponding to 
kp\R*\ ^ 1) the effective range is negative and becomes a 
relevant scale of the problem (Bruun and Pethick, 2004; 
Bruun, 2004; De Palo et al, 2004). 

Most experiments on ultracold fermions make use of 
broad Feshbach resonances. This is certainly the case for 
the 40 K resonance at Bo ~ 202 G used in the experiments 
at JILA (Loftus et al., 2002) and even more so for the 
extremely wide 6 Li resonance at Bq = 834 G used in the 
experiments at Duke (O'Hara et al., 2002), Paris (Bour- 
del et al, 2003), Innsbruck (Jochim et al, 2003), MIT 
(Zwierlein et al., 2003) and Rice (Partridge et al., 2005). 
In both cases the value of \R*\ close to the resonance 
is of the order of or smaller than a few nanometers and 
therefore kp\R*\ < 0.01 for typical values of the density. 
Different is the case of the resonance in 6 Li at Bo — 543 
G used at Rice (Strecker, Partridge and Hulet, 2003), 
where estimates give instead kp\R*\ > 1. In Fig. 2] we 
report the predicted behavior of the scattering length in 
6 Li as a function of the external magnetic field showing 
both the broad and the narrow resonance. Notice that 
for small values of the external magnetic field the scatter- 
ing length approaches the value a = 0, where 6 Li atoms 
are expected to behave as a non-interacting gas. 
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C. Interacting dimers 

The properties of shallow dimers formed near a Fesh- 
bach resonance are very important in the physics of ul- 
tracold Fermi gases. Consisting of fermionic atoms, these 
dimers are bosonic molecules and interact with each other 
as well as with single atoms. 

The scattering between atoms and weakly-bound 
dimers was first investigated by Skorniakov and Ter- 
Martirosian (1956) in connection with neutron-deuteron 
scattering and, more recently, by Petrov (2003) in the 
context of degenerate Fermi gases (see also Brodsky et 
ai, 2006; Levinsen and Gurarie, 2006; Taylor, Griffin and 
Ohashi, 2007; Iskin and Sa de Melo, 2007b). The solu- 
tion of the three-body Schrodinger equation for a pair of 
like fermions interacting with a third particle can be ob- 
tained exactly using the contact boundary condition (|17[) 
between particles with different spin. From the behav- 
ior of the scattering solution at large separation distance 
between the dimer and the free atom one can extract 
the atom-dimer scattering length, which is found to be 
proportional to a: 

a ad ~ 1.18a. (20) 

Fermi statistics plays here a crucial role since three-body 
bound states (Efimov states) are not permitted due to 
the Pauli principle. 

Petrov, Salomon and Shlyapnikov (2004) have also 
solved the problem of the collisions between two dimers. 
By using again the zero-range approximation they calcu- 
lated the dimcr-dimcr scattering length finding the value 

a dd ~ 0.60a . (21) 

The above result was later derived also using diagram- 
matic techniques by Brodsky et al. (2006) and by Levin- 
sen and Gurarie (2006). It is worth pointing out that by 
applying Born approximation one would find the results 
a ad = 8a/3 and add = 2a (Pieri and Strinati, 2006). 

The weakly-bound dimers formed near a Feshbach 
resonance are molecules in the highest roto-vibrational 
state. Due to collisions they can fall into deeper bound 
states of size on the order of the interaction range Rq. In 
this process a large energy of order h 2 /tuRq is released 
and converted into kinetic energy of the colliding atoms, 
which then leave the system. In the case of atom-dimer 
collisions one can estimate the probability for the three 
atoms to approach each other within distances ~ Rq. 
This probability is suppressed by the Pauli principle, be- 
cause two out of the three atoms have the same spin. A 
description of the relaxation process is provided by the 
equation: n a = —a a dn a nd, where n a and rid are respec- 
tively the densities of atoms and dimers and h a is the rate 
of atom losses. For the coefficcnt a a d the following de- 
pendence on a has been obtained (Petrov, Salomon and 
Shlyapnikov, 2004): 

a ad oc (hR /m){Ro/ay . (22) 



with s = 3.33. In the case of relaxation processes caused 
by dimer-dimer collisions, the coefficient entering the 
dimer loss equation fid — — ctddn 2 ! satisfies the same law 
(f2"2")l with s = 2.55. It is crucial that both a a d and add 
decrease with increasing a. This dependence ensures the 
stability of Fermi gases near a Feshbach resonance. It is 
the consequence of the fermionic nature of the atoms. In 
the case of bosons, instead, the relaxation time increases 
with increasing a and the system becomes unstable ap- 
proaching the resonance. 

According to the above results, the dimer-dimer relax- 
ation rate should dominate over the atom-dimer one in 
the limit Ro/a — + 0. Experiments on atom losses both in 
potassium (Regal, Greiner and Jin, 2004a) and in lithium 
(Bourdel et ai, 2004) close to the Feshbach resonance 
give relaxation rate constants add oc a~ s with values of 
the exponent s in reasonable agreement with theory. 

An interesting situation takes place in the case of het- 
eronuclear dimers, consisting of fermionic atoms with 
different masses mi and rri2, where mi > m-i (see 
Sec. IDTDT The theory (Petrov, 2003; Petrov, Sa- 
lomon and Shlyapnikov, 2005) predicts that at mass ratio 
m\/m2 — 12.3 the exponent s in the dimer-dimer relax- 
ation rate add oc a~ s changes its sign, violating thus the 
stability condition of the gas near resonance. For mass ra- 
tios larger than 13.6, short-range physics dominates and 
the universal description in terms of the scattering length 
a is lost. 



IV. THE MANY-BODY PROBLEM AT 
EQUILIBRIUM: UNIFORM GAS 

A. Hamiltonian and effective potential 

The ideal gas model presented in Sec. [TT] provides a 
good description of a cold spin polarized Fermi gas. In 
this case interactions are in fact strongly inhibited by the 
Pauli exclusion principle. When atoms occupy different 
spin states interactions instead deeply affect the solution 
of the many-body problem. This is particularly true at 
very low temperature where, as we will discuss in this 
Section, attractive interactions give rise to pairing effects 
responsible for the superfluid behavior. 

Let us consider a two-component system occupying 
two different spin states hereafter called, for simplicity, 
spin-up (<7 =t) and spin-down (a =!). We consider the 
grand canonical many-body Hamiltonian written in sec- 
ond quantization as 

+ J eMr'V(r- r')*}(r)*j(r')*i(r')*T(r) , (23) 

where the field operators obey the fermionic anticommu- 
tation relations (r), #*,(r')} = ^.^(r - r'). The 
one-body potential V a ^xt and the two-body potential V 
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account for the external confinement and for the inter- 
action between atoms of different spin, respectively. The 
number of atoms N a = J dr < \E'J-(r) ( I' cr (r) >, as well as 
the trapping potentials and the atomic masses of the two 
spin species can in general be different. In this Section 
we consider the uniform case {V-\. ex t — V^ext — 0) with 
= Ni = N/2 and m^ = mi = m. The densities of 
the two spin components are in this case uniform and the 
Fermi wavevector is related to the total density of the gas, 
n = 27i| = 2nj_, through the expression kp = {?>TT 2 n) 1 / 3 ' . 
The density n determines the Fermi energy 

h 2 

E F = k B T F = — (37r 2 n) 2 / 3 (24) 
2m 

of the non-interacting gas. In the following we will al- 
ways use the above definition of Ep. Notice that in the 
presence of interactions the above definition of Ep differs 
from the chemical potential at T = 0. 

In our discussion we are interested in dilute gases where 
the range of the interatomic potential is much smaller 
than the interparticle distance. Furthermore, we assume 
that the temperature is sufficiently small so that only 
collisions in the s-wave channel are important. Under 
these conditions interaction effects are well described by 
the s-wave scattering length a (see Sec. MI A|) . In this re- 
gard one should recall that the gaseous phase corresponds 
to a metastable solution of the many-body problem, the 
true equilibrium state corresponding in general to a solid- 
phase configuration where the microscopic details of the 
interactions are important. In Eq. (|23|) we have ignored 
the interaction between atoms occupying the same spin 
state, which is expected to give rise only to minor correc- 
tions due to the quenching effect produced by the Pauli 
principle. 

As we have already discussed in connection with two- 
body physics, a better theoretical understanding of the 
role played by the scattering length can be achieved by 
replacing the microscopic potential V with an effective 
short-range potential U e //. The regularized zero-range 
pseudo-potential has been introduced in Eq. (fT5|) . Sim- 
ilarly to the two-body problem, the regularization ac- 
counted for by the term (d/dr)r permits to cure the ul- 
traviolet divergences in the solution of the Schrodinger 
equation that arise from the vanishing range of the 
pseudo-potential. In general, this regularization is cru- 
cial to solve the many-body problem beyond lowest order 
perturbation theory, as happens, for example, in the BCS 
theory of superfluidity (Bruun et ai, 1999). 

The effect of the zero-range pseudo-potential is ac- 
counted for by the boundary condition (|17[) which, in 
the many-body problem, can be rewritten as (Bcthe and 
Peierls, 1935; Petrov, Salomon and Shlyapnikov, 2004) 

tf(r« -> 0) oc — - - , (25) 
Tij a 

where ry = |r$ — Tj\ is the distance between any pair 
of particles with different spin and the limit is 



taken for fixed positions of the remaining N — 2 par- 
ticles and of the center of mass of the pair For 
realistic potentials, the above short-distance behavior is 
expected to hold for length scales much larger than the 
effective range \R*\ of the interaction and much smaller 
than the mean interparticle distance: \R*\ <C <C kp 1 . 
This range of validity applies in general to atomic gas- 
like states both with repulsive (a > 0) and attractive 
(a < 0) interactions. If the many-body state consists in- 
stead of tight dimers of size a <C A;^ 1 described by the 
wavefunction (fTtJ)) the boundary condition (j2"5|) is valid 
in the reduced range: \R*\ <C -C a. In any case, at 
short distances, the physics of dilute systems is domi- 
nated by two-body effects. Under the conditions of di- 
luteness and low temperature discussed above, the solu- 
tion of the many-body problem with the full Hamilto- 
nian (|23[) is completely equivalent to the solution of an 
effective problem where the Hamiltonian only contains 
the kinetic energy term and the many-body wavefunc- 
tion satisfies the Bethe- Peierls boundary condition (j2"5|) . 

An approach which will be frequently considered in this 
review is based on microscopic simulations with quan- 
tum Monte Carlo techniques. In this case the contact 
boundary conditions ([25]) are difficult to implement and 
one must resort to a different effective interatomic po- 
tential. A convenient choice is the attractive square-well 
potential with range Rq and depth Vq defined in Eq. l|T8j) 
[other forms have also been considered in the literature 
(Carlson et ai, 2003)]. The interaction range Ro must 
be taken much smaller than the inverse Fermi wavevec- 
tor, kpRo -C 1, in order to ensure that the many-body 
properties of the system do not depend on its value. The 
depth Vb is instead varied so as to reproduce the actual 
value of the scattering length. 

The above approaches permit to describe the many- 
body features uniquely in terms of the scattering length 
a. These schemes are adequate if one can neglect the 
term in k 2 in the denominator of Eq. (|13p . When the 
effective range \R*\ of the interatomic potential becomes 
of the order of the inverse Fermi wavevector, as happens 
in the case of narrow Feshbach resonances, more complex 
effective potentials should be introduced in the solution 
of the many-body problem (see e.g. Gurarie and Radzi- 
hovsky, 2007). 



B. Order parameter, gap and speed of sound 

The phenomenon of superfluidity in 3D Fermi systems 
is associated with the occurrence of off-diagonal long- 
range order (ODLRO) according to the asymptotic be- 
havior (Gorkov, 1958) 

lim (*{(r 2 + r)*{(n+r)*;(ri)* T (r 2 )) = |F(n,r 2 )| 2 , 

(26) 

exhibited by the two-body density matrix. By assuming 
spontaneous breaking of gauge symmetry one can intro- 
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duce the pairing field 



F(R,s) = (^(R + s/2)^ T (R-s/2)) 



(27) 



[notice that ODLRO can be defined through Eq. |26|) 
also without using the symmetry breaking point of view 
(Yang, 1962)]. The vectors R = (n + r 2 )/2 and s = 
ri — Y2 denote, respectively, the center of mass and the 
relative coordinate of the pair of particles. In a Fermi 
superfluid ODLRO involves the expectation value of the 
product of two field operators instead of a single field 
operator as in the case of Bose-Einstein condensation. 
The pairing field in Eq. (f2"T|) refers to spin-singlet pairing 
and the spatial function F must satisfy the even-parity 
symmetry requirement F(R, — s) = F(R, s), imposed by 
the anticommutation rule of the field operators. 

The use of the Bethe-Peierls boundary conditions (123)) 
for the many-body wavefunction W implies that the pair- 
ing field (|27p is proportional to (1/s — 1/a) for small 
values of the relative coordinate s. One can then write 
the following short-range expansion 



F(R,s) = 



, ...A(R)(i-i 

Anh 2 Is a 



+ o(s), (28) 



which defines the quantity A(R), hereafter called the or- 
der parameter. The above asymptotic behavior holds 
in the same range of length scales as the Bethe-Peierls 
conditions (|25[) . For uniform systems the dependence on 
the center of mass coordinate R in Eqs. I|27 p -(|28 p drops 
and the order parameter A becomes constant. Moreover, 
in the case of s-wave pairing, the function (|27[) becomes 
spherically symmetric: F — F(s). The pairing field F(s) 
can be interpreted as the wavefunction of the macroscop- 
ically occupied two-particle state. The condensate frac- 
tion of pairs is then defined according to 



1 



^cond 



n/2 



ds \F{s) 



(29) 



where J rfs|F(s)| 2 is the density of condensed pairs. The 
quantity n CO nd is exponentially small for large, weakly 
bound Cooper pairs. It instead approaches the value 
n con d ~ 1 for small, tightly bound dimers that are al- 
most fully Bose-Einstein condensed at T = 0. 

Another peculiar feature characterizing superfluidity 
in a Fermi gas is the occurrence of a gap A gap in the 
single-particle excitation spectrum. At T = this gap 
is related to the minimum energy required to add (re- 
move) one particle starting from an unpolarized sys- 
tem according to the relationship E(N/2 ± l,N/2) = 
E(N/2, N/2) ±n + A gap . Here, £(iV T , N t ) is the ground- 
state energy of the system with iVjm particles of f (I) 
spin and n is the chemical potential defined by /i = 
±[E(N/2±l,N/2±l)-E(N/2,N/2)]/2 = dE/dN. By 
combining these two relations, one obtains the following 
expression for the gap (see, e.g., Ring and Schuck, 1980) 



A„ 



1 



[2E(N/2±l,N/2) 



(30) 



The gap corresponds to one half of the energy required 
to break a pair. The single-particle excitation spectrum, 
ejfcj is instead defined according to the relation Ek(N/2± 
1, N/2) = E(N/2,N/2)±ti+e k , where E k (N/2±l, N/2) 
denotes the energy of the system with one more (less) 
particle with momentum hk. Since A gap corresponds to 
the lowest of such energies E k , it coincides with the min- 
imum of the excitation spectrum. In general the order 
parameter A and the gap A gap are independent quanti- 
ties. A direct relationship holds in the weakly attractive 
BCS regime (see Sec. IIVDP where one finds A gap = A. 
The role of the gap in characterizing the superfluid be- 
havior will be discussed in Sec. IVII Dl 

Let us finally recall that a peculiar property of neu- 
tral Fermi superfluids is the occurrence of gapless den- 
sity oscillations. These are the Goldstone sound modes 
associated to the gauge symmetry breaking and are often 
referred to as the Bogoliubov- Anderson modes (Bogoli- 
ubov, Tolmachev and Shirkov, 1958; Anderson, 1958). 
These modes are collective excitations and should not be 
confused with the gapped single-particle excitations dis- 
cussed above. At small wavevectors they take the form 
of phonons propagating at T = with the velocity 



mc 2 = nd/i/dn 



(31) 



fixed by the compressibility of the gas. The description of 
these density oscillations will be presented in Sec. IVII Dl 
The Bogoliubov- Anderson phonons are the only gapless 
excitations in the system and provide the main contri- 
bution to the temperature dependence of all thermody- 
namic functions at low temperature (ksT <C A gap ). For 
example the specific heat C and the entropy S of the gas 
follow the free-phonon law C = 3S oc T 3 . 



C. Repulsive gas 

There are important cases where the many-body prob- 
lem for the interacting Fermi gas can be solved in an exact 
way. A first example is the dilute repulsive gas. Interac- 
tions are treated by means of the pseudo-potential ([15]) 
with a positive scattering length a. Standard pertur- 
bation theory can be applied with the small parameter 
kpa <C 1 expressing the diluteness condition of the gas. 
At T = 0, the expansion of the energy per particle up 
to quadratic terms in the dimcnsionless parameter fcjra 
yields the following expression (Huang and Yang, 1957; 
Lee and Yang, 1957) 



— — —Et? 
N 5 



1 



— kpa 
9vr 



4(11 -2 log 2) 



{k F af 



E(N/2 ± 1, N/2 ± 1) - E(N/2, N/2)] . 



21tt 2 

'(32) 

in terms of the Fermi energy (|24| . The above result 
is universal as it holds for any interatomic potential 
with a sufficiently small effective range \R*\ such that 
n|i?*| 3 <§; 1. Higher order terms in (f32|) will depend not 
only on the scattering length a, but also on the details of 
the potential (see Fetter and Walecka, 2003). In the case 
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of purely repulsive potentials, such as the hard-sphere 
model, the expansion (j52")) corresponds to the energy of 
the "true" ground state of the system (Lieb, Seiringer and 
Solovej, 2005; Seiringer, 2006). For more realistic poten- 
tials with an attractive tail, the above result describes 
instead the metastable gas-like state of repulsive atoms. 
This distinction is particularly important in the presence 
of bound states at the two-body level, since more stable 
many-body configurations satisfying the same condition 
k F a -C 1 consist of a gas of dimers (see Sec. IIV Ep . The 
weakly repulsive gas remains normal down to extremely 
low temperatures when the repulsive potential produces 
£ > pairing effects bringing the system into a super- 
fluid phase (Kohn and Luttinger, 1965; Fay and Layzer, 
1968; Kagan and Chubukov, 1988). In the normal phase 
the thermodynamic properties of the weakly repulsive gas 
are described with good accuracy by the ideal Fermi gas 
model. 



D. Weakly attractive gas 

A second important case is the dilute Fermi gas inter- 
acting with negative scattering length (A^M « 1). In 
this limit the many-body problem can be solved both at 
T = and at finite temperature and corresponds to the 
most famous BCS picture first introduced to describe the 
phenomenon of superconductivity (Bardeen, Cooper and 
Schrieffer, 1957). The main physical feature is the in- 
stability of the Fermi sphere in the presence of even an 
extremely weak attraction and the formation of bound 
states, the Cooper pairs, with exponentially small bind- 
ing energy. The many-body solution proceeds through 
a proper diagonalization of the Hamiltonian (|23[) by ap- 
plying the Bogoliubov transformation to the Fermi field 
operators (Bogoliubov, 1958). This approach is non per- 
turbative and predicts a second order phase transition 
associated with the occurrence of ODLRO. 

Exact results are available for the critical temperature 
and the superfluid gap (Gorkov and Melik-Barkhudarov, 
1961). For the critical temperature the result is 

T c = (-) lr F e 7 ' /2kFa « 0.28T F e 7r/2kFa , (33) 

where 7 = e c ~ 1.781, C being the Euler constant. The 
exponential, non analytical dependence of T c on the in- 
teraction strength k F \a\ is typical of the BCS regime. 
With respect to the original treatment by Bardeen, 
Cooper and Schrieffer (1957), the preexponential term 
in Eq. (|3"3"|) is a factor ~2 smaller as it accounts for the 
renormalization of the scattering length due to screening 
effects in the medium. A simple derivation of this result 
can be found in the book by Pethick and Smith (2002). 

The spectrum of single-particle excitations close to the 
Fermi surface, \k — k F \ <ti k F , is given by 

e k = ^A2 ap + [hv F (k - fc F )P , (34) 



where v F = hk F /m is the Fermi velocity, and is minimum 
at k — k F . The gap at T = is related to T c through 
the expression 

A gap = -k B T c a 1.76fc s T c . (35) 

7 

Furthermore, the ground-state energy per particle takes 
the form 

TP TP 3 A 2 

^ ^normal ° gap 

N = ~~ A ~8Ef~ ' [6b) 

where E norma i is the perturbation expansion (|32[) with 
a < and the term proportional to Ag Qp corresponds 
to the exponentially small energy gain of the superfluid 
compared to the normal state. 

Since the transition temperature T c becomes exponen- 
tially small as one decreases the value of fc,p|a|, the ob- 
servability of superfluid phenomena is a difficult task 
in dilute gases. Actually, in the experimentally rele- 
vant case of harmonically trapped configurations the pre- 
dicted value for the critical temperature easily becomes 
smaller than the typical values of the oscillator temper- 
ature huiho/ks- 

The thermodynamic properties of the BCS gas can 
also be investigated. At the lowest temperatures k B T <C 
Agap they are dominated by the Bogoliubov-Anderson 
phonons. However, already at temperatures k B T ~ A gap 
the main contribution to thermodynamics comes from 
fermionic excitations. For more details on the thermody- 
namic behavior of a BCS gas see, for example, the book 
by Lifshitz and Pitaevskii (1980). At T > T c the gas, 
due to the small value of T c , is still strongly degenerate 
and the thermodynamic functions are well described by 
the ideal-gas model. 



E. Gas of composite bosons 

Thanks to the availability of Feshbach resonances it 
is possible to tune the value of the scattering length in 
a highly controlled way (see Sec. IIII Bp . For example, 
starting from a negative and small value of a it is pos- 
sible to increase a, reach the resonance where the scat- 
tering length diverges and explore the other side of the 
resonance where a becomes positive and eventually small. 
One would naively expect to reach in this way the regime 
of the repulsive gas discussed in Sec. IIV Cl This is not the 
case because in the presence of a Feshbach resonance the 
positive value of the scattering length is associated with 
the emergence of a bound state in the two-body problem 
and the formation of dimers as discussed in Sec. IIII CI 
The size of the dimers is of the order of the scatter- 
ing length a and their binding energy is ef, ~ —h 2 /ma 2 . 
These dimers have a bosonic nature, being composed of 
two fermions, and if the gas is sufficiently dilute and cold 
they consequently give rise to the phenomenon of Bose- 
Einstein condensation. The size of the dimers can not 
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however be too small, as it should remain large com- 
pared to the size of the deeply-bound energy levels of the 
molecule. This requires the condition a 3> \R*\ which, 
according to the results of Sec. IIII C[ ensures that the 
system of weakly-bound dimers is stable enough and that 
the transition to deeper molecular states, due to collisions 
between dimers, can be neglected. 

The gas of dimers and the repulsive gas of atoms dis- 
cussed in Sec. IIV CI represent two different branches of 
the many-body problem, both corresponding to positive 
values of the scattering length (Pricoupenko and Castin, 
2004). The atomic repulsive gas configuration has been 
experimentally achieved by ramping up adiabatically the 
value of the scattering length, starting from the value 
a = (Bourdcl et al., 2003). If one stays sufficiently 
away from the resonance, losses are not dramatic and 
the many-body state is a repulsive Fermi gas. Conversely, 
the gas of dimers is realized by crossing adiabatically the 
Feshbach resonance starting from negative values of a, 
which allows for a full conversion of pairs of atoms into 
molecules, or by cooling down a gas with a fixed (posi- 
tive) value of the scattering length. 

The behavior of the dilute (kpa <C 1) gas of dimers, 
hereafter called BEC limit, is properly described by 
the theory of Bose-Einstein condensed gases available 
both for uniform and harmonically trapped configura- 
tions (Dalfovo et al, 1999). In particular, one can im- 
mediately evaluate the critical temperature T c . In the 
uniform case this is given by the text-book relationship 
T c = (2nh 2 /k B m)(n d /((3/2)) 2 / 3 , where n d is the density 
of dimers (equal to the density of each spin species) and 
£(3/2) ~ 2.612. In terms of the Fermi temperature (|24]) 
one can write 



T c = 0.218T F , 



(37) 



showing that the superfluid transition, associated with 
the Bose-Einstein condensation of dimers, takes place 
at temperatures of the order of Tp, i.e. at tempera- 
tures much higher than the exponentially small value 
(|3"3")) characterizing the BCS regime. The chemical po- 
tential of dimers, fid, is defined through the relationship 
2/i = eb + fJ-d involving the molecular binding energy e& 
and the atomic chemical potential \i. 

The inclusion of interactions between molecules, fixed 
by the dimer-dimer scattering lenght add according to 
the relationship add = 0.60a [see Eq. (|2"T1) ]. is provided 
to lowest order in the gas parameter rida dd by the Bo- 
goliubov theory for bosons with mass 2m and density 
rid = n a . At T = 0, the bosonic chemical potential is 
given by fid = 2irh 2 addnd/m. Higher order corrections 
to the equation of state are provided by the Lee-Huang- 
Yang expansion (Lee, Huang and Yang, 1957) 



E 
N 



£fc kFO-dd 

2 6tt 
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? (k F add) 3/2 



E F , (38) 



with small and positive scattering lengths, was proven by 
Leyronas and Combescot (2007). At very low tempera- 
tures the thermodynamics of the gas can be calculated 
using the Bogoliubov gapless spectrum ed{k) of density 
excitations 



e d (k) = hk (c| + h 2 k 2 /16m 2 ) 



>V2 



(39) 



where cb = y irfi 2 addnd/m 2 is the speed of Bogoliubov 
sound. The single-particle excitation spectrum is instead 
gapped and has a minimum at k — 0: e^o — A gap + 
k 2 /2m where 

A gap = TT + (3ew - a d d) • (40 

2 771 

In the above equation, the large binding-energy term 
|e&|/2 = h 2 /2ma 2 is corrected by a term which depends 
on both the dimer-dimer (add) and the atom-dimer (a a d) 
scattering length. This term can be derived from the 
definition (f3T)]) using for E(N^,Ni) an energy functional 
where the interactions between unpaired particles and 
dimers are properly treated at the mean-field level, the 
coupling constant being fixed by a ad and by the atom- 
dimer reduced mass. Since a a d — 1.18a [see Eq. (|2"0")) ]. the 
most important contribution to Eq. (|40p comes from the 
term proportional to the atom-dimer scattering length. 

In the BEC limit the internal structure of dimers can 
be ignored for temperatures higher than the critical tem- 
perature since in this regime one has 



\e b \/k B T c ~ l/(fc F a) 2 > 1 



(41) 



here expressed in units of the Fermi energy (|24p . The 
validity of the expression (|38p for a Fermi gas interacting 



The above condition ensures that at thermodynamic 
equilibrium the number of free atoms is negligible, be- 
ing proportional to e ei> / 2fcsT . 



F. Gas at unitarity 

A more difficult problem concerns the behavior of the 
many-body system when k F \a\ > 1, i.e. when the scat- 
tering length becomes larger than the interparticle dis- 
tance, which in turns is much larger than the range of the 
interatomic potential. This corresponds to the unusual 
situation of a gas which is dilute and strongly interact- 
ing at the same time. In this condition, it is not obvious 
whether the system is stable or collapses. Moreover, if the 
gas remains stable, does it exhibit superfluidity as in the 
BCS and BEC regimes? Since at present an exact solu- 
tion of the many-body problem for fcp|a| > 1 is not avail- 
able, one has to resort to approximate schemes or numer- 
ical simulations (see Sees IV AllVB[) . These approaches, 
together with experimental results, give a clear indica- 
tion that the gas is indeed stable and that it is superfluid 
below a critical temperature. The limit kp\a\ — > oo is 
called the unitary regime and has been already intro- 
duced in Sec. IIII Al when discussing two-body collisions. 
This regime is characterized by the universal behavior of 
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the scattering amplitude fo{k) — i/k which bears impor- 
tant consequences at the many-body level. As the scat- 
tering length drops out of the problem, the only relevant 
lcnght scales remain the inverse of the Fermi wavevector 
and the thermal wavelength. All thermodynamic quanti- 
ties should therefore be universal functions of the Fermi 
energy Ep and of the ratio T/Tp. 

An important example of this universal behavior is pro- 
vided by the T = value of the chemical potential: 

M = (1 + (3)E F , (42) 

where (3 is a dimensionless parameter. This relation 
fixes the density dependence of the equation of state, 
with non-trivial consequences on the density profiles and 
on the collective frequencies of harmonically trapped 
superfluids, as we will discuss in Sees. IVHIVIII The 
value of j3 in Eq. (j4"2"|) has been calculated using fixed- 
node quantum Monte-Carlo techniques giving the result 
(3 = -0.58 ±0.01 (Carlson et al, 2003; Astrakharchik et 
al., 2004a; Carlson and Reddy, 2005;). The most recent 
experimental determinations are in good agreement with 
this value (see Table [TT1 in Sec. IVII) . The negative value 
of [3 implies that at unitarity interactions are attractive. 
By integrating Eq. (|42p one finds that the same propor- 
tionality coefficient (1 + (3) also relates the ground-state 
energy per particle E/N and the pressure P to the cor- 
responding ideal gas values: E/N = (1 + f3)3Ep/5 and 
P = (1 + (3)2nEp /5, respectively. As a consequence the 
speed of sound (|3i"j) is given by 

c= (l+/3) 1/2 ^/V3 , (43) 

where vp/^/3 is the ideal Fermi gas value. The super- 
fluid gap at T = should also scale with the Fermi 
energy. Fixed-node quantum Monte Carlo simulations 
yield the result A gap = (0.50 ± 0.03)E F (Carlson et al, 
2003; Carlson and Reddy, 2005). A more recent QMC 
study based on lattice calculations (Juillet, 2007) gives 
for (3 a result consistent with the one reported above and 
a slightly smaller value for A gap (see also Carlson and 
Reddy, 2007). 

At finite temperature the most relevant problem, both 
theoretically and experimentally, is the determination of 
the transition temperature T c which is expected to de- 
pend on density through the Fermi temperature 

T c = aT F , (44) 

a being a dimensionless universal parameter. Quan- 
tum Monte-Carlo methods have been recently used to 
determine the value of a. Bulgac, Drut and Magierski 
(2006) and Burovski et al. (2006a) carried out simula- 
tions of fermions on a lattice where the sign problem, 
typical of fermionic quantum Monte-Carlo methods, can 
be avoided. These in principle "exact" studies require an 
extrapolation to zero filling factor in order to simulate 
correctly the continuum system and the reported value of 
the critical temperature corresponds to a — 0.157±0.007 



(Burovski et al., 2006a). Path integral Monte-Carlo sim- 
ulations, which work directly in the continuum, have also 
been performed by Akkinen, Ceperley and Trivedi (2006) 
using the restricted path approximation to overcome the 
sign problem. The reported value a ~ 0.25 is signifi- 
cantly higher compared to the previous method. 

Since at unitarity the gas is strongly correlated, one 
expects a significantly large critical region near T c . Fur- 
thermore, the phase transition should belong to the same 
universality class, corresponding to a complex order pa- 
rameter, as the one in bosonic liquid 4 He and should 
exhibit similar features including the characteristic A sin- 
gularity of the specific heat. 

The temperature dependence of the thermodynamic 
quantities is expected to involve universal functions of the 
ratio T/Tp (Ho, 2004). For example, the pressure of the 
gas can be written as P(n, T) — PT=o( n )fp(T/Tp) where 
Pt=o is the pressure at T — and fp is a dimension- 
less function. Analogously the entropy per atom takes 
the form S/Nkp — fg{T/Tp) involving the universal 
function fg related to fp by the thermodynamic relation 
dfp(x)/dx = xdfg(x)/dx. The above results for pressure 
and entropy imply, in particular, that during adiabatic 
changes the ratio T/Tp remains constant. This implies 
that the adiabatic processes, at unitarity, follow the law 
p n -5/3 _ cons t ; typical of non-interacting atomic gases. 

The scaling laws for the pressure and the entropy also 
hold at high temperatures, T ^> Tp, provided that the 
thermal wavelength is sti ll large c ompared to the effective 
range of interaction, h^J m/kpT \R*\. In this regime 
of temperatures the unitary gas can be described, to 
first approximation, by an ideal Maxwell-Boltzmann gas. 
Corrections to the equation of state can be determined 
by calculating the second virial coefficient B(T) defined 
from the expansion of the pressure P ~ nAflT[l+nB(T)] . 
Using the method of partial waves (Beth and Uhlenbeck, 
1937; Landau and Lifshitz, 1980) and accounting for the 
unitary contribution of the s-wave phase shift So = T 7r /2 
when a — > ±oo, one obtains the result 

It is worth noticing that the negative sign of the second 
virial coefficient corresponds again to attraction. Equa- 
tion (|45p takes into account effects of both Fermi statis- 
tics and interaction. The pure statistical contribution 
would be given by the same expression in brackets, but 
with the coefficient +1/4 of opposite sign. 

V. THE BCS-BEC CROSSOVER 

A. Mean- field approach at T — 

As discussed in Sec. IIVF1 there is not at present an 
exact analytic solution of the many-body problem along 
the whole BCS-BEC crossover. A useful approximation 
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is provided by the standard BCS mean-field theory of su- 
perconductivity. This approach was first introduced by 
Eagles (1969) and Leggett (1980) with the main motiva- 
tion to explore the properties of superconductivity and 
superfluidity beyond the weak-coupling limit kp\a\ <C 1. 
The main merit of this approach is that it provides a 
comprehensive, although approximate, description of the 
equation of state along the whole crossover, including 
the limit \jkpa — > and the BEC regime of small and 
positive a. At finite temperature the inclusion of fluc- 
tuations around the mean field is instead crucial to pro- 
vide a qualitatively correct description of the crossover 
(Nozieres and Schmitt-Rink, 1985; Sa de Melo, Randcria 
and Engelbrecht, 1993). In this Section we review the 
mean- field treatment of the crossover at T = 0, while 
some aspects of the theory at finite temperature will be 
discussed in Sec. IV CI 

The account of the BCS mean-field theory we give 
here is based on the use of the pseudo-potential (fT"5)) and 
follows the treatment of Bruun et al. (1999). Let us 
start considering a simplified Hamiltonian without ex- 
ternal confinement and where, in the interaction term 
of Eq. (|23p , only pairing correlations are considered and 
treated at the mean-field level 

#BCS=E«J**iW(-£-M)^(r) (46) 
-/dr{A(r) [^v)i>\{v)-\^\(v)^\{T))\+h.c] . 

We also restrict the discussion to equal masses m and 
to unpolarized systems: N\ = N\ = N/2. The direct 
(Hartree) interaction term proportional to the averages 
(*j(r)# T (r)) and (*t (r^fr)}, is neglected in Eq. gg) 
in order to avoid the presence of divergent terms in the 
theory when applied to the unitary limit X/a — > 0. The 
order parameter A is defined here as the spatial integral 
of the short-range potential V(s) weighted by the pairing 
field (123 



A(r) = - / dsF(s)(^(r + S /2)* T (r- S /2)) 

= -9(sF)' s=0 . (47) 

The last equality, which is obtained using the regular- 
ized potential (|15[) . is consistent with the definition of 
the order parameter given in Eq. (|28[) . The c- number 
term J drA(r)( 1 J r j(r)'^j(r))/2 + h.c. in the Hamiltonian 
(|46p avoids double counting in the ground-state energy, 
a typical feature of the mean-field approach. 

The Hamiltonian (|4"6"1) is diagonalized by the Bogoli- 

ubov transformation ^f(r) = J2i (ui( T )&i + v *{ r )Pi^j > 

^l( r ) — J2i ( w i( r )/3j — v *( r )oilj which transforms par- 
ticles into quasi-particles denoted by the operators dj 
and Pi. Since quasi-particles should also satisfy fermionic 
anti-commutation relations, {dj,d],} = = Sij', 

the functions Ui and Vi obey the orthogonality relation 
J dr[u*(r)iij(r) + v*(r)vj(r)] = Sij. As a consequence 



of the Bogoliubov transformations the Hamiltonian (|46]) 
can be written in the form 



which describes a system of independent quasi-particles. 
The corresponding expressions for the amplitudes Ui and 
Vi are obtained by solving the matrix equation 



Ho A(r) 
A*(r) -H 





(49) 



where Ho = — (h 2 /2m)V 2 — /i is the single-particle Hamil- 
tonian. The order parameter A(r) is in general a com- 
plex, position dependent function. Eq. (|4"9"| is known 
as the Bogoliubov-de Gennes equation (see de Gennes, 
1989). It can be used to describe both uniform and 
nonuniform configurations like, for example, quantized 
vortices (see Sec. I VIII CIO or solitons (Antezza, et al. 
2007). 

In the uniform case the solutions take the simple form 
of plane waves it, (r) 
with 



e lk - T u k /VV and v, -> e^ T v k /y/V 



l-vi 



1 



UkVk 



2e fc 



(50) 



i-2 . 2 

where rjk = %^ /i is the energy of a free particle 

calculated with respect to the chemical potential. The 
spectrum of quasi-particles (e, — * e^) which are the el- 
ementary excitations of the system, has the well-known 
form 



4 



(51) 



and close to the Fermi surface coincides with the result 
(|34p holding in the weak coupling limit with A gap = A. 
In this regime, the minimum of e k corresponds to the 
Fermi wavevector hp. The minimum is shifted towards 
smaller values of k as one approaches the unitary limit 
and corresponds to k = when the chemical potential 
changes sign on the BEC side of the resonance. Fur- 
thermore, one should notice that while the BCS the- 
ory correctly predicts the occurrence of a gap in the 
single-particle excitations, it is instead unable to describe 
the low-lying density oscillations of the gas (Bogoliubov- 
Anderson phonons). These can be accounted for by a 
time-dependent version of the theory (see, for example, 
Urban and Schuck, 2006). The vacuum of quasi-particles, 
denned by dk|0) = /3k|0) = 0, corresponds to the ground 
state of the system whose energy is given by 



h 2 k 2 

) 

2m 



A2 
2e fc 



(52) 



The above energy consists of the sum of two terms: the 
first is the kinetic energy of the two spin components, 
while the second corresponds to the interaction energy. 
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One should notice that both terms, if calculated sepa- 
rately, exhibit an ultraviolet divergence which disappears 
in the sum yielding a finite total energy. 

The order parameter A entering the above equations 
should satisfy a self-consistent condition determined by 
the short range behavior ([28]) of the pairing field ([27)) . 
This function takes the form 



F(s) 



dk 



-UkVke 



= A 



dk e iks 
(2tt) 3 2e k 



(53) 



By writing (Atts)^ 1 = J dke z]i ' s / [(2tt) 3 k 2 ] and comparing 
Eq. ([28]) with Eq. ([53| one straightforwardly obtains the 
important equation 



4irh 2 t 



dk f m 



(2tt) 3 \h 2 k 2 



1 

2~T k 



(54) 



where one is allowed to take the limit s — + since 
the integral of the difference in brackets is convergent. 
Eq. ([54"]) , through the expression (|5Tj) of the elementary 
excitations, provides a relationship between A and the 
chemical potential /i entering the single-particle energy 
7]k = — li. A second relation is given by the normal- 
ization condition 



V 



i 



Eft 



(55) 



which takes the form of an equation for the density. One 
can prove that the density dependence of the chemi- 
cal potential arising from the solution of Eqs. (|54[) and 
(|55p is consistent with the thermodynamic relation /i = 
dEo/dN, with E given by Eq. {52]). 

Result (|54|) can be equivalently derived starting from 
the contact potential g5(r) [rather than from the regu- 
larized form (|15p ] and using the renormalized value 



dk 



4irh 2 a J (2tt) 3 h 2 k 2 



(56) 



of the coupling constant, corresponding to the low-energy 
limit of the two-body T-matrix (Randeria, 1995). In this 
case one must introduce a cut-off in the calculation of the 
order parameter A = -g^^^} = -g J dkA/[2e fc (27r) 3 ], 
as well as in the integral in Eq. (|56p . in order to avoid 
the emergence of ultraviolet divergences. 

In the case of weakly attractive gases (hp \a\ <C 1 with 
a < 0) the chemical potential approaches the Fermi en- 
ergy ii ~ Ep and Eq. ([51)) reduces to the equation for 
the gap of standard BCS theory. In the general case 
the value of ii and A should be calculated by solving 
the coupled equations (|54|) and (|55|) . By expressing the 
energy in units of the Fermi energy Ep these equations 
only depend on the dimensionless parameter \jkpa which 
characterizes the interaction strength along the BCS- 
BEC crossover. In the following we will be referring to 
Eqs. ([54"| and {55]) as to the BCS mean-field equations. 

Analytical results for the energy per particle are ob- 
tained in the limiting cases l/kpa — » ±oo corresponding, 



respectively, to the BEC and BCS regimes. In the BCS 
limit the mean field equations give the result: 



N o 

while in the BEC limit one finds 

E _ h 2 3 fhkpa h{k F a) 

~N ~ ~~ 2ma 2 + 5 F \~9tt 54tt 2 



(57) 



(58) 



While the leading term in the energy per particle is cor- 
rectly reproduced in both limits (yielding, respectively, 
the non-interacting energy 3Ep/5 and half of the dimcr 
binding energy —h 2 /2ma 2 ), the higher order terms are 
wrongly predicted by this approach. In fact, in the BCS 
limit the theory misses the interaction-dependent terms 
in the expansion (f32|) . This is due to the absence of the 
Hartree term in the Hamiltonian (|4"6")l . In the BEC limit 
the theory correctly reproduces a repulsive gas of dimers. 
However, the term arising from the interaction between 
dimers corresponds, in the expansion (|58j) . to a molecule- 
molecule scattering length equal to add — 2a rather than 
to the correct value add = 0.60a (see Sec. lIIICj) . Further- 
more, the Lee-Huang- Yang correction in the equation of 
state of composite bosons [see Eq. (|38p ] is not accounted 
for by the expansion ([55]) . 

Finally, at unitarity (l/kpa — 0) one finds Eq/N ~ 
0.59(3-Ef/5) which is 40% larger than the value predicted 
by quantum Monte Carlo simulations (see Sec. IIVF]) . 

It is worth noticing that the energy per particle, as well 
as the chemical potential, change sign from the BCS to 
the BEC regime. This implies that there exists a value 
of kpa where fi = 0. This fact bears important con- 
sequences on the gap A gap characterizing the spectrum 
(fijTj) of single-particle excitations. If fx > 0, A gap coin- 
cides with the order parameter A. This is the case, in 
particular, of the BCS regime, where one finds the result 



A = -rE F exp 

e z 



2kpa 



(59) 



Notice that result ([59")) does not include the 
Gorkov— Melik-Barkhudarov [see Eqs. (|33p . 
([35)) ]. At unitarity one finds (Randeria, 1995) 



A, 



= A 



0.69E F . When n < 



instead given by A gap = y^A 2 



the gap is 
In particular 



in the BEC limit, where fi ~ —h 2 /2ma 2 , one finds 
A = (le/STr) 1 / 2 ^/^^ < \fj,\. In the same limit the 
gap is given by A gap = h 2 /2ma 2 + 3airh 2 n/m. 

The momentum distribution of either spin species — 
(a^akT) = (^kx^kj,) is another direct output of the BCS 
mean-field theory. For a given value of kpa along the 
crossover it is readily obtained from the corresponding 
values of fx and A through the expression 



1 



n k = v k = 



(60) 
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In the BCS regime, rik coincides approximately with 
the step function 8(1 — k/kp ) characteristic of the non- 
interacting gas, the order parameter A providing only a 
small broadening around the Fermi wavevector. By in- 
creasing the interaction strength the broadening of the 
Fermi surface becomes more and more significant. At 
unitarity, 1/kpa = 0, the effect is of the order of /cf, 
consistently with the size of A being proportional to the 
Fermi energy. In the BEC regime nk takes instead the 
limiting form 

4(fcFQ) 3 

" k = 3tt(1 + k 2 a 2 ) 2 ' (61) 

which is proportional to the square of the Fourier trans- 
form of the molecular wavefunction (|16p . 

It is important to notice that at large wavevec- 
tors the momentum distribution (|6"0)) decays as ~ 
m 2 A 2 /(h 4 k 4 ) for k > m\fi\/h. The large-fc 1/k 4 tail 
has important consequences for the kinetic energy of the 
system defined as E^i n — 2 n^ffk 2 /2m, which di- 
verges in 2D and 3D. This unphysical behavior arises 
from the use of the zero-range pseudopotential (|15p which 
describes correctly only the region of wavevectors much 
smaller than the inverse effective range of interactions, 
k -c l/|i?*|. It reflects the fact that the kinetic energy 
is a microscopic quantity that in general can not be ex- 
pressed in terms of the scattering length. 

Let us finally discuss the many-body structure of the 
ground-state wavefunction. The BCS ground state, de- 
fined as the vacuum state for the quasi-particles dk and 
/3k, can be written explicitly in terms of the amplitudes 
Uk, Vk giving the well known result 



|BCS) =Y[{u k + v k at ] al k ^j |0) , 



(62) 



where |0) is the particle vacuum. It is important to stress 
that the BCS mean-field Eqs. (j54"|) - (p35)) can be equiva- 
lently derived from a variational calculation applied to 
the state (|62[) where the grand-canonical energy is min- 
imized with respect to it/., v k , subject to the normaliza- 
tion constraint u\+v 2 = 1. The BCS state (|6"2"|) does not 
correspond to a definite number of particles. In fact, it 
can be decomposed into a series of states having an even 
number of particles |BCS) oc |0) + pt|0) + (P^) 2 \$) + 



where P t = Y^^ v k/ U k 



l kT U -ki 



is the pair creation op- 



erator. By projecting the state (f6"2"| onto the Hilbcrt 
space corresponding to N particles one can single out 
the component |BCS)jv oc (P i <) N / 2 \0) of the series. In 
coordinate space this iV-particle state can be expressed 
in terms of an antisymmetrized product of pair orbitals 
(Leggett, 1980) 

*bcs (ri, ...,r N ) = A [<M r ii') < X r 22')---<K r A r T jvJ] ■ (63) 

Here A is the antisymmetrizer operator and the func- 
tion <j)(r) = (2ir)~ 3 J dk(vk/uk)e l]si ' r depends on the rel- 
ative coordinate ru> = |r, — iy| of the pair of parti- 
cles, i (i') being labels for the spin-up (down) parti- 
cles. It is worth noticing that in the deep BEC regime, 



corresponding to ~ ti 2 /2ma 2 ^> A, the pair or- 
bital becomes proportional to the molecular wavefunction 
(HU): 4>{r) = (^Av/2)e- r /7(%/2W), and the many- 
body wavefunction (|63[) describes a system where all 
atoms are paired into bound molecules. Wave functions 
of the form (|63p are used in order to implement more 
microscopic approaches to the many-body problem (see 
Sec. EE). 

In conclusion, we have shown how BCS mean-field the- 
ory is capable of giving a comprehensive and qualitatively 
correct picture of the BCS-BEC crossover at T = 0. The 
quantitative inadequacies of the model will be discussed 
m more details in Sec. IVB21 



B. Quantum Monte Carlo approach at T = 

1. Method 

A more microscopic approach to the theoretical inves- 
tigation of the ground-state properties of the gas along 
the BCS-BEC crossover is provided by the fixed-node dif- 
fusion Monte Carlo (FN-DMC) technique. This method 
is based on the Hamiltonian (|2"3"|) where, as in Sec. IV A[ 
we consider uniform and unpolarized configurations of 
particles with equal masses. A convenient choice for the 
effective interatomic potential V e ff (r) consists of using 
the square-well model (I18|) where, in order to reduce 
finite-range effects, the value of Ro is taken as small 
as nRl = 10 -6 . The depth V = h 2 K 2 /m of the po- 
tential is varied in the range < KqRq < n to re- 
produce the relevant regimes along the crossover. For 
K R < 7r/2, the potential does not support a two-body 
bound state and the scattering length is negative. Vice- 
versa, for KqRq > 7t/2, the scattering length is positive 
and a molecular state appears with binding energy e&. 
The value KqRq = n/2 corresponds to the unitary limit, 
where \a\ = oo and = 0. 

In a diffusion Monte Carlo (DMC) simulation one in- 
troduces the function /(R, t) = '0t(R.) v E'(R, t), where 
^(R, r) denotes the wavefunction of the system and 
■0t(R) is a trial function used for importance sampling. 
The function /(R, r) is evolved in imaginary time, r = 
it/h, according to the Schrodinger equation 



df(R,r) 
dr 



£{V-h/CH.t) 
[£l(R) 



V r [F(R)/(R,t)]} 
E ref ]f(R,T). (64) 



In the above equation R denotes the position vec- 
tors of the N particles, R = (ri, rjv), El (TV) = 
^(R)" 1 HtJ)t(TV) denotes the local energy, F(R) = 
2t/>t(R)~ 1 Vri/'t(R) is the quantum drift force, and E re f 
is a reference energy introduced to stabilize the numer- 
ics. The various observables of the system are calculated 
from averages over the asymptotic distribution function 
/(R, t — > oo) (for more details see e.g. Boronat and Ca- 
sulleras, 1994). As an example, the DMC estimate of the 
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energy is obtained from 



Edmc = 



f dREx,(R)/(R, r — > oo) 
JdR/(R,r^ oo) 



(65) 



A crucial requirement, which allows for the solution of 
Eq. as a diffusion equation, is the positive definitc- 
ness of the probability distribution /(R, r). The con- 
dition /(R, t) > can be easily satisfied for the ground 
state of bosonic systems by choosing both <3> and ipT pos- 
itive definite, corresponding to the nodeless ground-state 
function. In this case the asymptotic distribution con- 
verges to /(R, r —> oo) = Vt(R)*o(R), where * (R) 
is the exact ground-state wavefunction, and the average 
(|65|) yields the exact ground-state energy Eq. The case of 
fermionic ground-states or of more general excited states 
is different, due to the appearance of nodes in the wave- 
function \&. In this case an exact solution is in general not 
available. An approximate treatment is provided by the 
FN-DMC method which enforces the positive definiteness 
of /(R, t) through the constraint that ipr and \& change 
sign together and thus share the same nodal surface. The 
nodal constraint is kept fixed during the calculation and 
the function /(R, r), after propagation in imaginary time 
according to Eq. ([M]) , reaches for large times the asymp- 
totic distribution /(R, r — > oo) = ?/'t(R)^ , f7v(R), where 
^fjv(R) is an approximation to the exact eigenfunction 
of the many-body Schrodinger equation. It can be proved 
that, due to the nodal constraint, the fixed-node energy 
obtained from Eq. (|65p is an upper bound to the ex- 
act eigenenergy for a given symmetry (Reynolds et al., 
1982). In particular, if the nodal surface of ipr were 
exact, then £dmc would also be exact. The energy cal- 
culated in a FN-DMC simulation depends crucially on a 
good parametrization of the many-body nodal surface. 

The calculations are carried out using a cubic simu- 
lation box of volume V — L 3, with periodic boundary 
conditions. The number of particles in the system ranges 
typically from N = 14 to N = 66 and finite-size analysis 
are performed to extrapolate the results to the thermo- 
dynamic limit. The most general trial wavefunction used 
in studies of ultracold fermionic gases has the form (Carl- 
son et al., 2003; Chang et al., 2004; Astrakharchik et al., 
2004a; Chang and Pandharipande, 2005; Astrakharchik 
et al, 2005a) 



V^t(R) = *j(R)*bcs(R) 



(66) 



where contains Jastrow correlations between all the 
particles, *j(R) = Hi j faa'(\r ia - r ja >), where r ia de- 
notes the position of the i-th particle with spin a, and 
the BCS-type wavefunction ^bcs is an antisymmetrized 
product of pair wavefunctions of the form (|63[) . The pair 
orbital </>(r) is chosen of the general form 



(r) 



< £ 

k a <km, 



Mr) 



(67) 



where a is a variational parameter and the sum is per- 
formed over the plane-wave states satisfying periodic 



boundary conditions, k Q = 2ir/L(£ ax x + £ ay y + t az z) 
(the ts are integer numbers), up to the largest closed 
shell k max = 2i:/L{t 2 maxx + i 2 naxy + l 2 maxz ) 1/2 occupied 
by the N/2 particles. A convenient functional form of the 
Jastrow correlation terms f aa i (r) and of the s-wave or- 
bital 4>s (r) is discussed by Astrakharchik et al. (2005a). 
The Jastrow function in Eq. (|66[) is chosen positive defi- 
nite, \&j(R) > 0, and therefore the nodal surface of the 
trial function is determined only by ^bcs- 



An important point concerns the wavefunction 'J'bcs 
which can be used to describe both the normal and the 
superfluid state. In fact, if one chooses in Eq. (|6"T|) 
4> s (r) = 0, it can be shown (Bouchaud, Georges and 
Lhuillier, 1988; Bouchaud and Lhuillier, 1989) that <f B cs 
coincides with the wavefunction of a free Fermi gas, i.e., 
the product of the plane-wave Slater determinants for 
spin-up and spin-down particles. In this case the trial 
wavefunction (f66|) is incompatible with ODLRO and de- 
scribes a normal Fermi gas similarly to the wavefunction 
employed in the study of liquid 3 He at low temperatures 
(Ceperley, Chester and Kalos, 1977). On the contrary, 
if </> s (r) ^ 0, the wavefunction (|66[) accounts for s-wave 
pairing and describes a superfluid gas. In particular, in 
the deep BEC limit 1/kpa ^> 1, the choice a = and 
<j) s {r) given by the molecular solution of the two-body 
problem in Eq. (|6Tj) reproduces the mean-field wavefunc- 
tion (|63p . As an example, at unitarity, the normal-state 
wavefunction [<f>s{r) — in Eq. (|6"T)) ] yields the value 
E/N = 0.56(1)3£'f/5 for the energy per particle, to be 
compared with the result E/N = 0.42(l)3i?F/5 obtained 
with the superfluid-state wavefunction. 



Another important remark concerns the gas-like nature 
of the many-body state described by the wave function 
^i?Ar(R). In the limit of a zero-range interatomic po- 
tential this state corresponds to the true ground-state of 
the system, because bound-states with more than two 
particles are inhibited by the Pauli exclusion principle 
(Baker, 1999). The situation is different for a finite- 
range potential. In the case of the square-well poten- 
tial (fTS)) one can easily calculate an upper bound to the 
ground-state energy using the Hartree-Fock state |HF) = 
U k <k F 4,t4j°>- One finds the result (HP|F|HF>/JV = 
3Ep/5 — ttVouRq/S, showing that at large n the inter- 
action energy is unbounded from below and the kinetic 
energy can not prevent the system from collapsing. How- 
ever, for realistic short-range potentials having a large 
depth Vq ~ h 2 /itiRq, this instability sets in at very high 
densities on order nR^ ~ 1. Such large density fluc- 
tuations are extremely unlikely, so that one can safely 
ignore this collapsed state when performing the simula- 
tions. Thus, for small enough values of i?o, the gas-like 
state corresponding to the wavefunction ^fn is expected 
to describe the ground state of the system with zero- 
range interactions. 
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2. Results 

The results of the FN-DMC calculations are reported 
in Figs. OH In Fig. [5] we show the energy per particle 
as a function of the interaction strength along the BCS- 
BEC crossover (Astrakharchik et at, 2004a; Chang et al., 
2004). In order to emphasize many-body effects on the 
BEC side of the resonance we subtract from E/N the 
two-body contribution e&/2 arising from the molecules. 
Notice that refers here to the dimer binding energy 
in the square- well potential (fT8|) . which, for the largest 
values of l/kpa, includes appreciable finite-range correc- 
tions compared to —h 2 /ma 2 . Nevertheless, no appre- 
ciable change is found in the difference E/N — e&/2 if 
the value of Rq is varied, demonstrating the irrelevance 
of this length scale for the many-body problem. Both 
in the BCS and in the BEC regime the FN-DMC ener- 
gies agree respectively with the perturbation expansion 
(|3"2")l and In particular, in the inset of Fig. \5\ we 

show an enlarged view of the results in the BEC regime 
indicating good agreement with the add — 0.60a result 
for the dimer-dimer scattering length as well as evidence 
for beyond mean-field effects. In Fig. [5] we compare the 
FN-DMC results with the ones from the BCS mean-field 
theory of Sec. IV Al The mean-field results reproduce the 
correct qualitative behavior, but are affected by impor- 
tant quantitative inadequacies. 

The quantum Monte Carlo and mean-field results for 
the momentum distribution and for the condensate 
fraction of pairs n con d are reported in Fig. [S] and re- 
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FIG. 5: (color online). Energy per particle along the BCS- 
BEC crossover with the binding-energy term subtracted from 
E/N. Symbols: quantum Monte Carlo results from As- 
trakharchik et al. (2004a). The dot-dashed line corresponds 
to the expansion (I32[l and the dashed line to the expansion 
(|38|) holding, respectively, in the BCS and in the BEC regime. 
The long-dashed line (red) refers to the result of the BCS 
mean-field theory. Inset: enlarged view of the BEC regime 
— 1/kFa < — 1. The solid line corresponds to the mean-field 
term in the expansion (138|) . the dashed line includes the Lee- 
Huang- Yang correction. 



spectively (Astrakharchik et al., 2005a). In both cases 
the mean-field predictions are in reasonable agreement 
with the findings of FN-DMC calculations. Significant 
discrepancies are found in the momentum distribution 
at unitarity (FigJH]) (see also Carlson et al., 2003), where 
the broadening of the distribution is overestimated by the 
mean-field theory consistently with the larger value pre- 
dicted for the pairing gap. The momentum distribution 
in harmonic traps is discussed in Sec. I VI CI 

An important remark concerns the condensate frac- 
tion n con d defined in Eq. (129|) which, in the BEC regime, 
should coincide with the Bogoliubov quantum depletion 
"cond = 1 — 8v n da d d/3\/Tr characterizing a gas of in- 
teracting composite bosons with density rid = n/2 and 
scattering length add = 0.60a. This behavior is indeed 
demonstrated by the FN-DMC results (Fig|7J), but is 
not recovered within the mean-field approximation. At 
— l/fci?a ~ 1 on the opposite side of the resonance, the 
FN-DMC results agree with the condensate fraction cal- 
culated using BCS theory including the Gorkov— Mclik- 
Barkhudarov correction. This result is expected to repro- 
duce the correct behavior of n con d in the deep BCS regime 
and is significantly lower as compared to the mean-field 
prediction. 

Condensation of pairs has been observed on both sides 
of the Feshbach resonance by detecting the emergence 
of a bimodal distribution in the released cloud after the 
conversion of pairs of atoms into tightly bound molecules 
using a fast magnetic-field ramp (Regal, Greiner and Jin, 
2004b; Zwierlein et al, 2004; Zwierlein et al, 2005a). 
The magnetic-field sweep is slow enough to ensure full 
transfer of atomic pairs into dimers, but fast enough to 
act as a sudden projection of the many-body wavefunc- 
tion onto the state of the gas far on the BEC side of the 
resonance. The resulting condensate fraction is an out- 
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FIG. 6: (color online). QMC results of the momentum dis- 
tribution rik for different values of 1/kpa (solid lines). The 
dashed lines (red) correspond to the mean-field results from 
Eq. (ESI). Inset: n k for l/k F a = 4. The dotted line (blue) 
corresponds to the momentum distribution of the molecular 
state Eq. (|6Tj) . 
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of-equilibrium property, whose theoretical interpretation 
is not straightforward (Perali, Pieri and Strinati, 2005; 
Altman and Vishwanath, 2005), but it strongly supports 
the existence of ODLRO in the gas at equilibrium also 
on the BCS side of the Feshbach resonance (negative a) 
where no stable molecules exist in vacuum. 

Another quantity that can be calculated using the FN- 
DMC method is the spin-dependent pair correlation func- 
tion (Astrakharchik et al., 2004a; Chang and Pandhari- 
pande, 2005) defined as 

W (|s - s'|) = ±($1(b')&M$„(b)$A*')) ■ (68) 

This function measures the relative probability of finding 
a particle with equal or opposite spin at distance |s — s'\ 
from a given particle. The results for g-\-\{r) are shown 
in Fig. [5] for different values of the interaction strength. 
Notice that in all cases g-\-\(r) must vanish at short dis- 
tances due to the Pauli exclusion principle. This ten- 
dency of fermions to avoid each other (anti-bunching) , 
as opposed to the bunching effect exhibited by bosons, 
has been recently revealed in the experiment by Jeltes 
et al. (2007) using helium atoms. In the BCS regime, 
where the effects of interaction in the t _ T spin chan- 
nel are negligible, one expects that the pair correlation 
function is well described by the non-interacting gas re- 
sult <7tT( r ) = 1 — 9/(kpr) 4 [sin(kpr)/kpr — cos(fcpr)] 2 . 
Quite unexpectedly this result holds true even at uni- 
tarity. Only when one approaches the BEC regime the 
effect of indirect coupling, mediated through interactions 
with the opposite spin component, becomes relevant and 
5tt( r ) deviates significantly from the behavior of the non- 
interacting gas. In particular, for the largest value of 
1/kpa reported in Fig. [8] (1/kpa = 4), we show the 
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FIG. 7: (color online). Condensate fraction of pairs n con( j 
[Eq. (f29|) ] as a function of the interaction strength: FN-DMC 
results (symbols), Bogoliubov quantum depletion of a Bose 
gas with add = 0.60a [dashed line (red)], BCS theory includ- 
ing the Gorkov— Melik-Barkhudarov correction [dot-dashed 
line (blue)] and mean-field theory using Eq. (|53[) [solid line 
(green)]. 



pair distribution function of a gas of weakly interacting 
bosons of mass 2m, density n/2 and scattering length 
a>dd = 0.60a calculated within Bogoliubov theory. For 
large distances, r 3> add, where the Bogoliubov approxi- 
mation is expected to hold, one finds a remarkable agree- 
ment. 

Finally, in Fig. [5J we report the results for the pair 
correlation function Qn(r). In the f-J, spin channel in- 
teractions are always relevant and give rise to a 1/r 2 di- 
vergent behavior at short distances, the coefficient be- 
ing determined by many-body effects on the BCS side 
of the resonance and at unitarity and by the molecular 
wavefunction in the deep BEC regime. In the latter case 
one finds (kpr) 2 g-\i(r) — > Zir/{kpa) while at unitarity 
one finds (kpr) 2 gn(r) — > 2.7. The divergent behavior 
of the pair-correlation function at short distances gives 
rise to a sizable bunching effect due to interactions in 
the f-J, spin channel, as opposed to the anti-bunching ef- 
fect due to the Pauli principle in the t _ T spin channel 
(Lobo et al., 2006a). The function <?fx( r ) nas been the 
object of experimental studies using spectroscopic tech- 
niques (Greiner et al., 2004; Partridge et al., 2005). In 
particular, Partridge et al. (2005) measured the rate of 
molecular photoexcitation using an optical probe sensi- 
tive to short-range pair correlations. They found that 
the rate is proportional to 1/kpa in the BEC regime and 
decays exponentially in the BCS regime. 

Following the proposal by Altman, Demler and Lukin 
(2004), pair correlations have also been detected using 
the atom shot noise in absorption images (Greiner et at, 




FIG. 8: (color online). QMC results for the pair correlation 
function of parallel spins, gn(r), for different values of the 
interaction strength: solid line (red) 1/kFa = 0, dashed line 
(green) 1/kpa = 1, dotted line (blue) 1/kpa = 4. The thin 
solid line (black) refers to the non-interacting gas and the 
thin dotted line (black) is the pair correlation function of a 
Bose gas with add = 0.60a and the same density as the single- 
component gas corresponding to 1/kpa = 4. 
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2005). The noise, related to the fluctuations of the col- 
umn integrated density, is extracted from 2D absorption 
images of the atom cloud by subtracting from each image 
pixel the azimuthal average of the signal. Crucial is the 
size of the effective image pixel (~15^tm) which should 
be small enough to be sensitive to atom shot noise. Us- 
ing this technique spatial |-| pair correlations have been 
observed on the BEC side of the resonance by comparing 
pixel-to-pixel the processed noise images relative to the 
two spin components. These images are obtained imme- 
diately after dissociating the molecules through a rapid 
sweep of the scattering length across the resonance. Even 
more spectacular is the observation of nonlocal pair cor- 
relations between atoms that have equal but opposite 
momenta and are therefore detected at diametrically op- 
posite points of the atom cloud in time-of-flight expan- 
sion. These correlations in momentum space are pro- 
duced by dissociating the molecules and by allowing the 
gas to expand before imaging. An important requirement 
here, which for the moment has limited the application 
of this technique to the BEC regime, is that the relative 
momentum should be significantly larger than the center 
of mass motion of the pairs, since this latter would de- 
grade the correlation signal due to image blurring. This 
method provides a novel, useful tool for detecting quan- 
tum correlations in many-body systems (Foiling et al, 
2005). 




FIG. 9: (color online). QMC results for the pair correlation 
function of antiparallel spins, gn(r), for different values of the 
interaction strength: dashed line (green) 1/kpa = — 1, solid 
line (red) 1/kpa = 0, dotted line (blue) 1/kpa = 4. The thin 
solid line (black) refers to the non-interacting gas. 



C. Other theoretical approaches at zero and finite 
temperature 

The extension of the BCS mean-field approach dis- 
cussed in Sec. IV Al to finite temperatures requires the in- 
clusion of thermal fluctuations in the formalism (see, e.g., 
Randeria, 1995). This can be accomplished by expanding 
the effective action determining the partition function of 
the system to quadratic terms in the order parameter 
A (Nozieres and Schmitt-Rink, 1985; Sa de Melo, Ran- 
deria and Engelbrecht, 1993). The method goes beyond 
the saddle-point approximation which corresponds to the 
mean- field theory at T = 0. The resulting predictions 
for T c are shown in Fig. O At unitarity one finds the 
value T c = 0.224Tp, which is not too far from the most 
reliable theoretical estimate based on quantum Monte 
Carlo simulations (Burovski et at, 2006a) discussed in 
Sec. IIVF1 In the BEC regime the region between the 
transition temperature T c and the higher temperature 
scale T* = ft 2 /kBiria 2 , fixed by the molecular binding 
energy, corresponds to the so called pseudogap regime, 
where bound pairs exist but are not "condensed" and 
form a normal phase (Chen et al., 2005). The presence 
of pseudogap effects near unitarity has been investigated 
by Stajic et al. (2004) and by Perali et al. (2004). 

Diagrammatic methods based on the T-matrix ap- 
proach have been proposed to extend the original treat- 
ment by Nozieres and Schmitt-Rink (1985) to the broken- 
symmetry phase below T c (Pieri, Pisani and Strinati, 
2004) as well as to improve on it by including higher 
order diagrams (Haussmann, 1993 and 1994; Chen et 
al., 2005; Combescot, Leyronas and Kagan, 2006; Hauss- 
mann et al, 2006). In particular, in the recent approach 
by Haussmann et al. (2006) based on a ladder approxi- 
mation, the fermionic degrees of freedom are accounted 
for using interacting Green's functions which are deter- 
mined in a self-consistent way. This approach applies 
to arbitrary temperatures and interaction strengths. At 
unitarity it predicts the value T c = 0.16Tp for the tran- 
sition temperature. 

Fully non-perturbative numerical techniques have also 
been applied to investigate the thermodynamics at fi- 
nite temperature in the unitary regime. They are based 
on the auxiliary field (Bulgac, Drut and Magierski, 2006 
and 2007) and diagrammatic determinant (Burovski et 
al., 2006a and 2006b) quantum Monte Carlo method 
on a lattice and on the restricted path-integral Monte 
Carlo method in the continuum (Akkineni, Ceperley and 
Trivedi, 2006). The results of these approaches for the 
critical temperature and the thermodynamic functions 
have already been discussed in Sec. IIVFI and the ones 
referring to harmonically trapped configurations are dis- 
cussed in Sec. lVI PI Lattice quantum Monte Carlo meth- 
ods have also been recently applied to investigate the 
ground-state properties at unitarity (Lee, 2006; Juillct, 
2007). 

An alternative approach to the theoretical treatment 
of the BCS-BEC crossover is provided by the two-channel 
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model (also called Bose-Fermi model). In this approach 
(Friedberg and Lee, 1989) the Hamiltonian includes in- 
teraction terms involving both fermionic and bosonic de- 
grees of freedom. A thorough account of the two-channel 
model can be found in the recent work by Gurarie and 
Radzihovsky (2007). This work discusses the compar- 
ison with the single-channel Hamiltonian (|2"5|) and the 
predictions of the model in the case of narrow Feshbach 
resonances, where the many-body problem can be exactly 
solved using a perturbative expansion (see also Dichl and 
Wetterich, 2006). In the case of broad resonances this ap- 
proach reduces to the single-channel Hamiltonian (f2"3"| . 
where interactions are accounted for by a contact po- 
tential fixed by the scattering length. The two-channel 
model is thus more general and can describe situations 
where the effective range plays an important role. The 
two-channel model was first introduced in the context of 
fcrmions with resonantly enhanced interactions by Hol- 
land et al. (2001) and by Timmermans et al. (2001) 
and was later developed to describe properties of the 
BCS-BEC crossover both at zero (Bruun and Pethick, 
2004; Romans and Stoof, 2006) and at finite temper- 
ature (Ohashi and Griffin, 2002 and 2003a; Diehl and 
Wetterich, 2006; Falco and Stoof, 2007) as well as in 
trapped configurations (Ohashi and Griffin, 2003b). 

Analytical methods based on an expansion around 
D = 4 — e spatial dimensions have also been recently 
applied to the unitary Fermi gas both at T — (Nishida 
and Son, 2006a and 2006b) and at finite temperature 
(Nishida, 2006). The starting point of the method is 
the observation (Nussinov and Nussinov, 2006) that a 
unitary Fermi gas in D = 4 behaves as an ideal Bose 
gas, i.e. that at T — the proportionality coefficient in 
Eq. (|4*2")l is 1 + /3 = 0. Results for the physical case of 
D = 3 are obtained by extrapolating the series expan- 
sion to e — 1. A similar approach is based on a l/N 
expansion, where N is the number of fermionic species 
(Nikolic and Sachdev, 2007; Veillette, Sheehy and Radzi- 
hovsky, 2007). For N — ► oo the field-theoretical problem 
can be solved exactly using the mean-field theory. Cor- 
rections to the mean-field predictions can be calculated 
in terms of the small parameter l/N and the results can 
be extrapolated to the relevant case of N = 2. 

Other theoretical approaches that have been applied 
to the physics of the BCS-BEC crossover include the 
dynamical mean-field theory (Garg, Krishnamurthy and 
Randeria, 2005) and the renormalization group method 
(Nikolic and Sachdev, 2007; Dichl et al, 2007) and the 
generalization of the BCS mean-field theory to include ef- 
fective mass and correlation terms within a density func- 
tional approach (Bulgac, 2007). 



VI. INTERACTING FERMI GAS IN 
HARMONIC TRAP 

The solution of the many-body problem for non- 
uniform configurations is a difficult task involving in 



most cases numerical calculations which are more com- 
plex than in uniform matter (an example of this type of 
numerical studies will be given when discussing the struc- 
ture of vortex configurations in Sec. IVIID|) . However, in 
the experimentally relevant case of quite large systems 
(N ~ 10 5 -10 7 ) confined by a harmonic potential, the lo- 
cal density approximation (LDA) provides a reliable and 
relatively simple description. This approximation, which 
is also often referred to as semiclassical or Thomas-Fermi 
approximation, profits of the knowledge of the equation 
of state of uniform matter to infer on the behavior of the 
system in traps. 

Let us also point out that on the BEC side of the res- 
onance, where the interacting Fermi gas behaves like a 
gas of weakly interacting dimers, systematic information 
is available from our advanced knowledge of the physics of 
dilute Bose gases in traps (see, e.g., Dalfovo et at, 1999). 
Although the deep BEC regime is not easily achieved 
in present experiments with ultracold Fermi gases, the 
corresponding predictions nevertheless provide a useful 
reference. 



A. Local density approximation at T — 0: density 
profiles 

In this Section we will consider systems at zero tem- 
perature where the equation of state of the uniform gas 
is provided by the density dependence e(n) of the energy 
density. The LDA consists of assuming that, locally, the 
system behaves like a uniform gas, so that its energy den- 
sity can be expressed as e(n) — nE(n)/N where E(n)/N 
is the energy per atom of uniform matter. The energy of 
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FIG. 10: Experimental results of the double integrated den- 
sity profiles along the BCS-BEC crossover for a gas of 6 Li 
atoms. The results at 850 G correspond to unitarity, while 
the ones at 809 G and 882 G correspond respectively to the 
BEC and BCS side of the resonance. The continuous curve at 
unitarity, is the best fit based on Eq. (J77J). The dashed lines 
correspond to the predictions for a non-interacting gas. From 
Bartenstein et al. (2004a). 
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the trapped system is then written in the integral form 

E = J dr{e[n{r)] + V ho (r)n(r)} , (69) 

and consists of the sum of the internal (also called release) 
energy 

E rei = J dr e [n(r)] , (70) 

and of the oscillator energy 

Eho = J drV ho (r)n(r) , (71) 

provided by the trapping potential Vho (r) , introduced in 
Eq. flj), which is assumed to be the same for both spin 
species. Furthermore, we also assume N-\ = N^. In 
Eqs. (|69I71|) . n(r) = n-f(r) + nj_(r) is the total density 
profile. Its value at equilibrium is determined by the 
variational relation 8(E — (XqN) / 'oYi(r) — 0, which yields 
the Thomas-Fermi equation 

Mo = /*[n(r)] + V ho (r) , (72) 

where //(n) = <9e(n) /9n is the local chemical potential de- 
termined by the equation of state of the uniform system 
and /io is the chemical potential of the trapped gas, fixed 
by the normalization condition J dr n(r) = N. Eq. Q72j) 
provides an implicit equation for n(r). 

The applicability of LDA is justified if the relevant en- 
ergies are much larger than the single-particle oscillator 
energy hu>i, i.e. if (j,q hiOi (i = x,y,z). While in the 
case of Bose-Einstein condensed gases this condition is 
ensured by the repulsive interaction among atoms, in the 
Fermi case the situation is different. In fact, thanks to 
the quantum pressure term related to the Pauli princi- 
ple, even in the non-interacting case one can apply the 
Thomas-Fermi relationship (|72p using the density depen- 
dence 

^) = M 2/ r" 2/3 . ( 73 ) 
2m 

for the local chemical potential yielding the equilibrium 
profile ©. 

Interactions modify the shape and the size of the den- 
sity profile. The effects are accounted for by Eq. (|72| once 
the local chemical potential n(n) is known. A simple re- 
sult is obtained at unitarity, where the equation of state 
has the same density dependence ([73")) as for the ideal 
gas, apart from the dimcnsionlcss rcnormalization factor 
(1+/3) [see Eq. (|32]l ]. By dividing Eq. (72]) by (1 + /3), one 
finds that the results at unitarity are obtained from the 
ones of the ideal Fermi gas by simply rescaling the trap- 
ping frequencies and the chemical potential according to 
uji^r ujij + (3 and fj,o —> Mo/(l + 0). In particular, the 
density profile at unitarity takes the same form ([5]) as in 
the ideal gas, the Thomas- Fermi radii being given by the 
rescaled law 

Ri = (l + /3) 1/4 i?° = (1 + /?) 1 / 4 a /lo (247V) 1 /6^ . (74) 

UJi 



From the above results one also finds the useful ex- 
pression E ho = (1 + Py^ 2 E^ for the oscillator energy 
(fTTjl of the trapped gas in terms of the ideal gas value 
E° ho = (3/8)NEg° (see Sec. EH]). 

In the BCS regime (a negative and small), the first 
correction to the non-interacting density profile |8j) can 
be calculated using perturbation theory. Interactions 
are treated at the Hartree level by adding the term 
— Anh 2 \a\n / 2m to the expression (|73|) for the local chemi- 
cal potential. The resulting density profile is compressed 
due to the effect of the attractive interaction and the 
Thomas-Fermi radii reduce according to the law 

V muif \ 3157^ / 

holding if kp\a\ <C 1, where kP F is the Fermi wavevector 
(9j of the non-interacting gas. 

Another interesting case is the BEC limit where one 
treats the interaction between dimers using the mean- 
field term [id = gdn/2 in the equation of state. The cou- 
pling constant = 2Trh 2 a,dd/ni is fixed by the molecule- 
molecule scattering length add — 0.60a. In this limit the 
density is given by the inverted parabola profile (Dalfovo 
et ai, 1999) and the Thomas- Fermi radii reduce to 

e> {15 a d d\ 1/5 Uho 

Ri = a ho -z-N . (76 

V 2 a ho J UJi 

In Fig. [TU] we show the results for the density profiles 
measured in situ in a harmonically trapped Fermi gas 
at low temperatures. The plotted profile is the double 
integrated density n^'(z) = f dxdy n(r), corresponding 
to the quantity measured in the experiment (Bartenstein 
et ai, 2004a). Very good agreement between experiment 
and theory is found at unitarity where one finds 

with R z given by Eq.fJTJ). The best fit to the experi- 
mental curve yields the value (3 = — 0.73^0 09 (Barten- 
stein et ai, 2004c). The attractive nature of interactions 
at unitarity is explicitly revealed in Fig. [TU] through the 
comparison with the density profile of a non-interacting 
gas. For a systematic comparison between experimental 
and theoretical results for the density profiles, see Perali, 
Pieri and Strinati (2004). 

A more recent experimental determination of /3, also 
based on the in situ measurement of the radius of the 
cloud, gives the value (3 = —0.54(5) (Partridge et ai, 
2006a). These measurements refer to a gas of 6 Li atoms. 
An important result consists in the agreement found with 
experiments on 40 K atoms, where j3 was determined by 
extrapolating to low temperature the measured values 
of the oscillator energy (Stewart et ai, 2006). The fact 
that the value of j3 does not depend on the atomic species 
considered is a further important proof of the universal 



24 



TABLE II: Experimental and theoretical values of the univer- 
sal parameter f3. The experimental results are taken from: (1) 
Tarruell et al. (2007); (2) Kinast et al. (2005); (3) Partridge 
et al. (2006a); (4) Bartenstein et al. (2004c), and (5) Stewart 
et al. (2006). The theoretical results are from: (6) Carlson et 
al. (2003); (7) Astrakharchik et al. (2004a); (8) Carlson and 
Reddy (2005) and (9) Perali, Pieri and Strinati (2004). 







P 


Exp. ( 6 Li) 


ENS (1) 


-0.59(15) 




Duke (2) 


-0.49(4) 




Rice (3) 


-0.54(5) 




Innsbruck (4) 


-0.73lg;S 


Exp. ( 40 K) 


JILA (5) 


-0.54^ 


Theory 


BCS mean field 


-0.41 




QMC (6,7,8) 


-0.58(1) 




T-matrix (9) 


-0.545 



behavior exhibited by these systems at unitarity (see Ta- 
ble [IT] for a list of available experimental and theoretical 
results of 0). 



B. Release energy and virial theorem 

In addition to the in situ density profiles a valuable 
source of information comes from the measurement of 
the energy after switching off the confining trap (release 
energy). This energy consists of the sum of the kinetic 
and the interaction term 

E re l — Ekin + E int , (78) 

and within LDA it is simply given by Eq. (|70p . 

The release energy was first measured along the 
crossover by Bourdel et al. (2004) on a gas of 6 Li atoms. 
The most recent experimental determination at unitar- 
ity yields the estimate f3 — —0.59(15) (Tarruell et al., 
2007). This value agrees with the one extracted from in 
situ measurements of the radii (see Table HH . 

A general relationship between the release energy ([70]) 
and the potential energy (|7l"j) can be derived with the 
help of the virial theorem. This theorem holds when 
the energy density e(n) has a polytropic (power law) de- 
pendence on the density: e(n) cx n 7+1 . The polytropic 
dependence characterizes the BEC limit, where 7 = 1, as 
well as the unitary limit, where 7 = 2/3. The theorem 
is derived by applying the number conserving transfor- 
mation n(r) — > (1 + a) 3 n[(l + a)r] to the density of the 
gas at equilibrium. By imposing that the total energy 
variation vanish to first order in a, one easily gets the 
result 

3jE rel = 2E ho , (79) 
which reduces to E re i = E^o at unitarity. 



C. Momentum distribution 

The momentum distribution of ultracold Fermi gases is 
an important quantity carrying a wealth of information 
on the role played by interactions. A simple theoretical 
approach for trapped systems is based on the BCS mean- 
field treatment introduced in Sec. IV Al and on the local 
density approximation (Viverit et al., 2004). The result 
is given by the spatial integral of the particle distribution 
function (|60|) of the uniform gas 

»00 = / ^a^W • ( 8 °) 

where the r-dependence enters through the chemical po- 
tential fi(r) and the order parameter A(r) which are de- 
termined by the local value of the density n(r). The 
momentum distribution n(k) is calculated for given val- 
ues of the interaction strength k° F a. Two limiting cases 
can be derived analitically: one corresponds to the non- 
interacting gas, where rik(r) = 0[1— fc/fcj?(r)] depends on 
the local Fermi wavevector A;f(i") = [37r 2 n(r)] 2 / 3 and the 
integral (|5D|) yields the result (fTTJjl . The other corresponds 
to the deep BEC regime, where, using Eq. (|rjl"j) . one finds 
the molecular result n(k) = (a 3 iV/27r 2 )/(l + k 2 a 2 ) 2 . A 
general feature emerging from these results is the broad- 
ening of the Fermi surface which, for trapped systems, 
is caused both by the confinement and by interaction ef- 
fects. 

The momentum distribution has been measured along 
the crossover in a series of studies (Regal et al., 2005; 
Chen et al., 2006 and Tarruell et al, 2007). The ac- 
cessible quantity in experiments is the column inte- 
grated distribution n co i(k^) = <ifc z n(k), where k±_ = 
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FIG. 11: (color online). Column integrated momentum dis- 
tributions of a trapped Fermi gas. The symbols correspond 
to the experimental results from Regal et al. (2005) and the 
lines to the mean-field results based on Eq. (|80[) for the same 
values of the interaction parameter l/k%a. From top to bot- 
tom l/k° F a = -71 (blue), l/k%a = (red) and l/k%a = 0.59 
(black). 



25 



■\Jk% + ky. These experiments are based on the technique 

of time-of-flight expansion followed by absorption imag- 
ing. A crucial requirement is that the gas must expand 
freely without any interatomic force. To this purpose the 
scattering length is set to zero by a fast magnetic-field 
ramp immediately before the expansion. The measured 
column integrated distributions along the crossover from 
Regal et al. (2005) are shown in Fig. [TT] together with 
the mean-field calculations of n co i based on Eq. (|80| for 
the same values of the interaction strength l/k F a (the 
value of a corresponds here to the scattering length before 
the magnetic- field ramp) . There is an overall qualitative 
agreement between the theoretical predictions and the 
experimental results. However, the dynamics of the ramp 
produces a strong quenching of n co i at large momenta 
k ^s> kp which greatly affects the released kinetic energy 
E re i — 2tt J dk±^k^n co i(k±^(ih 2 k\/ Am). A theoretical 
estimate of E re i, based on the equilibrium distributions 
would predict a very large value, on the order of the en- 
ergy scale fi 2 /mi?g associated with the interatomic po- 
tential [see discussion in Sec. IV Al after Eq. ([61]) ]. On the 
contrary, the measured energies from Regal et al. (2005) 
do not depend on the details of the interatomic potential 
because the magnetic-field ramp is never fast enough to 
access features on order of the interaction range Ro- A 
good quantitative agreement for the measured momen- 
tum distribution and kinetic energy is provided by the 
approach developed by Chiofalo, Giorgini and Holland 
(2006), which explicitly accounts for the time dependence 
of the scattering length produced by the magnetic-field 
ramp. 



D. Trapped gas at finite temperatures 

The local density approximation (|72[) for the density 
profile can be successfully used also at finite temperature, 
where the chemical potential depends on both the density 
and the temperature and is defined according to the ther- 
modynamic relationship /i(n,T) — (de(n,T)/dn) s where 
s = nS(n,T)/N is the entropy density. Since the transi- 
tion temperature decreases when the density decreases, 
the LDA predicts a sharp spatial boundary between a 
superfluid core and a normal external region. The ef- 
fect is particularly evident in the BEC regime where it 
gives rise to a typical bimodal distribution character- 
ized by a narrow Bose-Einstein condensate of molecules 
sorrounded by a broader cloud of thermal molecules. 
The experimental observation of this bimodal distribu- 
tion (Greincr, Regal and Jin, 2003; Bartenstein et al., 
2004a; Zwierlein et al, 2003; Bourdel et al, 2004; Par- 
tridge et al, 2005) represents the most spectacular and 
direct evidence of the emergence of Bose-Einstein con- 
densation from an interacting Fermi gas (see Fig. [I]) . In 
this BEC regime the critical temperature for the super- 
fluid transition is given by the well known expression 
k B T BEC = Hu> ho [N d /((3)] 1/3 , where N d = N/2 is the 
number of dimers and £(3) ~ 1.202. In terms of the 
Fermi temperature © this result reads 

7f E ° c = 0.52T#° , (81) 
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FIG. 12: Verifying the virial theorem at unitarity in a Fermi 
gas of 6 Li: (x 2 ) / (x 2 (0)} versus E/E(0) showing linear scaling. 
Here (x 2 ) is the transverse mean square radius, proportional 
to the oscillator energy. E is the total energy evaluated as in 
Kinast et al. (2005). -E(O) and (x 2 (0)} denote ground-state 
values. From Thomas, Kinast and Turlapov (2005). 
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FIG. 13: Measured entropy in a trapped 6 Li gas at unitarity 
vs. total energy. The data reveal the occurrence of a char- 
acteristic change of behavior in the slope at [E — E(0)]/N = 
0.41_Ef, which is interpreted as the signature of the phase 
transition to the superfluid regime. The solid lines are power- 
law fits below and above the critical point, while the dashed 
lines show the extended fits. From Luo et al. (2007). 
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and should be compared with the corresponding result 
(|37p holding in uniform gases. By including interaction 
effects to lowest order using a mean-field approach, one 
predicts a positive shift of the transition temperature 
(|8ip . proportional to the dimer-dimer scattering length 
add (Giorgini, Pitaevskii and Stringari, 1996). Due to 
large interaction effects in the condensate and in the ther- 
mal cloud, the bimodal structure becomes less and less 
pronounced as one approaches the resonance and conse- 
quently it becomes more difficult to reveal the normal- 
superfluid phase transition by just looking at the density 
profile. 

Once one knows the density profiles of the interact- 
ing Fermi g function of temperature and interac- 
tion strength, one can calculate the various thermody- 
namic functions along the crossover using the LDA. For 
example, the release energy and the entropy are given, 
respectively, by Eq. ([70]) and by the integral S(T) — 
f drs[T,n(r)]. At unitarity the relevant thermodynamic 
functions can be conveniently expressed in terms of the 
universal functions fp(x) and fs(x) defined in Sec lIVF"! 
where the ratio x = T/Tp is replaced by the local quan- 
tity /c B T2m/ft 2 [37r 2 n(r)] 2 / 3 (Ho, 2004; Thomas, Kinast 
and Turlapov, 2005). Furthermore, a useful result at uni- 
tarity, which keeps its validity at finite temperature, is 
the virial relationship E re \ = Eho [see Eq. ([79"]>]. Indeed, 
since the energy has a minimum at constant TV and S, by 
imposing the stationarity condition at constant T/n 2 / 3 
one immediately finds the above identity. 

The temperature dependence of the thermodynamic 
functions in a trapped gas along the crossover has been 
calculated in a series of papers using self-consistent 
many-body approaches (Perali et ai, 2004; Kinast et at, 
2005; Chen, Stajic and Levin, 2005; Stajic, Chen and 
Levin, 2005; Hu, Liu and Drummond, 2006a and 2006b). 
Other studies are based on quantum Monte Carlo tech- 
niques (Burovski et ai, 2006b; Bulgac, Drut and Magier- 
ski, 2007). In particular, the transition temperature has 
been determined using LDA and the Monte Carlo re- 
sults for uniform systems already discussed in connec- 
tion with Eq. igU) of Sec. IIVFI The reported value is 
T c = 0.20(2)S F (Burovski et ai, 2006b). Furthermore, 
the values of the energy and entropy per particle at the 
transition point have also been determined (Bulgac, Drut 
and Magierski, 2007), yielding respectively the results 
[E{T C ) - E(Q)]/N = 0.32£; F and S{T c )/Nk B = 2.15. 

From the experimental point of view a major difficulty 
in the study of the thermodynamic functions is the de- 
termination of the temperature when one cools down the 
system into the deeply degenerate regime. In fact, when 
the system is very cold the measurement of the density 
profile in the tails, which in general provides a natural 
access to T, is not accurate, the number of thermally 
excited atoms being small. This is not a severe prob- 
lem in the deep BCS regime, where the non-interacting 
Thomas- Fermi profile fitted to the whole spatial distribu- 
tion of the expanded cloud provides a reliable thermom- 
etry. Instead, in the most interesting strongly coupled 



regime, a model-independent temperature calibration is 
very difficult to obtain. Useful estimates of the tempera- 
tures achievable through adiabatic transformations along 
the BCS-BEC crossover can be obtained using entropy 
arguments. For example, starting from an initial config- 
uration in the molecular BEC regime with temperature 
T and changing adiabatically the value of the scattering 
length from positive to small negative values on the other 
side of the Feshbach resonance, one eventually reaches 
a final temperature in the BCS regime given by (Carr, 
Shlyapnikov, and Castin, 2004) 

( T -f) final 45 C(3) (t BE c) in i t i a i ^ ^ 

This relationship has been obtained by requiring that 
the entropies of the initial and final regimes, calculated 
using the predictions of the degenerate ideal Bose and 
Fermi gases respectively, be the same. The adiabatic 
transformation results in a drastic reduction of T. 

Many thermodynamic functions have already been 
measured in trapped Fermi gases. Measurements of the 
specific heat and a first determination of the critical tem- 
perature T c were reported by Kinast et al. (2005) who ex- 
tracted the value T c — 0.27Tp at unitarity. The value of 
T c was determined by identifying a characteristic change 
of slope of the measured energy as a function of tem- 
perature and relied on a model-dependent temperature 
calibration. 

In order to overcome the difficulties of the direct mea- 
surement of the temperature, recent experiments have 
also focused on the study of relevant thermodynamic re- 
lationships. One of these experiments concerns the verifi- 
cation of the virial relation at unitarity (Thomas, Kinast 
and Turlapov, 2005) that has been achieved by measur- 
ing independently the mean-square radius, proportional 
to the oscillator energy, and the total energy. The results 
shown in Fig. [T^] demonstrate that the virial relation is 
verified with remarkable accuracy, confirming the univer- 
sality of the equation of state at unitarity. 

Another experiment concerns the measurement of the 
entropy of the trapped gas (Luo et ai, 2007). In this 
experiment one starts from a strongly interacting config- 
uration (for example at unitarity) and switches off adia- 
batically the interactions bringing the system into a final 
weakly interacting state, where the measurement of the 
radii (and hence of the oscillator energy) is expected to 
provide a reliable determination of the entropy. Since en- 
tropy is conserved during the transformation, the mea- 
surement of the radii of the final (weakly interacting) 
sample determines the entropy of the initial (strongly in- 
teracting) configuration. This procedure has been used 
to measure the entropy as a function of the energy at 
unitarity in a gas of 6 Li atoms. The results are shown 
in Fig. [13J These measurements give access to the tem- 
perature of the gas through the relationship dE/dS = T. 
They have been used to study the thermodynamic be- 
havior near the critical point which can be identified 
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as a change in the energy dependence S(E) of the en- 
tropy. This method provides the experimental estimate 
T c = 0.29Tp for the critical temperature in good agree- 
ment with the analysis of the specific heat (Kinast et al, 
2005). Furthermore, the extracted values of the critical 
entropy, S(T c )/NkB = 2.7(2), and energy per particle, 
[E(T C ) - E(0)]/N = 0A1(5)E F , reasonably agree with 
the theoretical estimates reported above. 



VII. DYNAMICS AND SUPERFLUIDITY 



be equal (iVj — N^). In the same picture the current is 
given by j = riv. The equations take the form 

iU + V-(nv) = 0, (84) 
at 

for the density (equation of continuity) and 

m J; v + v ( \ mw2 +M n ) + V ho ) = , ( 85 ) 



Superfluidity is one of the most striking properties ex- 
hibited by ultracold Fermi gases, analog to superconduc- 
tivity in charged Fermi systems. Among the most no- 
ticeable manifestations of superfluidity one should recall 
the absence of shear viscosity, the hydrodynamic nature 
of macroscopic dynamics even at zero temperature, the 
Josephson effect, the irrotational quenching of the mo- 
ment of inertia, the existence of quantized vortices and 
the occurrence of pairing effects. The possibility of ex- 
ploring these phenomena in ultracold gases is providing 
a unique opportunity to complement our present knowl- 
edge of superfluidity in neutral Fermi systems, previously 
limited to liquid 3 He, where pairing occurs in a p-w&ve 
state. In this Section we will first discuss the hydrody- 
namic behavior exhibited by superfluids and its implica- 
tions on the dynamics of trapped Fermi gases (expansion 
and collective oscillations). The last part is devoted in- 
stead to a discussion of pair-breaking effects and of Lan- 
dau's critical velocity. Other manifestations of superflu- 
idity will be discussed in Sec. IVIIII 



A. Hydrodynamic equations of superfluids at T = 

The macroscopic behavior of a neutral superfluid is 
governed by the Landau equations of irrotational hydro- 
dynamics. The condition of irrotationality is the conse- 
quence of the occurrence of off-diagonal long-range order, 
characterized by the order parameter A(r) [see Eqs. (|26p- 
(|28|) ] . In fact the velocity field is proportional to the gra- 
dient of the phase S of the order parameter according 
to 

v = i V5 - (83) 

At zero temperature the hydrodynamic equations of 
superfluids consist of coupled and closed equations for 
the density and the velocity field. Actually, due to the 
absence of the normal component, the superfluid density 
coincides with the total density and the superfluid cur- 
rent with the total current. The hydrodynamic picture 
assumes that the densities of the two spin components are 
equal and move in phase [ri|(r, t) = ni(r,t) = n(r,t)/2 
and V|(r, t) = v^(r,t) = v(r, t)]. This implies, in par- 
ticular, that the number of particles of each spin species 



for the velocity field (Euler's equation). Here fj,(n) is 
the atomic chemical potential, fixed by the equation of 
state of uniform matter. At equilibrium (v = 0) Eu- 
ler's equation reduces to the LDA equation (fT2")) for the 
ground-state density profile. The hydrodynamic equa- 
tions (|84H85p can be generalized to the case of superfluids 
of unequal masses and unequal trapping potentials (see 
Scc. lDCCl) . 

Despite the quantum nature underlying the superfluid 
behavior, the hydrodynamic equations of motion have 
a classical form and do not depend explicitly on the 
Planck constant. This peculiarity raises the question 
whether the hydrodynamic behavior of a cold Fermi gas 
can be used to probe the superfluid regime. In fact, as 
we will see, Fermi gases above the critical temperature 
with resonantly enhanced interactions enter a collisional 
regime where the dynamic behavior is governed by the 
same equations ([84 |) - (|85|) . In this respect it is impor- 
tant to stress that collisional hydrodynamics admits the 
possibility of rotational components in the velocity field 
which are strictly absent in the superfluid. A distinc- 
tion between classical and superfluid hydrodynamics is 
consequently possible only studying the rotational prop- 
erties of the gas (see Sec. IVIII|) . Let us also emphasize 
that the hydrodynamic equations of superfluids have the 
same form both for Bose and Fermi systems, the effects 
of statistics entering only the form of fJ.(n). 

The applicability of the hydrodynamic equations is in 
general limited to the study of macroscopic phenomena, 
characterized by long wavelength excitations. In partic- 
ular, the wavelengths should be larger than the so-called 
healing length £. The value of £ is easily estimated 
in the BEC limit where the Bogoliubov spectrum ([55]) 
approaches the phonon dispersion law for wavelengths 
larger than £ ~ (na^) -1 / 2 . In the opposite BCS regime 
the healing length is instead fixed by the pairing gap. In 
fact, the Bogoliubov- Anderson phonon mode is defined 
up to energies of the order hck ~ A gQp , corresponding 
to £ - Hv F /A gap . In both the BEC and BCS limits the 
healing length becomes larger and larger as kp \a\ — > 0. 
Near resonance instead the only characteristic length is 
fixed by the average interatomic distance and the hydro- 
dynamic theory can be applied for all wavelengths larger 
than kp 1 . In Sec. lVIIDI we will relate the healing length 
to the critical Landau velocity and discuss its behavior 
along the crossover. 
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B. Expansion of a superfluid Fermi gas 



scaling ansatz 



In most experiments with ultracold atomic gases im- 
ages are taken after expansion of the cloud. These mea- 
surements provide information on the state of the gas in 
the trap and it is consequently of crucial importance to 
have an accurate theory describing the expansion. As 
discussed in Sec. IIIB| in the absence of interactions the 
expansion of a Fermi gas is asymptotically isotropic also if 
before the expansion the gas is confined by an anisotropic 
potential. Deviations from isotropy are consequently an 
important signature of the role of interactions. In the ex- 
periment of O'Hara et al. (2002) the first clear evidence 
of anisotropic expansion of a cold Fermi gas at unitar- 
ity was reported (see Fig. I14[) . opening an important de- 
bate in the community aimed to understand the nature 
of these novel many-body configurations. Hydrodynamic 
theory has been extensively used in the past years to 
analyze the expansion of Bose-Einstein condensed gases. 
More recently (Menotti, Pedri and Stringari, 2002) this 
theory was proposed to describe the expansion of a Fermi 
superfluid. The hydrodynamic solutions are obtained 
starting from the equilibrium configuration, correspond- 
ing to a Thomas-Fermi profile, and then solving Eq. (|85p 
by setting V/ lo — for t > 0. In general the hydrody- 
namic equations should be solved numerically. However, 
for an important class of configurations the spatial de- 
pendence can be determined analytically. In fact, if the 
chemical potential has the power-law dependence /lan 1 
on the density, then the Thomas-Fermi equilibrium pro- 
files have the analytic form oc (/xq — Vho) 1 ^ 1 and the 
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FIG. 14: (color online). Aspect ratio as a function of time 
during the expansion of an ultracold Fermi gas at unitarity. 
Upper (red) symbols: experiment; upper (red) line: hydrody- 
namic theory. For comparison the figure also shows the results 
in the absence of interactions. Lower (blue) symbols: experi- 
ment; lower (blue) line: ballistic expansion. From O'Hara et 
al. (2002). 



n(x,y,z,t) = (b x b y b z ) ln o(y, y, y) , (86) 

x y z 

provides the exact solution. The scaling parameters bi 
obey the simple time-dependent equations: 



k(b x b y b z )t 



. 



(87) 



The above equation generalizes the scaling law previously 
introduced in the case of an interacting Bose gas (7 = 1) 
(Castin and Dum, 1996; Kagan, Surkov and Shlyapnikov, 
1996). From the solutions of Eq. ([57)1 one can calculate 
the aspect ratio as a function of time. For an axially 
symmetric trap {uj x = uj y = lu±) this is defined as the 
ratio between the radial and axial radii. In terms of the 
scaling parameters bi it can be written as 

R±(t) _ b x (t) uj z 



Z(t) b z {t) u x ' 

For an ideal gas the aspect ratio tends to unity, while 
the hydrodynamic equations yield an asymptotic value 
^ 1. Furthermore, hydrodynamics predicts a peculiar 
inversion of shape during the expansion caused by the 
hydrodynamic forces which are larger in the direction of 
larger density gradients. As a consequence an initial cigar 
shaped configuration is brought into a disk shaped profile 
at large times and vice-versa. One can easily estimate 
the typical time at which the inversion of shape takes 
place. For a highly elongated trap (cu± ^S> UJ Z ) the axial 
radius is practically unchanged for short times since the 
relevant expansion time along the z axis is fixed by 1 /u) z . 
Conversely, the radial size increses fast and, for ui±t ^> 1 
one expects R±(t) ~ R±(0)ui±t. One then finds that the 
aspect ratio is equal to unity when uj z t ~ 1. 
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FIG. 15: Aspect ratio after expansion from the trap along the 
crossover. The solid lines are the hydrodynamic predictions 
at unitarity and for a gas in the deep BEC regime. The 
Fermi momentum kp in the figure corresponds to the value 
©. From Altmeyer et al. (2007c). 
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In Fig. [14] we show the predictions for the aspect ra- 
tio given by Eqs. (|87H88[) at unitarity, where 7 = 2/3, 
together with the experimental results of O'Hara et al. 
(2002) and the predictions of the ideal Fermi gas. The 
configuration shown in the figure corresponds to an ini- 
tial aspect ratio equal to R±/Z = 0.035. The compar- 
ison strongly supports the hydrodynamic nature of the 
expansion of these ultracold Fermi gases. The experi- 
ment was repeated at higher temperatures and found to 
exhibit a similar hydrodynamic behavior even at temper- 
atures of the order of the Fermi temperature, where the 
system can not be superfluid. One then concludes that 
even in the normal phase the hydrodynamic regime can 
be ensured by collisional dynamics. This is especially 
plausible close to unitarity where the scattering length 
is very large and the free path of atoms is of the order 
of the interatomic distances. The anisotropic expansion 
exhibited by ultracold gases shares important analogies 
with the expansion observed in the quark gluon plasma 
produced in heavy ion collisions (see, for example, Kolb 
et al., 2001; Shuryak, 2004). 

In more recent experiments (Tarruell et al., 2007; Alt- 
meyer et al., 2007c) the aspect ratio was measured at 
the coldest temperatures along the BCS-BEC crossover, 
at a fixed expansion time. At unitarity the results (see 
Fig. I15[) well agree with the hydrodynamic predictions 
but, as one moves towards the BCS regime, the aspect 
ratio becomes closer and closer to unity thereby revealing 
important deviations from the hydrodynamic behavior. 
In the deep BCS regime the hydrodynamic theory pre- 
dicts indeed the same behavior for the time dependence 
of the aspect ratio as at unitarity, the equation of state 
being characterized by the same power law coefficient 
7 = 2/3. The deviations from the hydrodynamic pic- 
ture follow from the fact that, since the density becomes 
smaller and smaller during the expansion, the BCS gap 
becomes exponentially small and the hydrodynamic pic- 
ture can not longer be safely applied. Superfluidity is ex- 
pected to be instead robust at unitarity since the energy 
gap is much larger than the typical oscillator frequency 
(whose inverse fixes the time scale of the expansion) and 
that, as a consequence, pairs can not easily break during 
the expansion. 

A case which deserves special attention is the expan- 
sion of the unitary gas for isotropic harmonic trapping. 
In this case an exact solution of the many-body problem 
is available (Castin, 2004) without invoking the hydro- 
dynamic picture. One makes use of the following scaling 
ansatz for the many-body wavefunction 



the kinetic energy, the Schrodinger equation is exactly 
satisfied by Eq. ([59")) provided b — uj 2 /b 3 . One conse- 
quently finds that the expansion of the unitary gas is ex- 
actly given by the non-interacting law b(t) — y/oj 2 a t 2 + 1, 
which coincides with the prediction of the hydrodynamic 
equations (|57|) after setting u>i = w; lo . 



C. Collective oscillations 

The collective oscillations of a superfluid gas pro- 
vide a unique source of information, being at the same 
time of relatively easy experimental access and of rele- 
vant theoretical importance. At zero temperature they 
can be studied by considering the linearized form of 
the time dependent hydrodynamic equations (|84|) - (|85|) . 
corresponding to small oscillations of the density, n — 
rio + Sn exp(— iwt), with respect to the equilibrium pro- 
file tiq{v). The linearized equations take the form 



- u> 5n 



1 



-V • 



?i V(— bn) 
an 



(90) 



In the absence of harmonic trapping the solutions cor- 
respond to sound waves with dispersion lu — cq and, in a 
Fermi superfluid, coincide with the Bogoliubov- Anderson 
phonons. In the BEC limit one recovers the Bogoliubov 
result c = \J nh 2 addnd/ni 2 for the sound velocity, where 
add is the dimer-dimer scattering length and rid = n/2 
is the molecular density. In the BCS limit one instead 
approaches the ideal gas value c = vf/V3. Finally, at 
unitarity one has the result (|4"3"]) . 

In the presence of harmonic trapping the propagation 
of sound is affected by the inhomogeneity of the medium. 
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*(r 1 ,...r Ar ,t) = N{ty^ r > h ' 2hb t! {v 1 /b,...r N /b) , 

(89) 

where ^0 is the many-body wavefunction at t = 0, b(t) is 
a time-dependent scaling variable and N(i) is a normal- 
ization factor. Eq. (|89|) exactly satisfies the Bethe-Peierls 
boundary condition (|25p if one works at unitarity where 
1/a = 0. At the same time, for distances larger than the 
range of the force where the Hamiltonian coincides with 



FIG. 16: Sound velocity in the center of the trap measured 
along the BCS-BEC crossover (symbols). The theory curves 
are obtained using Eq. (I9ip and are based on different equa- 
tions of state: BCS mean field (dotted line), quantum Monte 
Carlo (solid line) and Thomas-Fermi molecular BEC with 
add = 0.60a (dashed line). The Fermi momentum fej? in the 
figure corresponds to the value ((9} . From Joseph et al. (2006). 
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In a cylindrical configuration (lj z — 0) the sound veloc- 
ity can be easily calculated using the ID result rnc\ D — 
nmdnin/dniD, where n\r> = J dxdyn is the ID den- 
sity and the ID chemical potential /iiu is determined 
by the Thomas-Fermi relation + V/ lo (r^) = for 
the radial dependence of the density profile. One finds 
(Capuzzi el al, 2006) 

1D J dxdy(dfj,/dn)~ 1 ^) ^ ^ 

On the BEC side, where (J, oc n, one recovers the result 
c\o = c/v2 in terms of the value (|3"Tj) of uniform sys- 
tems. This result was first derived by Zaremba (1998) 
in the context of Bose-Einstein condensed gases. At uni- 
tarity, where [i oc n 2 / 3 , one instead finds c\o — y^3/5c. 
The propagation of sound waves in very elongated traps 
has been recently measured by Joseph et al. (2006). The 
experimental results are in good agreement with the es- 
timate of cid based on the QMC equation of state and 
clearly differ from the prediction of the BCS mean-field 
theory (see Fig. [TC]) . 

In the presence of 3D harmonic trapping the lowest 
frequency solutions of the hydrodynamic equations are 
discrctized, their frequencies being of the order of the 
trapping frequencies. These modes correspond to wave- 
lengths of the order of the size of the cloud. Let us 
first consider the case of isotropic harmonic trapping 
(lo x = ujy = lo z = uJho)- A general class of diver- 
gency free (also called surface) solutions is available in 
this case. They are characterized by the velocity field 
v oc V(r^Y£ m ), satisfying V • v = and corresponding to 
density variations of the form (dfi/dn)Sn oc r e Yi m . Y? m 



is here the spherical harmonic function. By using the 
identity {d^/dn)S/n Q = — Who, holding at T — for the 
density profile at equilibrium, one finds that these solu- 
tions obey the dispersion law uj(£) — ^/ItUho, independent 
of the form of the equation of state, as expected for the 
surface modes driven by an external force. This result 
provides a useful characterization of the hydrodynamic 
regime. The result in fact differs from the prediction 
uj(£) = Icuho of the ideal gas model, revealing the impor- 
tance of interactions accounted for by the hydrodynamic 
description. Only in the dipole case (£ — 1), correspond- 
ing to the center of mass oscillation, the interactions do 
not affect the frequency of the mode. 

In addition to surface modes an important solution is 
the I — 0, to = breathing mode whose frequency can 
be found in analytic form if the equation of state is poly- 
tropic (/i oc n 1 ). In this case the velocity field has the 
radial form v oc r and one finds cu(m = 0) = ^/S'y + 2u!f lo - 
For 7 = 1 one recovers the well known BEC result 
\/5w/io (Stringari, 1996b) while, at unitarity, one finds 
2ujho- It is worth stressing that the result at unitar- 
ity keeps its validity beyond the hydrodynamic approx- 
imation. The proof follows from the same arguments 
used for the free expansion in Sec. IVII Bl In fact, the 
scaling ansatz solves the Schrddinger equation also 
in the presence of isotropic harmonic trapping (Castin, 
2004). The scaling function in this case obeys to the 
equation b — ujf lo /b 3 — w^b, which admits undamped so- 
lutions of the form b(t) = A cos(2ujh t + 4>) + B. Here, 
A = \/B 2 — 1 and B (> 1) is proportional to the energy 
of the system (at equilibrium B = 1). In a BEC gas 
a similar result holds for the radial oscillation in cylin- 
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FIG. 17: Frequency of the radial quadrupole mode for an elon- 
gated Fermi gas in units of the radial frequency (upper panel) . 
The dashed line is the hydrodynamic prediction \^2loj_. The 
lower panel shows the damping of the collective mode. The 
Fermi momentum k f in the figure corresponds to the value 
©. From Altmeyer et al. (2007b). 



FIG. 18: (color online). Frequency of the radial compression 
mode for an elongated Fermi gas in units of the radial fre- 
quency (blue symbols). The curves refer to the equation of 
state based on BCS mean-field theory [lower (black) line] and 
on Monte Carlo simulations [upper (red) line] . The Fermi mo- 
mentum kF in the figure corresponds to the value ([9|. From 
Altmeyer et al. (2007a). 
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drical geometry (Pitaevskii, 1996; Kagan, Surkov and 
Shlyapnikov, 1996) and the corresponding mode was ex- 
perimentally investigated by Chevy et al. (2002). 

In the case of axisymmetric trapping (lu x = uj y = lo±_ =/= 
u>z) the third component hm of the angular momentum 
is still a good quantum number and one also finds simple 
solutions of Eq. (|90[) . The dipole modes, corresponding 
to the center of mass oscillation, have frequencies luj_ 
(to = ±1) and io z (to = 0). Another surface solution is 
the radial quadrupole mode (m = ±2) characterized by 
the velocity field v oc V(x ± iy) 2 and by the frequency 

w(m = ±2) = V2oj±_ , (92) 

independent of the equation of state. This collective 
mode has been recently measured by Altmeyer et al. 
(2007b). The experimental results (see Fig. [17]) show 
that while the agreement with the theoretical prediction 
(f9"2"| is very good on the BEC side, including unitarity, as 
soon as one enters the BCS regime strong damping and 
deviations from the hydrodynamic law take place, sug- 
gesting that the system is leaving the superfluid phase. 
For small and negative a the system is actually expected 
to enter a collisionless regime where the frequency of the 
to = ±2 quadrupole mode is 2wj_ apart from small mean- 
field corrections (Vichi and Stringari, 1999). Another 
class of surface collective oscillations, the so called scis- 
sors mode, will be discussed in Sec. IVIII Al due to its 
relevance for the rotational properties of the system. 

Differently from the surface modes, the to = com- 
pressional modes depend instead on the equation of state. 
For a polytropic dependence of the chemical potential 
(/i cx n 7 ) the corresponding solutions can be derived an- 
alytically. They are characterized by a velocity field of 
the form v cx V[a(x 2 + y 2 ) + bz 2 }, resulting from the 
coupling between the £ = 2 and £ = modes caused 
by the deformation of the trap (Cozzini and Stringari, 
2003). In the limit of elongated traps (ui z <C u>±) the 
two solutions reduce to the radial [lu = + 1)w_l] 

and axial [lu — (37 + 2)/(7 + l)w z ] breathing modes. 
Experimentally both modes have been investigated (Ki- 
nast et al., 2004; Bartenstein et al., 2004b; Altmeyer et 
al, 2007a). In Fig. [TH] we show the most recent exper- 
imental results (Altmeyer et al., 2007a) for the radial 
breathing mode. At unitarity the agreement between 
theory (yl0/3cj_i_ = 1.83o;j_) and experiment is remark- 
ably good. It is also worth noticing that the damping of 
the oscillation is smallest near unitarity. 

When we move from unitarity the collective oscilla- 
tions exhibit other interesting features. Theory predicts 
that in the deep BEC regime (7 = 1) the frequencies of 
both the axial and radial modes are higher than at unitar- 
ity. Furthermore, the first corrections with respect to the 
mean-field prediction can be calculated analytically, by 
accounting for the first correction to the equation of state 
Md = 9d n d produced by quantum fluctuations. This is the 
so called Lee-Huang- Yang (LHY) correction [see Eq. (I3"8"l) ] 
first derived in the framework of Bogoliubov theory of 



interacting bosons. The resulting shifts in the collective 
frequencies are positive (Pitaevskii and Stringari, 1998; 
Braaten and Pearson, 1999). As a consequence, when 
one moves from the BEC regime towards unitarity, the 
dispersion law exhibits a typical non monotonic behav- 
ior, since it first increases, due to the LHY effect, and 
eventually decreases to reach the lower value y/10/3tu± 
at unitarity (Stringari, 2004). 

In general, the collective frequencies can be calculated 
numerically along the crossover by solving the hydro- 
dynamic equations once the equation of state is known 
(Combescot and Leyronas, 2002; Hu et al., 2004; Heisel- 
berg, 2004; Kim and Zubarev, 2004; Bulgac and Bertsch, 
2005; Manini and Salasnich, 2005; Astrakharchick et al., 
2005b) . Fig. [TH] shows the predictions ( Astrakharchik et 
al., 2005b) obtained using the equation of state of the 
Monte Carlo simulations discussed in Sec. IV Bl and of the 
BCS mean-field theory of Sec. E3 The Monte Carlo 
equation of state accounts for the LHY effect while the 
mean-field BCS theory misses it, providing a monotonic 
behavior for the compressional frequency as one moves 
from the BEC regime to unitarity. The accurate mea- 
surements of the radial compression mode reported in 
Fig.HUconfirm the Monte Carlo predictions, providing an 
important test of the equation of state and first evidence 
for the LHY effect. Notice that the LHY effect is visible 
only at the lowest temperatures (Altmeyer et al., 2007a) 
where quantum fluctuations dominate over thermal fluc- 
tuations (Giorgini, 2000). This explains why previous 
measurements of the breathing mode failed in revealing 
the enhancement of the collective frequency above the 
value 2lo j_. 

The behavior of the breathing modes on the BCS side 
of the resonance exhibits different features. Similarly to 
the case of the quadrupole mode (see Fig [T7|) one expects 
that the system soon looses superfluidity and eventually 
behaves like a dilute collisionless gas whose collective fre- 
quencies are given by 2wj_ and 2lu z for the radial and 
axial modes respectively. Experimentally this transition 
is clearly observed for the radial mode (Kinast et al., 
2004, Bartenstein et al., 2004b), where it occurs at about 
kp\a\ ~ 1. It is also associated with a strong increase of 
the damping of the collective oscillation. 

The temperature dependence of the collective oscilla- 
tions has also been the object of experimental investiga- 
tions. Kinast, Turlapov and Thomas (2005) have shown 
that the frequency of the radial compression mode, mea- 
sured at unitarity, remains practically constant when one 
increases the temperature, suggesting that the system is 
governed by the same hydrodynamic equations both at 
the lowest temperatures, when the gas is superfluid, and 
at higher temperature when it becomes normal. Con- 
versely, the damping exhibits a significant temperature 
dependence, becoming smaller and smaller as one lowers 
the temperature. This behavior strongly supports the 
superfluid nature of the system in the low temperature 
regime. In fact a normal gas is expected to be less and 
less hydrodynamic as one decreases T with the conse- 



.32 



quent increase of the damping of the oscillation. The 
accurate determination of the damping of the collective 
oscillations can provide useful information on the viscos- 
ity coefficients of the gas whose behavior at unitarity has 
been the object of recent theoretical investigations us- 
ing universality arguments [see Son (2007) and references 
therein] . 



D. Phonons vs pair- breaking excitations and 
Landau's critical velocity 

In Sec. I VII O we have described the discretized modes 
predicted by hydrodynamic theory in the presence of har- 
monic trapping. This theory describes correctly only 
the low-frequency oscillations of the system, correspond- 
ing to sound waves in a uniform body. When one con- 
siders higher excitations energies the dynamic response 
should also include dispersive corrections to the phononic 
branch and the breaking of pairs into two fermionic exci- 
tations. The general picture of the excitations produced 
by a density probe in the superfluid state can then be 
summarized as follows (for simplicity we consider a uni- 
form gas): at low frequency the system exhibits a gapless 
phononic branch whose slope is fixed by the sound veloc- 
ity; at high frequency one expects the emergence of a 
continuum of excitations starting from a given thresh- 
old, above which one can break pairs. The value of the 
threshold frequency depends on the momentum hq car- 
ried by the perturbation. A first estimate is provided 
by the BCS mean-field theory discussed in Sec. IV Al 
according to which the threshold is given by the line 
LUthr = Mmk(e+ + e_), where e± is the energy ([5Tj) of a 
quasi-particle excitation carrying momentum fi,(q/2±k). 
The pair then carries total momentum hq. It is easy to 
see that the minimum is obtained for k • q = which 
gives e + = e_. One can distinguish two cases: fj, > and 
(i < 0. In the first case (including the unitary as well as 
the BCS regimes) and for hq < 2 v / 2m/7, the minimum is 
at h 2 k 2 /2m — (j, — h 2 q 2 /8m, so that huj t hr — 2A. For 
fi > and hq > 2y/2mfi, as well as for fi < (including 
the BEC regime), the minimum is at k = 0, which leads 
to hjJthr = 2y / (h 2 q 2 /8m - fl) 2 + A 2 . 

The interplay between phonon and pair breaking 
excitations gives rise to different scenarios along the 
crossover. In the BCS regime the threshold occurs at 
low frequencies and the phonon branch very soon reaches 
the continuum of single-particle excitations. The behav- 
ior is quite different in the opposite BEC regime where 
the phonon branch extends up to high frequencies. At 
large momenta this branch actually looses its phononic 
character and approaches the dispersion H 2 q 2 /4m, typi- 
cal of a free molecule. In the deep BEC limit it coincides 
with the Bogoliubov spectrum of a dilute gas of bosonic 
molecules. At unitarity the system is expected to exhibit 
an intermediate behavior, the discretized branch surviv- 
ing up to momenta of the order of the Fermi momentum. 
A detailed calculation of the excitation spectrum, based 



on a time-dependent formulation of the BCS mean-field 
equations is reported in Fig. [15] (Combescot, Kagan and 
Stringari, 2006). 

The results for the excitation spectrum provide a useful 
insight on the superfluid behavior of the gas in terms of 
the Landau criterion according to which a system can not 
give rise to energy dissipation if its velocity, with respect 
to a container at rest, is smaller than Landau's critical 
velocity defined by the equation 

v c = min q (omega q /q) , (93) 

where huj q is the energy of an excitation carrying momen- 
tum hq. According to this criterion the ideal Fermi gas 
is not superfluid because of the absence of a threshold 
for the single-particle excitations, yielding v c = 0. The 
interacting Fermi gas is instead superfluid in all regimes. 
By inserting the results for the threshold frequency de- 
rived above into Eq. (f9"3"|) one can calculate the critical 
value v c due to pair breaking. The result is 

m(v s c p f = VA 2 + M 2 - M - (94) 

and coincides, as expected, with the critical velocity cal- 
culated by applying directly the Landau criterion (|93l) 
to the single-particle dispersion law e q of Eq. (|5T|) . In 
the deep BCS limit kp\a\ — > 0, corresponding to A <C fi, 
Eq. (|94[) approaches the exponentially small value v c = 
A/hkp- On the BEC side of the crossover the value ([94]) 
becomes instead larger and larger and the relevant exci- 
tations giving rise to Landau's instability arc no longer 
single-particle excitations but phonons and the critical 
velocity coincides with the sound velocity: v c — c. A 
simple estimate of the critical velocity along the whole 
crossover is then given by the expression 

v c = Min(c,v s c p ) . (95) 

Remarkably one sees that v c has a maximum near uni- 
tarity (see Fig. [20]) , further confirming the robustness 
of superfluidity in this regime. This effect has been re- 
cently demonstrated experimentally by moving a one di- 
mensional optical lattice in a trapped superfluid Fermi 
at tunable velocity (see Fig. [21]) (Miller et al, 2007). 

The critical velocity permits to provide a general def- 
inition of the healing length according to £ = h/mv c . 
Apart from an irrelevant numerical factor it coincides 
with the usual definition hj v / 2ro/id °f t ne healing length 
in the BEC regime and with the size of Cooper pairs in 
the opposite BCS limit as already discussed in Sec. lVII Al 
The healing length provides the typical length above 
which the dynamic description of the system is safely 
described by the hydrodynamic picture. It is smallest 
near unitarity (Pistolesi and Strinati, 1996), where it is 
on the order of the interparticle distance. 

E. Dynamic and static structure factor 

The dynamic structure factor provides an important 
characterization of quantum many-body systems (see, 



33 




FIG. 19: Spectrum of collective excitations of the superfiuid Fermi gas along the BCS-BEC crossover (thick lines). Energy is 
given in units of the Fermi energy. Left: BCS regime [kpa — —1). Center: unitarity. Right: BEC regime {kpa = +1). The 
thin lines denote the threshold of single-particle excitations. From Combescot, Kagan and Stringari (2006). 



for example, Pines and Nozieres, 1966; Pitaevskii and 
Stringari, 2003). At low energy transfer it gives infor- 
mation on the spectrum of collective oscillations, includ- 
ing the propagation of sound, while at higher energies 
it is sensitive to the behavior of single-particle excita- 
tions. In general, the dynamic structure factor is mea- 
sured through inelastic scattering experiments in which 
the probe particle is weakly coupled to the many-body 
system so that the scattering may be described within the 
Born approximation. In dilute gases it can be accessed 
with stimulated light scattering by using two photon 
Bragg spectroscopy (Stamper-Kurn et at, 1999). The 
dynamic structure factor is defined by the expression 

S(q,a;)=Q- 1 ^e-^|(0|5p g |n)| 2 5(^;-^ ron ) , 

m,n 

(96) 



where fiq and fko are, respectively, the momentum and 
energy transferred by the probe to the sample, r5/5 q = 
Pi ~ (Pi) 1S t ne fluctuation of the Fourier component 
p q = 52 . exp(— iq • Tj) of the density operator, u> mn = 
Emn/fr = {E m — E n )/h are the usual Bohr frequencies 
and Q is the partition function. The definition of the 
dynamic structure factor is immediately generalized to 
other excitation operators like, for example, the spin den- 
sity operator. 

The main features of the dynamic structure factor are 
best understood in uniform matter, where the excitations 
are classified in terms of their momentum. From the re- 
sults of the Sec. I VHP I one expects that, for sufficiently 
small momenta, the dynamic structure factor be char- 
acterized by a sharp phonon peak and by a continuum 
of single-particle excitations above the threshold energy 
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FIG. 20: Landau's critical velocity (in units of the Fermi ve- 
locity) calculated along the crossover using BCS mean field 
theory. The figure clearly shows that the critical velocity is 
largest near unitarity. The dashed line is the sound velocity. 
From Combescot, Kagan and Stringari (2006). 
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FIG. 21: Measured critical velocity along the BCS-BEC 
crossover. The solid line is a guide to the eye. From Miller et 
al. (2007). 
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hjjthr- Measurements of the dynamic structure factor 
in Fermi superfluids can then provide unique informa- 
tion on the gap parameter. Theoretical calculations of 
the dynamic structure factor in the small q regime were 
carried out using a proper dynamic generalization of the 
BCS mean-field approach (Minguzzi, Ferrari and Castin, 
2001; Biichler, Zoller and Zwerger, 2004). At higher mo- 
mentum transfer the behavior will depend crucially on 
the regime considered. In fact, for values of q of the 
order of the Fermi momentum the discretized branch 
no longer survives on the BCS side of the resonance. 
At even higher momenta the theoretical calculations of 
Combescot, Giorgini and Stringari (2006) have revealed 
that on the BEC side of the resonance the response is 
dominated by a discretized peak correponding to the ex- 
citation of free molecules with energy h 2 q 2 /4m. This 
molecular-like peak has been shown to survive even at 
unitarity. On the BCS side of the resonance, the molec- 
ular signatures are instead completely lost at high mo- 
menta and the response is very similar to the one of an 
ideal Fermi gas. 

From the knowledge of the dynamic structure factor 
one can evaluate the static structure factor, given by 

S(q) = 1 f°° cLuS(q, u) = l + n[ dv[g{r) - l]e^ r 

(97) 

showing that the static structure factor is directly related 
to the two-body correlation function g = (g^ + g^)/2 
discussed in Sec. IV Bl The measuremnt of S(q) would 
then provide valuable information on the correlation ef- 
fects exhibited by these systems. The behavior of the 
static structure factor along the BCS-BEC crossover has 
been calculated by Combescot, Giorgini and Stringari 
(2006). 

F. Radiofrequency transitions 

In Sccs. lVII DtiVIIEI we have shown that pair-breaking 
transitions characterize the excitation spectrum associ- 
ated to the density-density response function and are 
directly visible in the dynamic structure factor. Infor- 
mation on pair-breaking effects is also provided by tran- 
sitions which outcouple atoms to a third internal state 
(Torma and Zoller, 2000). Experiments using radiofre- 
quency (RF) excitations have already been performed in 
Fermi superfluids (Chin et ai, 2004) in different condi- 
tions of temperature and magnetic fields. The basic idea 
of these experiments is the same as for the measurement 
of the binding energy of free molecules (see Sec. IIII A|l . 

The structure of the RF transitions is determined by 
the Zeeman diagram of the hyperfine states in the pres- 
ence of an external magnetic field. Starting from a sam- 
ple where two hyperfine states (hereafter called 1 and 
2) are occupied, one considers single-particle transitions 
from the state 2 (with energy E2) to a third, initially 
unoccupied, state 3 with energy £3. The typical excita- 
tion operator characterizing these RF transitions has the 



form Vrf = A / dr[^|(r)^2( r ) + h.c] and does not carry 
any momentum. The experimental signature of the tran- 
sition is given by the appearence of atoms in the state 3 
or in the reduction of the number of atoms in the state 
2. In the absence of interatomic forces the transition is 
resonant at the frequency 1^3 = {E3 — E2)/h. If instead 
the two atoms interact and form a molecule the frequency 
required to induce the transition is higher since part of 
the energy carried by the radiation is needed to break 
the molecule. The threshold for the transition is given 
by the frequency v = V23 + \eb\/h, where |ef,| ~ h 2 /ma 2 is 
the binding energy of the molecule (we are assuming here 
that no molecule is formed in the 1-3 channel). The ac- 
tual dissociation line shape is determined by the overlap 
of the molecular state with the continuum and is affected 
by the final-state interaction between the atom occupy- 
ing the state 3 and the atom occupying the state 1. In 
particular, the relationship between the value of the fre- 
quency where the signal is maximum and the threshold 
frequency depends on the value of the scattering length 
031 characterizing the interaction between the state 3 and 
the state 1 (Chin and Julienne, 2005). 

For interacting many-body configurations along the 
BCS-BEC crossover a similar scenario takes place. Also 
in this case breaking pairs costs energy so that the fre- 
quency shifts of the RF transitions provide information 
on the behavior of the gap, although this information is 
less direct than in the case of free molecules. The ideal 
situation would take place if the interaction between the 
final state 3 and the states 1 and 2 were negligible. In 
this case the threshold frequency is provided by the sim- 
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FIG. 22: RF spectra for various magnetic fields along the 
BCS-BEC crossover and for different temperatures. The RF 
offset 8p — v — V23 is given relative to the atomic transition 
2^3. The molecular limit is realized for B=72Q G (first 
column). The resonance regime is studied for B—822 G and 
837 G (second and third column). The data at 875 G (fourth 
column) correspond to the BCS side of the crossover. Upper 
row, signals of unpaired atoms at high temperature (larger 
than Tf); middle row, signals for a mixture of unpaired and 
paired atoms at intermediate temperature (fraction of Tf); 
lower row, signals for paired atoms at low temperature (much 
smaller than 2>). From Chin et al. (2004). 
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pic formula 

hdu = h(y — V23) — £fc=o — M ! (98) 

where we have considered the most favourable case where 
the RF transition excites single-particle states with zero 
momentum. Here [L is the chemical potential while Ck is 
the quasi-particle excitation energy defined in Sec. IIVBI 
Eq. reproduces the result hdv = \ei,\ in the deep 
BEC limit of free molecules where /1 approaches the value 
— I ef, 1/2. In the opposite BCS regime, where fi — > Ep and 
the single-particle gap A gap is much smaller than the 
Fermi energy, one has instead h8v = A^ ap /2Ep. The 
proper inclusion of final-state interactions in the calcu- 
lation of the RF spectra is a difficult problem that has 
been the object of several papers (Kinnunen, Rodriguez 
and Torma, 2004; He, Chen and Levin, 2005; Ohashi 
and Griffin, 2005; Yu and Baym, 2006; Perali, Pieri and 
Strinati, 2007). 

Typical experimental results on 6 Li are shown in 
Fig. [221 where the observed lineshapes are presented for 
different values of the temperature along the BCS-BEC 
crossover. In 6 Li the relevant scattering lengths 043 (< 0) 
and ai2 are both large in modulus and final-state inter- 
actions can not be ignored. On the BEC side of the res- 
onance (lowest magnetic field in Fig. [22)) one recognizes 
the emergence of the typical molecular lineshape at low 
temperature with the clear threshold effect for the RF 
transition. 

In the most interesting region where kp\a\ ~ 1 many- 
body effects become important and change the scenario 
of the RF transitions. While at high temperature (up- 
per row in Fig. [22[) the measured spectra still reveal the 
typical feature of the free atom transition, at lower tem- 
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FIG. 23: Frequency shift 5v ma x measured in Li at low tem- 
perature as a function of the magnetic field for two different 
configurations corresponding to Tf= 1.2 /iK (filled symbols) 
and 3.6 (open symbols). The solid line shows the value 
predicted in the free molecular regime where it is es- 
sentially given by the molecular binding energy. From Chin 
et al. (2004). 



peratures the lineshapes are modified by interactions in 
a nontrivial way. This is shown in Fig. |2"31 where the shift 
Svmax — v max — V23, defined in terms of the frequency 
Vmax where the RF signal is maximum, is reported for 
two different values of Tp. In the deep BEC regime the 
value of SiSmax is independent of the density and directly 
related to the binding energy of the molecule. Vicev- 
ersa, at unitarity and on the BCS side it shows a clear 
density dependence. From these data one extracts the 
relationship hSv max ~ 0.2Ep at unitarity. The depen- 
dence on the density is even more dramatic in the BCS 
regime, due to the exponential decrease of pairing effects 
as fcp|a| — * 0. 

Since pairing effects are density dependent and become 
weaker and weaker as one approaches the border of the 
atomic cloud, one can not observe any gap in Fig. 1221 ex- 
cept on the BEC side of the resonance where the gap is 
density independent. Spatially resolved RF spectroscopy 
has been recently become available at unitarity (Shin et 
at, 2007) revealing in a clear way the occurrence of the 
gap and hence opening new perspectives for a direct com- 
parison with the theoretical predictions in uniform mat- 
ter. 



VIII. ROTATIONS AND SUPERFLUIDITY 

Superfluidity shows up in spectacular rotational prop- 
erties. In fact a superfluid can not rotate like a rigid body, 
due to the irrotationality constraint (|83[) imposed by the 
existence of the order parameter. At low angular velocity 
an important macroscopic consequence of superfluidity is 
the quenching of the moment of inertia. At higher an- 
gular velocities the superfluid can instead carry angular 
momentum via the formation of vortex lines. The cir- 
culation around these vortex lines is quantized. When 
many vortex lines are created a regular vortex lattice is 
formed and the angular momentum acquired by the sys- 
tem approaches the classical rigid-body value. Both the 
quenching of the moment of inertia and the formation 
of vortex lines have been the object of fundamental in- 
vestigations in the physics of quantum liquids and have 
been recently explored in a systematic way also in dilute 
Bose-Einstein condensed gases. In this Section we sum- 
marize some of the main rotational features exhibited by 
dilute Fermi gases where first experimental results are al- 
ready available. In Secs. lVIII Al and [VIIIBI we discuss the 
consequences of the irrotationality constraint on the mo- 
ment of inertia, on the collective oscillations and on the 
expansion of a rotating gas, while Sec. IVIII CI is devoted 
to discuss some key properties of quantized vortices and 
vortex lattices. 



A. Moment of inertia and scissors mode 

The moment of inertia relative to the z axis is de- 
fined as the response of the system to a rotational field 
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— flL z according to (L z ) — HQ where L z is the z compo- 
nent of the angular momentum operator and the average 
is taken on the stationary configuration in the presence 
of the perturbation. For a non-interacting gas trapped 
by a deformed harmonic potential the moment of inertia 
can be calculated explicitly. In the small f2 limit (linear 
response) one finds the result (Stringari 1996a) 



e = 



mN 



0J„(X 



2 »] . 
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(99) 



where the expectation values should be evaluated in the 
absence of rotation. This result applies both to bosonic 
and fermionic ideal gases. It assumes uj x ^ uo y , but ad- 
mits a well defined limit when uj x — > oj y . In the ideal 
Fermi gas, when the number of particles is large, one 
can use the Thomas- Fermi relationships (x 2 ) cx 1/oj 2 and 
{y 2 ) (x 1/uj 2 for the radii. In this case Eq. (|99j) reduces 
to the rigid value of the moment of inertia: 



e ng = Nm(x 2 +y 2 ) . 



For a non-interacting Bose-Einstein condensed gas at 
T = 0, where the radii scale according to (a; 2 ) cx 1/oj x and 
{y 2 } oc 1/uj x , one instead finds that Q — > as u) x — ► u> y . 

Interactions change the value of the moment of inertia 
of a Fermi gas in a profound way. To calculate in the 
superfluid phase one can use the equations of irrotational 
hydrodynamics developed in Sec. IVII1 by considering a 
trap rotating with angular velocity f2 and looking for the 
stationary solution in the rotating frame. The equations 
of motion in the rotating frame are obtained by includ- 
ing the term — flL z in the Hamiltonian. In this frame 
the trap is described by the time independent harmonic 
potential V% of Eq. §Q and the hydrodynamic equations 
admit stationary solutions characterized by the velocity 
field 



v = -(JfiV(xy) 



(101) 



where v is the irrotational superfluid velocity in the labo- 
ratory frame, while 5 — (y 2 — x 2 ) / (y 2 +x 2 ) is the deforma- 
tion of the atomic cloud in the xy plane. By evaluating 
the angular momentum (L z ) = m J dr(r x v)n one finds 
that the moment of inertia takes the irrotational form 



Q = 6 2 Q 



rig , 



(102) 



which identically vanishes for an axisymmetric configura- 
tion, pointing out the crucial role played by superfluidity. 
The result for the moment of inertia holds for both Bose 
and Fermi superfluids rotating in a harmonic trap. 

The behavior of the moment of inertia at high angular 
velocity was investigated in the case of BEC's by Recati, 
Zambelli and Stringari (2001) who found that, for angu- 
lar velocities larger than where lo 2 ± = (uj 2 +oj 2 )/2, 
the adiabatic increase of the rotation can sizably affect 
the value of the deformation parameter 6 yielding large 



deformations even if the deformation of the trap is small. 
Physically this effect is the consequence of the energetic 
instability of the quadrupole oscillation. At even higher 
angular velocites a dynamic instability of the rotating 
configuration was predicted by Sinha and Castin (2001) 
suggesting a natural route to the nucleation of vortices. 
These theoretical predictions were confirmed experimen- 
tally (Madison et at, 2001). Similar predictions have 
been made also for rotating Fermi gases (Tonini, Werner 
and Castin, 2006). 

The irrotational nature of the moment of inertia has 
important consequences on the behavior of the so called 
scissors mode. This is an oscillation of the system caused 
by the sudden rotation of a deformed trap which, in 
the superfluid case, has frequency (Guery-Odelin and 
Stringari, 1999) 



CO 



(103) 



This result should be compared with the prediction of the 
normal gas in the collisionless regime where two modes 
(100) with frequencies ui± = \oj x ± oj y \ are found. The oc- 

j„\ reflects 



currence of the low frequency mode \u x — 
the rigid value of the moment of inertia in the normal 
phase. The scissors mode, previously observed in a Bose- 
Einstein condensed gas (Marago et ai, 2000), has been 
recently investigated also in ultracold Fermi gases (see 
Fig. [MD along the BCS-BEC crossover (Wright et at, 
2007) . At unitarity and on the BEC side of the resonance 
one clearly observes the hydrodynamic oscillation, while 
when a becomes negative and small the beating between 
the frequencies lu± = \lu x ± uj y \ reveals the transition 
to the normal collisionless regime. If the gas is normal, 
but deeply collisional as happens at unitarity above the 
critical temperature, classical hydrodynamics predicts an 
oscillation with the same frequency (| 1 03[) in addition to a 
low frequency mode of diffusive nature caused by the vis- 
cosity of the fluid. This mode, however, is located at too 
low frequencies to be observable. The persistence of the 
scissors frequency f| 103[) has been observed at unitarity in 
the recent experiment of Wright et at. (2007) even above 
T c . This result, together with the findings for the aspect 
ratio of the expanding gas and for the radial compression 
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FIG. 24: Time evolution of the angle characterizing the scis- 
sors mode in a Fermi gas at unitarity (left panel) and on the 
BCS side of the resonance (right panel). The measured fre- 
quencies well agree with the theoretical predictions (see text). 
From Wright et al. (2007). 
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mode (see discussion in Sees. IVIIBllVII Cj) . confirms that 
near resonance the gas behaves hydrodynamically in a 
wide range of temperatures below and above the critical 
temperature. This makes the distinction between the su- 
pcrfluid and the normal phase based on the study of the 
collective oscillations, a difficult task. 

Promising perspectives to distinguish between super- 
fluid and classical (collisional) hydrodynamics are pro- 
vided by the study of the collective oscillations excited 
in a rotating gas. In fact, in the presence of vorticity 
V x v 7^ 0, the equations of collisional hydrodynam- 
ics contain an additional term depending on the curl of 
the velocity field, which is absent in the equations of su- 
pcrfluid hydrodynamics. The resulting consequences on 
the behavior of the scissors mode have been discussed by 
Cozzini and Stringari (2003). 



B. Expansion of a rotating Fermi gas 

Another interesting consequence of the irrotational na- 
ture of the superfluid motion concerns the expansion of 
a rotating gas. Let us suppose that a trapped super- 
fluid Fermi gas is initially put in rotation with a given 
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FIG. 25: (color online). Aspect ratio versus expansion time 
of a unitary gas for different values of the initial angular ve- 
locity fio (in units of the frequency lo z of the long axis of the 
trap) and of the temperature of the gas (parametrized by the 
energy per particle E/N). Squares (purple): no angular ve- 
locity; solid circles (blue): Q /u) z = 0.40, E/(NE F ) = 0.56; 
open circles (red): Q,q/ijJ z = 0.40, E/(NEf) = 2.1; triangles 
(green): Q. /tu z = 1.12, E/(NE F ) = 0.56. The dashed, solid 
and dotted lines are the results of calculations based on ir- 
rotational hydrodynamics. From Clancy, Luo and Thomas 
(2007). 



value of angular velocity (in practice this can be realized 
through a sudden rotation of a deformed trap which ex- 
cites the scissors mode). The gas is later released from 
the trap. In the plane of rotation the expansion along 
the short axis of the cloud will initially be faster than 
the one along the long axis, due to larger gradients in 
the density distribution. However, differently from the 
hydrodynamic expansion of a non rotating gas where the 
cloud takes at some time a symmetric shape, in the ro- 
tating case the deformation of the cloud can not vanish. 
In fact, due to the irrotational constraint, this would re- 
sult in a vanishing value of the angular momentum and 
hence in a violation of angular momentum conservation. 
The consequence is that the angular velocity of the ex- 
panding cloud will increase as the value of the deforma- 
tion is reduced, but can not become too large because 
of energy conservation. The result is that the deforma- 
tion parameter will acquire a minimum value during the 
expansion but will never vanish. This non trivial conse- 
quence of irrotationality was first predicted (Edwards et 
al., 2002) and observed (Hcchcnblaikner et al., 2002) in 
Bose-Einstein condensed atomic gases. Very recently the 
experiment has been repeated in a cold Fermi gas at uni- 
tarity (Clancy, Luo and Thomas, 2007). Fig. [25] reports 
the measured aspect ratio as a function of the expansion 
time for different initial values of the angular velocity. It 
shows that, if the gas is initially rotating (lower curves), 
the aspect ratio never reaches the value 1 corresponding 
to a vanishing deformation. The experimental data are 
very well reproduced by the solutions of the equations 
of irrotational hydrodynamics (solid and dotted line) . A 
remarkable feature is that the same behavior for the as- 
pect ratio is found not only in the superfluid regime but 
also above the critical temperature, revealing that even 
in the normal phase the dynamics of the expansion is de- 
scribed by the equations of irrotational hydrodynamics. 
The reason is that viscosity effects are very small in the 
normal phase and that the relevant time scales in this 
experiment are too short to generate a rigid component 
in the velocity field. 



C. Quantized vortices 

The existence of quantized vortices is one of the most 
important predictions of superfluidity. Recent experi- 
ments have confirmed their existence also in ultracold 
Fermi gases along the BCS-BEC crossover (see Fig. [26l). 
In these experiments vortices are produced by spinning 
the atomic cloud with a laser beam and are observed by 
imaging the released cloud of molecules, which are sta- 
bilized through a rapid sweep of the scattering length to 
small and positive values during the first ms of the ex- 
pansion. This technique, which is similar to the one em- 
ployed to measure the condensate fraction of pairs (see 
Sec. IV Bj) . increases the contrast of the vortex cores and 
therefore their visibility. 

A quantized vortex along the z axis is associated with 
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the appearence of a phase in the order parameter [see 
Eqs. (|27 p -(|28 p ] of the form exp(z<^) where cf> is the az- 
imuthal angle. This yields the complex form 



A(r) = A(r±, z) exp(i<j 



(104) 



for the order parameter A where, for simplicity, we have 
assumed that the system exhibits axial symmetry and we 
have used cylindrical coordinates. The velocity field v = 
(h/2m)'V(p of the vortex configuration has a tangential 
form (v • rj_ = 0) with modulus v = h/2mr±, which 
increases as one approaches the vortex line, in contrast to 
the rigid body dependence v = f2 x r characterizing the 
rotation of a normal fluid. The circulation is quantized 
according to the rule 



vd£ = — 
m 



(105) 



which is smaller by a factor 2 with respect to the case of 
a Bose superfluid with the same atomic mass m. Vortices 
with higher quanta of circulation can also be considered. 
The value of the circulation is independent of the radius 
of the contour. This is a consequence of the fact that 
the vorticity is concentrated on a single line and hence 
deeply differs from the vorticity Vxv = 2CI of the rigid 
body rotation. 

The angular momentum carried by the vortex is given 
by the expression 



(L z ) = m I dr(r x v)n(r) = N 



(106) 



holding if the vortex line coincides with the symmetry 
axis of the density profile. If the vortex is displaced to- 
wards the periphery of a trapped gas the angular mo- 
mentum takes a smaller value. In this case the axial 
symmetry of the problem is lost and the order parameter 
can not be written in the form (|104| . 

A first estimate of the energy of the vortex line is ob- 
tained using macroscopic arguments based on hydrody- 
namics and considering, for simplicity, a gas confined in 
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FIG. 26: Experimental observation of quantized vortices in a 
superfluid Fermi gas along the BCS-BEC crossover (Zwierlein 
et al, 2005b). 



a cylinder of radial size R. The energy E v acquired by 
the vortex is mainly determined by the hydrodynamic 
kinetic energy (m/2)n J drv 2 which yields the following 
estimate for the vortex energy: 



E v = 



Nh R 
AmR 2 U J 



(107) 



where we have introduced the radius £ of the core of the 
vortex of the order of the healing length (see Sec. lVIIDj) . 
The need for the inclusion of the core size £ in Eq. Q107P 
follows from the logarithmic divergence of the integral 
n J drv 2 at short radial distances. Eq. (|107|) can be 
used to evaluate the critical angular velocity J7 C for the 
existence of an energetically stable vortex line. This 
value is obtained by imposing that the change in the 
energy E — tt c (L z ) of the system in the frame rotat- 
ing with angular velocity f2 c be equal to E v . One finds 
il c = (h/2mR 2 )\n(R/£). By applying this estimate to 
a harmonically trapped configuration with Rtf ~ R 
and by neglecting the logarithmic term which provides 
a correction of order of unity, we find Q c /uj± ~ hw±/Eh 
where is the radial frequency of the harmonic po- 
tential and Eho ~ rnio\R\ F is the harmonic oscillator 
energy. The above estimate shows that in the Thomas- 
Fermi regime, Eho 3> fiuj±, the critical frequency is much 
smaller than the radial trapping frequency, thereby sug- 
gesting that vortices should be easily produced in slowly 
rotating traps. This conclusion, however, does not take 
into account the fact that the nucleation of vortices is 
strongly inhibited at low angular velocities by the occur- 
rence of a barrier. For example in rotating Bose-Einstein 
condensates it has been experimentally shown that it is 
possible to increase the angular velocity of the trap up to 
values higher than f2 c without generating vortical states. 
Under these conditions the response of the superfluid is 
governed by the equations of irrotational hydrodynamics 
(see Sec. lVIIIAjl . 

A challenging problem concerns the visibility of the 
vortex lines. Due to the smallness of the healing length, 
especially at unitarity, they can not be observed in situ, 
but only after expansion. Another difficulty is the re- 
duced contrast in the density with respect to the case of 
Bose-Einstein condensed gases. Actually, while the order 
parameter vanishes on the vortex line the density does 
not, unless one works in the deep BEC regime. In the op- 
posite BCS regime the order parameter is exponentially 
small and the density profile is practically unaffected by 
the presence of the vortex. 

The explicit behavior of the density near the vortex 
line as well as the precise calculation of the vortex en- 
ergy requires the implementation of a microscopic calcu- 
lation. This can be carried out following the lines of the 
mean- field BCS theory developed in Sec. lVA"! While this 
approach is approximate, it nevertheless provides a useful 
consistent description of the vortical structure along the 
whole crossover. The vortex is described by the solution 
of the Bogoliubov- de Gennes Eq. (|4T))) corresponding to 
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the ansatz 



Ui(r) = u n {r ± )e~ lm *e lk > z /^/2^L 

Vi (r) = v n (r ± )e- t( - m+1)4, e^ z /V2^L , (108) 

for the normalized functions u% and vi where (n, m, k z ) 
are the usual quantum numbers of cylindrical symmetry 
and L is the length of the box in the z direction. The 
ansatz (|108|) is consistent with the dependence (|104[) of 
the order parameter A on the phase <f>. Calculations of 
the vortex structure based on the above approach, have 
been carried out by several authors (Nygaard et ai, 2003; 
Machida and Koyama, 2005; Chien et ai, 2006; Sen- 
sarma, Randeria and Ho, 2006). A generalized version of 
the Bogoliubov-de Gennes equations based on the den- 
sity functional theory was used instead by Bulgac and Yu 
(2003). An important feature emerging from these calcu- 
lations is that, near the vortex line, the density contrast 
is reduced at unitarity with respect to the BEC limit and 
is absent in the BCS regime. 

At higher angular velocities more vortices can be 
formed giving rise to a regular vortex lattice. In this 
limit the angular momentum acquired by the system ap- 
proaches the classical rigid-body value and the rotation is 
similar to the one of a rigid body, characterized by the law 
V x v = 20. Using result (|105j) and averaging the vortic- 
ity over several vortex lines one finds Vxv = (Trh/m)n v z, 
where n v is the number of vortices per unit area, so that 
the density of vortices n v is related to the angular veloc- 
ity by the relation n v = 2mQ/Tth showing that the dis- 
tance between vortices (proportional to 1/ y/ri^) depends 
on the angular velocity but not on the density of the gas. 
The vortices form thus a regular lattice even if the av- 
erage density is not uniform as happens in the presence 
of harmonic trapping. This feature, already pointed out 
in the case of Bose-Einstein condensed gases, has been 
confirmed in the recent experiments on Fermi gases (see 
Fig. [26j) . It is worth noticing that, due to the repulsive 
quantum pressure effect characterizing Fermi gases, one 
can realize trapped configurations with a large size R± 
hosting a large number of vortices N v — T:R 2 L n V) even 
with relatively small values of fi. For example choosing 
VL = oj±/3, ui z ~ uj± and N — 10 6 , one predicts N v ~ 130 
at unitarity which is significantly larger than the number 
of vortices that one can produce in a dilute Bose gas with 
the same angular velocity. 

At large angular velocity the vortex lattice is responsi- 
ble for a bulge effect associated with the increase of the 
radial size of the cloud and consequently with a modifi- 
cation of the aspect ratio. In fact, in the presence of an 
average rigid rotation, the effective potential felt by the 
atoms is given by Vho — (m/2)£l 2 r 2 L . The new Thomas- 
Fermi radii satisfy the relationship 



Rl 



Q 2 



(109) 



showing that at equilibrium the angular velocity can not 
overcome the radial trapping frequency. This formula 



can be used to determine directly the value of f2 by just 
measuring the in situ aspect ratio of the atomic cloud. 

Important consequences of the presence of vortex 
lines concern the frequency of the collective oscillations. 
For example, using a sum-rule approach (Zambclli and 
Stringari, 1998) it is possible to show that the splitting 
Aw between the to = ±2 quadrupole frequencies is given 
by the formula 



Alu = uj(m = +2) - uj(m = -2) = 2 



e z 



m(rj_) 



(110) 



where l z = (L z )/N is the angular momentum per par- 
ticle carried by the vortical configuration. For a single 
vortex line l z is equal to H/2, while for a vortex lattice 



l z is given by the rigid-body value f2m( 



The split- 



ting (|110|) . and hence the angular momentum i z can be 
directly measured by producing a sudden quadrupole de- 
formation in the x — y plane and imaging the correspond- 
ing precession <f> = Alu /4 of the angle (f> of the symmetry 
axis of the deformation during the quadrupole oscillation. 
This precession phenomenon has been already observed 
in the case of Bose-Einstein condensates containing quan- 
tized vortex lines (Chevy, Madison and Dalibard, 2000). 
For a single vortex line this experiment gives direct ac- 
cess to the quantization of the angular momentum (jl06|) 
carried by the vortex, in analogy with the famous Vinen 
experiment of superfluid helium (Vinen, 1961). 

In the presence of many vortex lines the collective os- 
cillations of the system can be calculated using the equa- 
tions of rotational hydrodynamics (Cozzini and Stringari, 
2003). In fact, from a macroscopic point of view, the vor- 
tex lattice behaves like a classical body rotating in a rigid 
way. For example, in the case of axisymmetric configu- 
rations one finds the result 



a; (to = ±2) 



2u! 2 ± -n 2 ±n 



(111) 



for the frequencies of the two to = ±2 quadrupole modes, 
which is consistent with the sum-rule result (|110|) for the 
splitting in the case of a rigid rotation. Other impor- 
tant modes affected by the presence of the vortex lattice 
are the compressional to = oscillations resulting from 
the coupling of the radial and axial degrees of freedom. 
These modes were discussed in Sec. I VII Cl in the absence 
of rotation. The effect of the rotation in a Fermi gas has 
been recently discussed by Antezza, Cozzini and Stringari 
(2007). In the centrifugal limit — ► lu± the frequency of 
the radial mode approaches the value u> = 2lu± indepen- 
dent of the equation of state, while the frequency of the 
axial mode approaches the value 



= y / J + 2u) z 



(112) 



where 7 is the coefficient of the polytropic equation of 
state (see Sec. lVII C|) . 

It is finally worth noticing that the achievement of the 
centrifugal limit for a superfluid containing a vortex lat- 
tice can not be ensured on the BCS side of the resonance. 
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In fact, due to the bulge effect (]109[) . the centrifugal limit 
is associated with a strong decrease of the density and 
hence, for a < 0, with an exponential decrease of the or- 
der parameter A. It follows that the superfluid can not 
support rotations with values of Q too close to lu±_ and 
that the system will exhibit a transition to the normal 
phase (Zhai and Ho, 2006; Veillette et ai, 2006) [see also 
Moller and Cooper (2007) for a recent discussion of the 
new features exhibited by the rotating BCS Fermi gas]. 
If f2 becomes too close to u>±, superfluidity will be even- 
tually lost also at resonance and on the BEC side of the 
resonance because the system enters the Quantum Hall 
regime (for the case of bosons see, e.g., Cooper, Wilkin 
and Gunn, 2001; Regnault and Jolicoeur, 2003) 



IX. SPIN POLARIZED FERMI GASES AND 
FERMI MIXTURES 



The description of Fermi superfluidity presented in the 
previous Sections was based on the assumption that the 
gas has an equal number of atoms occupying two dif- 
ferent spin states. One can also consider more complex 
configurations of spin imbalance where the number of 
atoms in the two spin states is different (JVj ^ iVj,) as 
well as mixtures of atomic species with different masses 
(m-\ =/= m-i), including Fermi- Fermi and Bose- Fermi mix- 
tures. Recent realizations of these novel configurations 
are opening new perspectives for both experimental and 
theoretical research. 




FIG. 27: Qualitative phase diagram as a function of the in- 
teraction strength — l/fc^a and of the polarization P. The 
circle at unitarity corresponds to the critical value P c = 0.39 
discussed in Sec. IIX Bl The Fermi wavevector corresponds 



here to the total average density: k F = [37r (n-j + n{)] 



1/3 _ 

Notice that the possible occurrence of the FFLO phase is not 
considered on this diagram. 



A. Equation of state of a spin polarized Fermi gas 

The problem of spin imbalance has an old story in the 
context of BCS theory of superconductivity. In super- 
conductors, due to the fast relaxation between different 
spin states leading to balanced spin populations, the only 
possibility to create the asymmetry is to add an exter- 
nal magnetic field. However, in bulk superconductors 
this field is screened by the orbital motion of electrons 
(Meissner effect). The situation in superfluid Fermi gases 
is more favourable. In fact, in this case, the relaxation 
time is very long and the numbers of atoms occupying 
different spin states can be considered as independent 
variables. 

Let us recall that the mechanism of BCS superfluidity, 
in the regime of small and negative scattering lengths, 
arises from the pairing of particles of different spin occu- 
pying states with opposite momenta, close to the Fermi 
surface. This mechanism is inhibited by the presence 
of spin imbalance, since the Fermi surfaces of the two 
components do not coincide and pairs with zero total 
momentum are difficult to form. Eventually, if the gap 
between the two Fermi surfaces is too large, superfluid- 
ity is broken and the system undergoes a quantum phase 
transition towards a normal state. The existence of a 
critical value for the polarization is easily understood by 
noticing that, at zero temperature, the unpolarized gas 
is superfluid, while a fully polarized gas is normal due 
to the absence of interactions. The occurrence of such a 
transition was first suggested in the papers by Clogston 
(1962) and Chandrasekhar (1962) who predicted the oc- 
currence of a first order transition from the normal to 
the superfluid state. This transition takes place when 
the gain in the grand-canonical energy associated with 
the finite polarization of the normal phase is equal to 
the energy difference between the normal and the super- 
fluid unpolarized states. In the BCS regime one finds the 
critical condition [see Eq. (|36[)] 



h 



Ml - Ml 



V2 



(113) 



where A gap = (2/e) 1 ^E F e" T ' 2kFa is the BCS gap and 
we have assumed the spin-down particles as the minority 
component. The chemical potential difference h in the 
above equation plays the role of an effective magnetic 
field. In terms of the polarization 



P = 



(114) 



it can be expressed through the relationship h = 
2E F P/3, which holds if P < 1. The condition (fm]) 
then immediately yields the critical value of P at which 
the system phase separates: 



7/3 



7T Ilk-pa 



(115) 
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For P > P c the system is normal and corresponds to 
a uniform mixture of the two spin components well de- 
scribed by the non-interacting model. For P < P c the 
system is instead in a mixed state, where the unpolarizcd 
BCS superfluid coexists with the normal phase which ac- 
commodates the excess polarization. In this mixed state, 
the chemical potential difference of the normal phase re- 
tains the critical value (|113p irrespective of polarization, 
a decrease in P being accounted for by an increase in 
the volume fraction of the superfluid phase which even- 
tually occupies the entire volume for P = 0. An im- 
portant remark concerning the Clogston-Chandrasekhar 
condition (|113|) is that the critical effective magnetic field 
h is smaller than the superfluid gap A gQp . If one had 
h > Ag a p the above scenario would not apply because 
the system would prefer to accommodate the excess po- 
larization by breaking pairs and creating quasi-particles. 
The gapless superfluid realized in this way would be ho- 
mogeneous and the transition to the normal state would 
be continous. Such a uniform phase is indeed expected 
to occur in the deep BEC regime (see below). 

The physical understanding of polarized Fermi gases 
became more complicated when exotic superfluid phases 
were proposed such as the inhomogeneous Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state (Fulde and Ferrell, 
1964; Larkin and Ovchinnikov, 1964). Other alternative 
states include the breached pair or Sarma state (Sarma, 
1963; Liu and Wilczek, 2003) and states with a deformed 
Fermi surface (Sedrakian et ai, 2005). In the FFLO 
state Cooper pairs carry a finite momentum resulting in 
a spontaneous breaking of translational symmetry with 
a periodic structure of the order parameter of the form 
A(x) oc cos(g:r), where the direction of the x axis is arbi- 
trary. The wavevector q is proportional to the difference 
of the two Fermi wavevectors q oc kp-\ — kp^ with a pro- 
portionality coefficient of order unity. The excess spin-up 
atoms are concentrated near the zeros of the order pa- 
rameter A(i). One can show that in the BCS limit the 
FFLO phase exists for h < 0.754A gap , corresponding to 
the tiny interval of polarization < P < 1.13A gap / Ep, 
and at P > 1.13A gap /Ep the system is normal (see, e.g., 
Takada and Izuyama, 1969). In the deep BCS regime 
this scenario is more energetically favourable compared 
to the Clogston-Chandrasekhar transition which would 
take place at P c = 1.06A gap /E F . The FFLO state 
is of interest both in condensed matter physics and in 
elementary particle physics, even though direct experi- 
mental evidences of this phase are still lacking [for re- 
cent reviews, see Casalbuoni and Nardulli (2004) and 
Combescot (2007)]. 

In ultracold gases the BCS regime is not easily achieved 
due to the smallness of the gap parameter and one is 
naturally led to explore configurations with larger val- 
ues of /c_F|a|, where the concept of Fermi surface looses 
its meaning due to the broadening produced by pair- 
ing. A major question is whether the FFLO phase sur- 
vives when correlations are stronger and if it can be real- 
ized in trapped configurations (Mizushima, Machida and 



Ichioka, 2005; Sheehy and Radzihovsky 2007; Yoshida 
and Yip, 2007). As a function of the interaction strength, 
parametrized by 1/kpa where kp = [37r 2 (n-f + tt-jJ] 1 / 3 is 
fixed by the total density, different scenarios can take 
place as schematically shown in Fig. [27] (Sheehy and 
Radzihovsky, 2006 and 2007; Son and Stephanov, 2006; 
Pao, wu and Yip, 2006; Hu and Liu, 2006; Iskin and 
Sa de Melo, 2006b; Parish et ai, 2007a). An important 
region of the phase diagram is the deep BEC regime of 
small and positive scattering lengths, where the energet- 
ically favourable phase consists of a uniform mixture of 
a superfluid gas of bosonic dimers and of a normal gas 
of spin polarized fermions. In this regime one expects 
that the normal uniform gas exists only for P = 1, cor- 
responding to the fully polarized ideal Fermi gas. The 
general problem of an interacting mixture of bosons and 
fermions was investigated by Viverit, Pethick and Smith 
(2000) who derived the conditions of miscibility in terms 
of the values of the densities and masses of the two com- 
ponents and of the boson-boson and boson-fermion scat- 
tering lengths (see Sec. HXEp . In particular, for P ~ 1 
corresponding to a small number of bosonic dimers in a 
fully polarized Fermi sea, the relevant condition for the 
solubility of the mixture reads [see Eq. (|130|) below] 



By using the values add = 0.60a and a a d = 1.18a for 
the dimer-dimer and atom-dimer scattering length re- 
spectively, one finds that the uniform phase exists for 
1/kpa > 2.1. This Bose-Fermi picture, however, looses 
its validity as one approaches the resonance region and 
more detailed analyses are needed to understand the 
phase diagram of the system close to unitarity (Pilati and 
Giorgini, 2007). In the presence of harmonic trapping the 
conditions of phase separation in the BEC regime change 
significantly due to the non-uniform effective potentials 
felt by the two components. The density profiles of the 
bosonic dimers and of the unpaired fermions have been 
investigated by Pieri and Strinati (2006) within the local 
density approximation. 

The determination of the energetically favourable con- 
figuration in the unitary regime of infinite scattering 
length is a difficult problem (Carlson and Reddy, 2005). 
Will the unpolarized superfluid and the polarized normal 
gas co-exist like in the BEC regime or will they separate? 
Will the FFLO phase play any role? First experiments 
carried out with spin imbalanced trapped Fermi gases 
close to unitarity suggest the occurrence of a phase sepa- 
ration between an unpolarized superfluid and a polarized 
normal phase (Partridge et ai, 2006a; Shin et ai, 2006). 

Other important questions that will not be addressed 
in this review concern the phases of these spin polarized 
Fermi gases at finite temperature (Gubbels, Romans and 
Stoof, 2006; Parish et ai, 2007a; Chien et ai, 2007) and 
the occurrence of p-wave superfluid phases at very low 
temperatures (Bulgac, Forbes and Schwenk, 2006). 
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B. Phase separation at unitarity 

A possible scenario for the equation of state of spin 
imbalanced configurations at unitarity is based on the 
phase separation between an unpolarized superfiuid 
and a polarized normal gas similar to the Clogston- 
Chandrasekhar transition discussed in the Sec. IIX Al (De 
Silva and Mueller, 2006a; Haque and Stoof 2006; Chevy 
2006a). An important ingredient of this scenario, that 
is not accounted for by the mean-field description, is 
the proper inclusion of interaction effects in the normal 
phase (Chevy 2006b; Lobo et al., 2006b; Bulgac and 
Forbes, 2007). While there is not at present a formal 
proof that the phase separated state is the most ener- 
getically favourable, the resulting predictions well agree 
with recent experimental findings (see Sec. IIX Cj) . 

The unpolarized superfiuid phase was described in de- 
tails in Sees. IIVI and W\ and is characterized, at unitarity, 
by the equation of state 

^- = \E F {l + f3), (117) 

where E$ is the energy of the system, N is the to- 
tal number of atoms and Ep = (h 2 /2m)(3ir 2 ns) 2 / 3 is 
the Fermi energy, where rig — 2n^ = In^ is the total 
density of the gas. Starting from (|117[) one can derive 
the pressure Ps = —dE/dV and the chemical potential 
Lis = dE s /dN. 

Differently from the superfiuid the normal phase is po- 
larized and consequently its equation of state will depend 



1.15 




FIG. 28: Equation of state of a normal Fermi g function 
of the concentration x. The solid line is a polynomial best fit 
to the QMC results (circles). The dashed line corresponds 
to the expansion (|119[) . The dot-dashed line is the coexis- 
tence line between the normal and the unpolarized superfiuid 
states and the arrow indicates the critical concentration x c 
above which the system phase separates. For x — 1, both the 
energy of the normal and of the superfiuid (diamond) states 
are shown. 



also on the concentration 

x = ni/ n 1 > ( 118 ) 

which, in the following, will be assumed to be smaller or 
equal to 1, corresponding to iV-j- > iVj_. A convenient way 
to build the x-dependence of the equation of state is to 
take the point of view of a dilute mixture where a few 
spin-down atoms are added to a non-interacting gas of 
spin-up particles. When i<1, the energy of the system 
can be written in the form (Lobo et al., 2006b) 

^^ = 3 -EMl-A, + ^ + ...], (119) 

A<j 5 to* 

where Ep-\ — (h 2 /2m)(6TT 2 n^) 2 / 3 is the Fermi energy of 
the spin-up particles. The first term in Eq. (|119[) corre- 
sponds to the energy per particle of the non-interacting 
gas, while the term linear in x gives the binding energy 
of the spin-down particles. Eq. (|119[) assumes that, when 
we add spin-down particles creating a small finite density 
ni , they form a Fermi gas of quasi-particles with effective 
mass to* occupying, at zero temperature, all the states 
with wavevector k up to km — {&'K 2 n\ i ) 1 / 3 ' and contribut- 
ing to the total energy (| 1 1 9[) with the quantum pressure 
term proportional to a; 5 / 3 . The interaction between the 
spin-down and spin-up particles is accounted for by the 
dimensionless parameters A and to/to*. The expansion 
(I119P should in principle include additional terms origi- 
nating from the interaction between quasi-particles and 
exhibiting a higher order dependence on x. 

The values of the coefficients entering (|119l) can be 
calculated using fixed-node diffusion Monte Carlo sim- 
ulations where one spin-down atom is added to a non- 
interacting Fermi gas of spin-up particles. The results of 
these calculations are A — 0.97(2) and to/to* = 1.04(3) 
(Lobo et al., 2006b). The same value of A has been ob- 
tained from an exact diagrammatic Monte Carlo calcula- 
tion (Prokof 'ev and Svistunov, 2007) and, quite remark- 
ably, also employing a simple variational approach based 
on a single particle-hole wavefunction (Chevy, 2006b, 
Combescot et al, 2007). The prediction of Eq. (TTTg|) 
for the equation of state is reported in Fig. [2H1 where we 
also show the FN-DMC results obtained for finite values 
of the concentration x. It is remarkable to see that the 
expansion (|1 19[) reproduces quite well the best fit to the 
FN-DMC results up to the large values of x where the 
transition to the superfiuid phase takes place (see discus- 
sion below). In particular, the repulsive term in x 5 ^ 3 , 
associated with the Fermi quantum pressure of the mi- 
nority species, plays a crucial role in determining the x 
dependence of the equation of state. Notice that when 
x = 1 (JVf = Ni) the energy per particle of the nor- 
mal phase is smaller than the ideal gas value 6Ep^/5, 
reflecting the attractive nature of the force, but larger 
than the value in the superfiuid phase (QA2)6Ep^ /5 (see 
Sec. lVBj) . 

We can now determine the conditions of equilibrium 
between the normal and the superfiuid phase. A first 
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condition is obtained by imposing that the pressures of 
the two phases be equal 



dE s dE. 



dV s dVi 



N 



(120) 



where Vs and Vn are the volumes occupied by the two 
phases respectively. A second condition is obtained by 
requiring that the chemical potential of each pair of spin 
up-spin down particles be the same in the two phases. In 
order to exploit this latter condition one takes advantage 
of the thermodynamic identity /is = (/if +/xj.)/2 for the 
chemical potential in the superfluid phase yielding the 
additional relation 



8E S _ 1 fdEN_ + 9En 
OX 1 \ ON 



ON 



(121) 



where we have used the expression Mt(J.) = 9En /dN^n\ 
for the chemical potentials of the spin-up and spin-down 
particles calculated in the normal phase. Eq. (|121[) . com- 
bined with (|120p . permits to determine the critical val- 
ues of the thermodynamic parameters characterizing the 
equilibrium between the two phases. For example, if 
applied to the BCS regime where En = 3/5(NiEf^ + 
NtEpt) and E s = E N ~3NA 2 gap /8E F [see Eq. (|56]>]. this 
approach reproduces the Clogston-Chandrasekhar condi- 
tion (|113|) (Bcdaque, Caldas and Rupak, 2003). At uni- 
tarity instead, a calculation based on the QMC values of 
En and Es yields the values x c = 0.44, corresponding 
to P c = (1 - 0/(1 + x c ) = 0.39, and (n N /n s ) c = 0.73, 
where tin = + is the density of the normal phase 
in equilibrium with the superfluid (Lobo et at, 2006b). 
For values of the polarization larger than P c = 0.39 the 
stable configuration is the uniform normal phase, while if 
we increase the number of the spin down particles (corre- 
sponding to a reduction of P) there will be a phase sep- 
aration between a normal phase with the concentration 
x c = 0.44 and a superfluid unpolarized phase. The phase 
transition has first order character consistently with the 
critical value h — 0.8lA gap being smaller than the su- 
perfluid gap. In particular, while the spin-up density 
is practically continuous, the density of the spin-down 
particles exhibits a significant jump at the transition. 
Let us finally point out that the parameters characteriz- 
ing the transition between the superfluid and the normal 
phases depend in a crucial way on the many-body scheme 
employed for the calculation. For example, if instead 
of the Monte Carlo results we use the BCS mean-field 
theory of Sec. IV Al and the non-interacting expression 
En = 3/5(N-[Epi + NiEpi) for the energy of the normal 
phase, we would find the very different value x c — 0.04 
for the critical concentration. 



trapping. We make use of the local density approxi- 
mation which permits to express the local value of the 
chemical potential of each spin species as 



MT(l)( r ) =M?(|) - Vh (r) 



(122) 



The chemical potentials are fixed by imposing the 

proper normalization to the spin-up and spin-down den- 
sities and the condition M^m — V/ lo (r = R^{i)) = de- 
fines the Thomas-Fermi radii R](\) of the two species. In 
our discussion we assume isotropic trapping in order to 
simplify the formalism. A simple scaling transformation 
permits to apply the results also to the case of anisotropic 
trapping. 

For small concentrations of the spin-down particles 
(N± <C N{) only the normal state is present and one 
can ignore the change in /i-f due to the attraction of 
the spin-down atoms. In this case n-\ reduces to the 
Thomas-Fermi profile ([8]) of an ideal gas, whereas nj, is a 
Thomas-Fermi profile with a modified harmonic poten- 
tial Vho — ► (1 + 3A/5)Vho- The confining potential felt 
by the spin-down atoms is stronger due to the attraction 
produced by the spin-up atoms. One can also calculate 
the frequency of the dipole oscillation of the spin-down 
particles. This is given by (i — x, y, z) 



dipole 



3 m 
(1 + -A) — ~ 1.23^ 

5 m 



(123) 



and is affected by both the interaction parameter A and 
the effective mass parameter m/m* entering the expan- 
sion pro]) . 

When in the center of the trap the local concentration 
of spin-down particles reaches the critical value x c — 0.44 
(see Sec. IPC Bp a superfluid core starts to nucleate in 
equilibrium with a polarized normal gas outside the su- 
perfluid. The radius Rsf of the superfluid is determined 
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C. Phase separation in harmonic traps at unitarity 

The results presented in Sec. IPC Bl can be used to cal- 
culate the density profiles in the presence of harmonic 



FIG. 29: Radii of the three phases in the trap in units of the 
radius itf = a ho (48iV T ) 1/6 of a non-interacting fully polarized 
gas. 
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by the equilibrium conditions between the two phases dis- 
cussed in the previous Section. Towards the periphery of 
the polarized normal phase the density of the spin-down 
particles will eventually approach a vanishing value cor- 
responding to R±. For even larger values of the radial 
coordinate only the density of the spin-up particles will 
be different from zero up to the radius In this pe- 
ripherical region the normal phase corresponds to a fully 
polarized, non-interacting Fermi gas. 

By using the Monte Carlo equation of state discussed 
in Sec. lLXBl one finds that the superfluid core disappears 
for polarizations P > P c tra P = 0.77 (Lobo et al, 2006b). 
The difference between the critical value P*™p and the 
value P c = 0.39 obtained for uniform gases reflects the 
inhomogeneity of the trapping potential. In Fig. [55] we 
show the calculated radii of the superfluid and of the spin- 
down and spin-up components as a function of P. The 
radii are given in units of the non-interacting radius of 
the majority component = a? lo (487V|) 1 ' 6 . The figure 
explicitly points out the relevant features of the prob- 
lem. When P — > one approaches the standard unpo- 
larized superfluid phase where the three radii (Rsf, R] 
and i? x ) coincide with the value (1 + (3)^ 4 R^ ~ 0.81^ 
[see Eq. (|74|) ]. By increasing P one observes the typical 
shell structure with R$ f < R\ < R*\ ■ The spin-up radius 
increases and eventually approaches the non-interacting 
value when P — > 1. The spin-down radius instead de- 
creases and eventually vanishes as P — ► 1. Finally the 
radius Rsf of the superfluid component decreases and 
vanishes for P* ra P = 0.77, corresponding to the disap- 
pearance of the superfluid phase. 

The predicted value for the critical polarization is in 
good agreement with the findings of the MIT experiments 
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FIG. 30: Phase separation in an interacting Fermi gas at 
unitarity. Normalized central density difference (circles) and 
condensate fraction (triangles). From both sets of data one 
extract the same value ~ 0.75 for the critical polarization. 
From Shin et al. (2006). 




FIG. 31: Double integrated density difference measured at 
unitarity in an ultracold trapped Fermi gas of 6 Li with po- 
larization P — 0.58 (from Shin et al, 2006). The theoretical 
curves correspond to the theory of Sees. IIX BlIlXCl based on 
the local density approximation (red solid line) and to the 
predictions of a non-interacting gas with the same value of P 
(dashed line). (From Lobo, Recati and Stringari, 2007). 



(Zwierlein et al., 2006a; Shin et al., 2006), where the 
interplay between the superfluid and the normal phase 
was investigated by varying the polarization P of the 
gas. The experimental evidence for superfluidity in a 
spin imbalanced gas emerges from measurements of the 
condensate fraction and of the vortex structure in fast 
rotating configurations (Zwierlein et al., 2006a). These 
time-of-flight measurements were performed by rapidly 
ramping the scattering length to small and positive val- 
ues after opening the trap, in order to stabilize fermion 
pairs and increase the visibility of the bimodal distribu- 
tion and of vortices (see discussion in Sec. IV B 21 and in 
Sec. IVIIICp . In another paper (Shin et al., 2006) the 
in situ density difference n-\ (r) — (r) was directly mea- 
sured with phase contrast techniques. Phase separation 
was observed by correlating the presence of a core region 
with n-f — n j = with the presence of a condensate of 
pairs (see Fig. |3"0|) . At unitarity these results reveal that 
the superfluid core appears for P < 0.75. An interest- 
ing quantity that can be directly extracted from these 
measurements is the double integrated density difference 

n ( d 1) (z) = J dxdy[n^{r) -n l (r)\ , (124) 

which is reported in Fig. [31] for a unitary gas with 
P = 0.58. The figure reveals the occurrence of a char- 
acteristic central region where n^\z) is constant. The 
physical origin of this plateau can be understood using 
the local density approximation (De Silva and Mueller, 
2006a; Haque and Stoof, 2006): it is a consequence of 
the existence of a core region with (r) = n^(r) which 
is naturally interpreted as the superfluid core. Further- 
more, the value of z where the density exhibits the cusp 
corresponds to the Thomas- Fermi radius Rsf of the su- 
perfluid. In the same figure we show the theoretical pre- 
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dictions, based on the Monte Carlo results for the equa- 
tion of state of the superfluid and normal phases, which 
well agree with the experimental data (Lobo, Recati and 
Stringari, 2007). 

The theoretical predictions discussed above are based 
on a zero temperature assumption and on the local den- 
sity approximation applied to the various phases of the 
trapped Fermi gas. While the applicability of the LDA 
seems to be adequate to describe the MIT data, the Rice 
experiments (Partridge et al., 2006a), carried out with a 
very elongated trap, show that in this case surface tension 
effects, not accounted for by LDA, play a major role. It 
is possible to prove (De Silva and Mueller, 2006a; Haque 
and Stoof, 2006) that the double integrated density dif- 
ference, when evaluated within the LDA, should exhibit 
a monotonic non increasing behavior as one moves from 
the trap center. The non monotonic structure of the den- 
sity observed in the Rice experiment can be explained 
through the inclusion of surface tension effects. For a re- 
cent discussion of surface tension and thermal effects in 
spin imbalanced configurations see De Silva and Mueller 
(2006b), Partridge et al. (2006b) and Haque and Stoof 
(2007). 

Other important questions concern the dynamic be- 
havior of these spin polarized configurations. The prob- 
lem is interesting since the dynamic behavior of the su- 
perfluid is quite different from the one of the normal gas, 
the latter being likely governed, at very low temperature, 
by collisionlcss kinetic equations rather than by the equa- 
tions of hydrodynamics. As a consequence the role of the 
boundary separating the two phases should be carefully 
taken into account for a reliable prediction of the dy- 
namic properties, such as the collective oscillations and 
the expansion. 



D. Fermi superfluids with unequal masses 

When one considers mixtures of Fermi gases belong- 
ing to different atomic species and hence having differ- 
ent masses new interesting issues emerge. One should 
point out that even if the masses are different, config- 
urations where the atomic densities of the two compo- 
nents are equal result in equal Fermi momenta: fcpf = 
hp i = (37T 2 ?!) 1 / 3 , where n = nj + n± is the total den- 
sity. This means that the mechanism of Cooper pairing, 
where two atoms of different spin can couple to form a 
pair of zero momentum, is still valid. At T = the 
BCS mean-field theory predicts a simple scaling behav- 
ior of the equation of state in terms of the reduced mass 
of the two atoms, which holds in the whole BCS-BEC 
crossover. However, in the BCS regime, correlations be- 
yond the mean-field approximation introduce a non triv- 
ial dependence on the mass ratio in the superfluid gap 
(Baranov, Lobo and Shlyapnikov, 2007). Furthermore, in 
the BEC regime, theoretical studies of four-fermion sys- 
tems (Petrov, 2003; Petrov, Salomon and Shlyapnikov, 
2005) emphasize the crucial role of the mass ratio on the 



interaction between dimers, resulting in instabilities if 
the mass ratio exceeds a critical value (see Sec. MI CI) . 
First quantum Monte Carlo results have also become 
available in the crossover, exploring the dependence of 
the equation of state and of the superfluid gap on the 
mass ratio (von Stecher, Greene and Blume 2007; As- 
trakharchik, Blume and Giorgini, 2007). These studies 
point out significant deviations from the predictions of 
the BCS mean-field theory. Important questions emerge 
also at finite temperature, where one can immediately 
conclude that the simple scaling in terms of the reduced 
mass can not hold since, for example, the BEC transition 
temperature of composite bosons (see Sec. IIVE[) should 
depend on the molecular mass, M — m-\ + mi, rather 
than on the reduced mass. Other interesting scenarios, 
not addressed in this review, refer to the interplay be- 
tween unequal masses and unequal populations of the 
two spin components (see, e.g., Wu, Pao and Yip, 2006; 
Iskin and Sa de Melo, 2006b; Parish et al, 2007b; Iskin 
and Sa de Melo, 2007a). The experimental realization 
of fermionic heteronuclear mixtures in the strongly cor- 
related regime is at present actively pursued in different 
laboratories (see, e.g., Wille et al, 2007). 



Equation of state along the crossover 



The BCS mean-field approach developed in Sec. IV Al 
can be generalized in a straightforward manner to the 
case of unequal masses. Starting from the BCS Hamilto- 
nian (|46[) . where the single-particle energy term contains 
now the mass WjQ) and the chemical potential fJ-fil) °f 
the two spin components, one follows the same steps lead- 
ing to the gap and number equations (|5"4"1) - ([5"5"|) . These 
equations, as well as Eqs. (|5"0"|) - (f5"Tj) for the quasi-particle 
amplitudes and excitation energies, read exactly the same 
as in the equal mass case, with the only difference that 
the atomic mass m should be replaced by twice the re- 
duced mass 



2m r = 2- 



(125) 



and that the single-particle energy should be replaced by 
r/k — ► h 2 k 2 /Am r — li, where [i = (/ij+/i|)/2 is the average 
chemical potential. In units of the reduced Fermi energy 
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h 2 k F 



(126) 



the resulting values of the chemical potential /i and of 
the order parameter A, for a given value 1/kpa of the 
interaction strength, are then independent of the mass 
ratio mj /wj. 

At unitarity dimensionality arguments permit to write 
the energy of the system in the general form E$/N = 
3/5 [1 + 0{m^/m i ))E F , where N = iV T + N t is the to- 
tal number of particles. In the BCS mean-field approach 
the value of f3 is given by = —0.41 (see Table [TT] in 
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Sec lVIl) . First results based on QMC simulations, sug- 
gest that the dependence of (3 on the mass ratio is very 
weak (Astrakharchik, Blume and Giorgini, 2007). On the 
contrary, the corresponding value of the superfluid gap 
A gap is strongly reduced by increasing the mass ratio. 

In the BEC regime, the mean-field approach predicts 
that the binding energy of the dimers is given by the 
formula ef, = — h 2 /2m r a 2 and that these molecules inter- 
act with the same scattering length as in the symmetric 
mass case (add = 2a). While the result for the binding 
energy is correct, the actual relationship between add and 
a is sensitive to the value of the mass ratio, approaching 
the value 0.60 when — m x (Petrov, 2003; Petrov, 
Salomon and Shlyapnikov, 2005). 

2. Density profiles and collective oscillations 

Also in the case of unequal masses the equation of state 
of uniform systems can be used to evaluate the density 
profiles of the harmonically trapped configurations. In 
the local density approximation the chemical potential 
of each atomic species varies in space according to the 
law (|122|) where, however, the confining potential is now 
spin dependent being related to different magnetic and 
optical properties of the two atomic species and should 
be replaced by V ho (r) -> V^ l) (r). Since the chemical 
potential of the superfluid phase is given by the average 
/is = (a*T + Mi)/2> the corresponding density is deter- 
mined at equilibrium by the Thomas-Fermi relation 

M») = Mo - V ho (r) , (127) 

where 

M 

V ho (r) = — (£o 2 x x 2 + Q 2 y y 2 + Co 2 z z 2 ) , (128) 

is the average of the trapping potentials. The effective 
frequencies are given by 

- 2 = m T K T ) 2 + m x K X ) 2 (129) 
mi + mj, 

where u>J are the oscillator frequencies relative to 

, whereas M = m-f + denotes the mass of the 
pair. If the oscillator lengths of the two atomic com- 
ponents coincide, i.e. if m^ojj — mj,w|, the effective 

frequencies (|129|) take the simplified form cDj = J u>Ju>^. 

In the superfluid phase the densities of the spin-up 
and spin-down atoms are equal, even if the trapping po- 
tentials, in the absence of interatomic forces, would give 
rise to different profiles. In principle, even if one of the 
two atomic species does not feel directly any external 
potential (for example u>i =0), it can be trapped due 
to the interaction with the other species. At unitar- 
ity, one can show (Orso, Pitaevskii and Stringari, 2007) 
that the trapped superfluid configuration is energetically 



favourable if the condition (1 + /3)M/mi < 1 is satisfied. 
This condition is easily fulfilled if mj < m^. 

The dynamic behavior of Fermi superfluids with un- 
equal masses can be studied by properly generalizing the 
equations of hydrodynamics which take the form (|54"|) - 
([55)1 with Vho replaced by the effective trapping poten- 
tial (|128[) and where m is replaced by half of the pair 
mass m — > M/2. In uniform systems (V^ c = 0) the 
equations give rise to the propagation of sound with 
the sound velocity fixed by the thermodynamic relation 
Mc 2 = 2nsdjjLs I 'dns- In the BCS limit, where the equa- 
tion of state approaches the ideal gas expression the 
sound velocity takes the value c = Hkp (Sm^mi) 
where hp — (ii^n) 1 ^. At unitarity the sound veloc- 
ity is given by the above ideal gas value multiplied by 
the factor yT + f3. 

For harmonically trapped configurations the hydrody- 
namic equations can be used to study the expansion of 
the gas after opening the trap as well as the collective 
oscillations (Orso, Pitaevskii and Stringari, 2007). The 
effect of the mass asymmetry is accounted for through the 
effective frequencies (II 29[) as well as through the changes 
in the equation of state /xg. At unitarity, where the den- 
sity dependence of the chemical potential keeps the n 2 / 3 
power law, the same results discussed in Sec. I VII Cl hold 
with the oscillator frequencies replaced by £>,•. Notice 
that the effective frequencies cD, differ in general from ei- 
ther uij and ojj. Actually in the superfluid phase the 
two atomic species can not oscillate independently, but 
always move in phase as a consequence of the pairing 
mechanism produced by the interaction. 

E. Fermi-Bose mixtures 

The problem of quantum degenerate mixtures of a spin 
polarized Fermi gas and a Bose gas has been the object 
of considerable experimental and theoretical work. Par- 
ticularly interesting are the mixtures where the interac- 
tions between the fermions and the bosons are tunable 
by means of a Feshbach resonance. This is the case, for 
example, of the 40 K- 87 Rb system which has been exten- 
sively investigated by the groups in Florence (Modugno et 
at, 2002; Ferlaino et at, 2006; Zaccanti et al, 2006; Mod- 
ugno, 2007) and in Hamburg (Ospelkaus et al, 2006a; 
Ospelkaus et al., 2006b). In these configurations the mix- 
ture is characterized by a repulsive boson-boson interac- 
tion with scattering length aBB — 5.3nm and a magneti- 
cally tunable boson-fermion interaction parametrized by 
the scattering length asF- Both the observation of an 
induced collapse of the mixture for large negative values 
of as f and of phase separation for large positive values 
of obf have been reported. 

From the theoretical point of view these experimental 
findings can be understood using a mean-field approach. 
The conditions of stability of a Bose-Fermi mixture in 
uniform systems at T = have been investigated by 
Viverit, Pethick and Smith (2000). Depending on the 
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sign of cibf different scenarios can apply. 

If a,BF < 0, the only relevant condition is determined 
by mechanical stability, i. e. by the requirement that the 
energy of the mixture must increase for small fluctua- 
tions in the density of the two components. Starting 
from a mean-field energy functional, which includes to 
lowest order the interaction effects between bosons and 
between fermions and bosons, the linear stability require- 
ment fixes an upper limit on the fermionic density nj? 
irrespective of the value of the bosonic density 
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aBF\ < — 
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(130) 



where tub and mp denote, respectively, the masses of 
bosons and fermions and ass is assumed to be positive 
to ensure the stability of the configurations where only 
bosons are present. If Uf exceeds the upper bound (|130|) 
the system collapses. 

If instead clbf > 0, the uniform mixture can become 
unstable against phase separation into a mixed phase 
and a purely fermionic one. While the fermionic density 
can not in any case violate the condition (|130p . yet for 
each value of np there exists a critical bosonic density 
n c B (np) above which the system phase separates. The 
function n c B (np) is nontrivial and for a given pair of den- 
sities {np and ns) one has to check whether the phase 
separated configuration is in equilibrium and is energet- 
ically favourable compared to the uniform mixture. For 
vanishingly small bosonic densities Ub, the relevant con- 
dition for np coincides with (|130|) . It is worth noticing 
that the repulsive as f scenario describes the mixtures of 
fermions and composite bosons in the deep BEC regime 
considered in Sec. IPC At 

In the presence of harmonic trapping the conditions for 
the collapse and for the phase separation change. The 
density profiles of the two components have been investi- 
gated by M0lmer (1998) using the same mean-field energy 
functional described above in the local density approxi- 
mation. The results are consistent with the scenario of 
a collapsed state if olbf is large and negative and of a 
phase separated state (a core of bosons surrounded by 
a shell of fermions) in the opposite regime of large and 
positive cibf- These two scenarios are in agreement with 
the features observed in 40 K- 87 Rb mixtures (Ospelkaus 
et ai, 2006b; Zaccanti et al, 2006; Modugno, 2007). 

The collective oscillations in harmonically trapped 
Bose-Fermi mixtures have also been the object of a 
considerable number of theoretical investigations (see, 
e.g., the recent work by Maruyama and Bertsch, 2007, 
and references therein). In particular, the monopole 
(Maruyama, Yabu and Suzuki, 2005) and the dipole 
mode (Maruyama and Bertsch, 2007) have been studied 
at zero temperature using a dynamic approach based on 
the solution of a coupled system of time-dependent equa- 
tions: the Gross-Pitaevskii equation for the bosons and 
the collisionless Vlasov equation for the fermions. These 
studies point out the existence of a characteristic damp- 



ing in the motion of the fermionic component affecting 
both types of oscillations. 

An important aspect of Bose-Fermi mixtures concerns 
the boson-induced interactions experienced by the other- 
wise non- interacting fermionic atoms. The physical ori- 
gin of the induced interactions is the polarization of the 
bosonic medium which acts as an effective potential be- 
tween the fermions. The density-density response func- 
tion of the bosons is the relevant quantity to describe this 
effect and the induced interaction is thus frequency and 
wavevector dependent. At low frequencies and long wave- 
lenghts it is always attractive, irrespective of the sign of 
asF, is independent of the density of bosons and repro- 
duces the mechanism of instability discussed in the para- 
graph after Eq. (fT3T))) (Bijlsma, Heringa and Stoof, 2000; 
Viverit, Pethick and Smith, 2000). The physical picture 
is similar to the effective attraction between 3 He atoms 
in solution in superfluid 4 He (Edwards et ai, 1965) and 
to the most famous phonon-induced attraction between 
electrons in ordinary superconductors (see, for example, 
de Gennes, 1989). 

An additional major interest of exploiting heteronu- 
clear Feshbach resonances (including the Fermi-Fermi 
mixtures considered in Sec. ILXDp is the possibility of 
creating polar molecules characterized by long-range 
anisotropic interactions which are expected to have a pro- 
found impact on the many-body physics (Baranov et al, 
2002). Furthermore, the creation of such heteronuclear 
molecules in optical lattices (Ospelkaus et ai, 2006c) 
might lead to important applications in quantum infor- 
mation processing (Micheli, Brennen and Zoller, 2006). 



X. FERMI GASES IN OPTICAL LATTICES 

The availability of optical lattices has opened new fron- 
tiers of research in the physics of ultracold atomic gases 
(for a recent review, see Morsch and Oberthaler, 2006, 
and Bloch, Dalibard and Zwerger, 2007). In the case of 
Fermi gases the problem is closely related to the physics 
of electrons in metals or semiconductors. However, op- 
tical lattices differ favorably from traditional crystals in 
many important aspects. The period of the optical lattice 
is macroscopically large, which simplifies the experimen- 
tal observation. The lattice can be switched off, and its 
intensity can be tuned at wish. Atoms, differently from 
electrons, are neutral and furthermore their interaction is 
tunable thanks to the existence of Feshbach resonances. 
The lattices are static and practically perfect, being free 
of defects. Disorder can be added in a controllable way by 
including random components in the optical field. Lastly, 
it is easy to produce one- and two-dimensional structures. 

An atom in a monochromatic electric field feels a time 
averaged potential proportional to the square of the field 
amplitude (see, e.g., Pethick and Smith, 2002; Pitaevskii 
and Stringari 2003) . It is useful to work sufficiently close 
to the frequency loq of the absorption line of an atom, 
where the force on the atom becomes strong. The atoms 
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are pulled into the strong- field region for lu < ujq ("red 
detuning") and pushed out of it for uj > uj ("blue de- 
tuning" ) . 

One-dimensional periodic potentials can be produced 
by a standing light wave. In this case, the potential en- 
ergy is conveniently written as V op t (z) — S-Ersui 2 (Kz), 
where K is the wavevector of the laser, Er — h 2 K 2 /2m is 
the "recoil energy" and s is the dimensionless parameter 
proportional to the laser field intensity. Typical values 
of s in experiments range from 1 to 20. The potential 
has a period d = A/2 = tt/ K , where A is the wavelength 
of the laser. If the two counterpropagating laser beams 
interfere under an angle less than 180° the period is 
increased by the factor sin.(0/2) -1 . By using three mutu- 
ally orthogonal laser beams one can generate a potential 
of the form 

V opt (r) - sE R [sin 2 (Kx) + sin 2 (Ky) + sin 2 (Kz)] . 

(131) 

Notice that in experiments atoms are trapped by ad- 
ditional confining potentials, in most cases of harmonic 
form. To preserve the effective periodicity of the problem 
the harmonic potential should vary slowly with respect 
to the period of the lattice. This demands the condition 
Er 3> hujho, which is easily satisfied in experiments. 



A. Ideal Fermi gases in optical lattices 

1. Fermi surface and momentum distribution 

An advantage of experiments with cold fermions is the 
possibility of realizing a non-interacting gas by creating a 
spin-polarized configuration. As explained in Sec. IIII Al 
the interaction between fermionic atoms with parallel 
spins is in fact negligible. 

The quantum mechanical description of the motion of 
a particle in a periodic external field was developed by 
Bloch (Bloch, 1928). In ID the wavefunctions have the 
Bloch form i\) q ^ n (z) = exp(iq z z)u qzn (z) and are classified 
in terms of the quasi- momentum p z = Hq z , while n is 
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FIG. 32: Time-of-flight images obtained after adiabatically 
ramping down the optical lattice. Image (a) is obtained with 
A CT = 3500 and s = 5Er. Images (b)-(e) are obtained with 
N a — 15000 and correspond to s = 5Er (b), s = 6En (c), 
s = 8Er (d) and s = V2Er (e). The images show the optical 
density (OD) integrated along the vertically oriented z axis 
after 9 ms of ballistic expansion. From Kohl et al. (2005). 



a discrete number labelling the Bloch band. The values 
of the wavevector q z differing by the reciprocal lattice 
vector 2K are physically equivalent and it is therefore 
enough to restrict the values of q z to the first Brillouin 
zone: —K < q z < K . The function u q , n (z) is a periodic 
function of z with period d. The formalism is straightfor- 
wardly generalized to 3D lattices where the eigenstates 
of the Hamiltonian with the potential (|131|) are prod- 
ucts of wavefunctions of the Bloch form along the three 
directions and are hence classified in terms of the 3D 
quasi-momcntum p. In a Fermi gas at zero temperaure 
all the states with excitation energy e(p) = E(p) — E(0) 
are occupied up to values such that e(p) = Ep, where 
Ep is the Fermi energy. The corresponding values of p 
characterize the Fermi surface. 

Experimentally one can measure the momentum distri- 
bution of these non-interacting configurations by imag- 
ing the atomic cloud after release from the trapping po- 
tential. In fact, the spatial distribution n(r) of a non- 
interacting expanding gas reproduces asymptotically the 
initial momentum distribution n(p) according to the law 
n(r,t) — > (m/t) 3 n(p = mr/t). 

In order to have access to the quasi-momentum distri- 
bution a practical procedure consists of switching off the 
lattice potential in an adiabatic way so that each state 
in the lowest energy band with quasi-momcntum p is 
adiabatically transformed into a state with momentum 
p. The condition of adiabaticity requires that the lattice 
potential be switched off in times longer than the inverse 
of the energy gap between the first and second band. 

The Bloch problem can be solved analytically in the 
tight binding approximation, holding for large lattice 
heights. For a 3D cubic lattice one finds the result 

(p = ftq) 

e p = 25 [sin 2 (fed/2) + sin 2 (q y d/2) + sin 2 (q z d/2)](132) 

for the dispersion law of the single-particle excita- 
tions in the lowest band. The band width 25 de- 
creases exponentially for large s according to 5 = 
8E RS V 4 exp(-2^)/^ (Zwerger, 2003), and is in- 
versely proportional to the tunneling rate through the 
barriers. The energy gap between the first and second 
band coincides with the energy splitting huj opt = 2^/sEr 
between the states in the harmonic potential produced 
by the optical potential (|131[) around each local mini- 
mum. When p — * 0, Eq. (|132p takes the simple form 
e p = p 2 /2m*, the effective mass being related to the 
band width by the relation m* = h 2 /5d 2 . The explicit 
dependence of m* on s, including the case of small laser 
intensities, was calculated by Kramer et al. (2003). 

The dispersion law (|132|) results from the removal of 
the degeneracy between the lowest energy states of each 
well produced by the intcrwcll tunneling. This moans 
that the number of levels in the band is equal to the 
number of the lattice cells. Then, if one works within 
the first Bloch band, the maximum achievable density is 

n max = l/ff. 

The inclusion of harmonic trapping can be accounted 
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for by introducing the semiclassical distribution function 
which, at zero temperature, takes the form / (p, r) = 
O [Ep — H (p, r)] where Ep is the Fermi energy. In the 
above equation p is the quasi-momentum variable and 
H = e p + Vfjo(r), with e p given by (|132|) . Starting from 
the distribution function, one can evaluate the quasi- 
momentum (qm) distribution by integrating over r: 

where u>ho "C ui op t is the usual geometrical average of 
the harmonic frequencies of the potential Vto (r) and p is 
restricted to the first Brillouin zone. 

If the Fermi energy is much smaller than the band 
width 25, one can expand the dispersion law up to terms 
quadratic in p. In this case one recovers the same 
Thomas-Fermi form (|10[) holding for the momentum dis- 
tribution in the absence of the optical potential, the only 
difference being the presence of an effective-mass term 
which renormalizes the trapping frequencies. In the op- 
posite regime Ep ^> 25, but still Ep < fiw opt , the quasi- 
momentum distribution becomes flat within the first Bril- 
louin zone, giving rise to a characteristic cubic shape for 
the Fermi surface. In this limit we find 

The experimental investigation of the Fermi surface in a 
3D optical lattice was carried out by Kohl et al. (2005), 
who observed the transition from the spherical to the cu- 
bic shape by increasing the intensity of the laser generat- 
ing the optical lattice (see Fig. [321) ■ F° r s — \2 [panel (e) 
in the figure], corresponding to 5 ~ 10 nK, the values of 
the relevant parameters were Er = 348 nK, fkt)ho = 2~kYl 
191 Hz=9.2 nK and N a = 15000, so that the conditions 
hu) pt S> Ep ^> 25, needed to reach the cubic shape for 
the Fermi surface, were well satisfied in this experiment. 

In the same limit Ep 3> 25 the coarse-grained density 
distribution takes the constant value n™ ax — 1 /d 3 within 
the ellipsoid fixed by the radii (i = x, y, z) 

R i =(±y\ys d ^. (135) 

where d is the periodicity of the laser field. The constant 
value of the density reflects the insulating nature of the 
system. 

2. Block oscillations 

Atomic gases confined by optical lattices are well suited 
to study Bloch oscillations (these oscillations are difficult 
to observe in natural crystals because of the scattering 
of electrons by the lattice defects). Let us consider a 
ID optical lattice in the presence of a gravitational field 



V ex t — —mgz oriented along the direction of the lat- 
tice. The dynamics of atoms in the lowest band can be 
described in the semiclassical approximation using the 
effective single-particle Hamiltonian 

H (p, r) = + e + Vext (r) . (136) 
2m 

The Hamilton equation dp z /dt = —dH/dz then yields 
the obvious solution p z = mgt. In this problem it is con- 
venient not to restrict p z to the first Brillouin zone, but 
to allow p z to take also larger values. Then all the ob- 
servable physical quantities must be periodic functions 
of p z with period 2KH. From the time dependence of p z 
it follows that these quantities oscillate in time with the 
frequency lob — mgd/h. The periodicity of these Bloch 
oscillations is ensured by the periodicity of the optical 
lattice and the theory is applicable if <C u> pt- 

Bloch oscillations have been observed in Bose gases 
both above (Ben Dahan et al., 1996; Clade et al., 2006; 
Ferrari et al., 2006) and below (Anderson and Kasevich, 
1998; Morsch et al., 2001) the critical temperature for 
Bose-Einstein condensation. A general limitation in the 
study of Bloch oscillations with bosons is due to instabil- 
ities and damping effects produced by the interactions. 
Precise measurements have been recently achieved in a 
dense BEC gas of 39 K, near a Feshbach resonance which 
permits to tune the scattering length to vanishing values 
(Roati et al, 2007). 

The use of polarized fermions permits to work with 
relatively dense gases because of the absence of s-wave 
collisions. In the experiment of Roati et al. (2004) a 
Fermi gas of 40 K atoms was initially confined in a har- 
monic trap so that the external potential takes the form 
Vext(r) = V^ (r) — mgz. The quasi-momentum distri- 
bution is obtained by integrating the T = distribution 
function in all the variables except p z : 

n^ m \ Pz ) oc (Ep - e p J 5/2 Q(E F - e p J . (137) 

As long as Ep is smaller than 25 the quasi-momentum 
distribution is localized in a narrow region around p z = 0, 
while the contrast deteriorates for larger values of Ep . 

At t = the vertical harmonic confinement is sud- 
denly switched off and the atoms evolve in the presence 
of the lattice and gravitational potentials. At the ini- 
tial time the quasi-momentum distribution is centered 
at p z = 0. It will later move according to the law 
n {i m ) (p z ,t) — n^ qm ^ (p z — mgt). When the cloud reaches 
the edge of the Brillouin zone K, it reappears on the 
opposite side and the quasi-momentum distribution ac- 
quires a two-peak character. At t = 2t: /lub it regains its 
initial shape. 

After a given evolution time the lattice potential is 
adiabatically switched off in order to transfer the quasi- 
momentum distribution into the momentum one. The 
cloud is then imaged after a given time of free expan- 
sion. In this experiment it was possible to observe about 
100 Bloch periods. The high precision achievable in the 
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measurement of Bloch frequencies opens new perspec- 
tives in sensitive measurements of weak forces, like the 
Casimir-Polder force between atoms and a solid substrate 
(Carusotto et al, 2005). 



3. Center of mass oscillation 

In addition to Bloch oscillations it is also interesting 
to study the consequence of the lattice on the oscillations 
of the gas occurring in coordinate space. In this Section 
we focus on the dipole oscillation which can be excited 
by a sudden shift of the confining harmonic trap. Ac- 
cording to Kohn's theorem the dipole oscillation in the 
absence of the lattice has no damping and its frequency 
is equal to the trap frequency. In a superfluid these os- 
cillations exist also in the presence of the lattice thanks 
to the coherent tunneling of atoms through the barriers 
separating consecutive wells. The main consequence of 
the lattice is a renormalization of the collective frequency 
determined by the effective mass of the superfluid. These 
oscillations have already been observed in Bose-Einstein 
condensates (Cataliotti et al., 2001) and are expected to 
occur also in Fermi superfluids (Pitaevskii, Stringari and 
Orso, 2005). 

The behavior of a non-interacting Fermi gas is very 
different (Pezze et at, 2004). To understand the origin of 
the differences let us consider the simplest case of a one- 
dimensional Fermi gas characterized by the dispersion 
law e Pz — 26 sin 2 (p z d/2H) and trapped by the harmonic 
potential mui^z 2 /2. Atoms with energy smaller than 26 
can perform closed orbits in the z—p z phase plane. These 
atoms oscillate around the center of the trap. Vice-versa, 
atoms with energy higher than 26 perform open orbits, 
being unable to fully transfer the potential energy into 



the Bloch energy e Pz . They consequently perform small 
oscillations in space, remaining localized on one side of 
the harmonic potential. As a consequence if Ep > 26 
the cloud no longer oscillates around the new center of 
the trap but is trapped out of the center and performes 
small oscillations around an offset point, reflecting the 
insulating nature of the system. 

In order to investigate a three-dimensional case, one 
can use the semiclassical collisionless kinetic equation 
for the distribution function with the Hamiltonian (|136p . 
The results of the calculations show that, if Ep > 28, 
also in 3D the cloud is not able to oscillate around the 
new equilibrium position, but exhibits damped oscilla- 
tions around an offset point, similarly to the ID case. 
The damping is due to the fact that different atoms os- 
cillate with different frequencies as a consequence of the 
non harmonic nature of the Hamiltonian. These phenom- 
ena were investigated experimentally using a Fermi gas 
of 40 K atoms (Pezze et al, 2004). In Fig. [33] we show 
the observed time dependence of the z-coordinate of the 
cloud at T = 0.3 Tp, both without and in the presence of 
the lattice. One clearly sees the offset of the oscillations 
as well as their damping. 



4- Anti-bunching effect in the correlation function 

The discussion has so far concerned one-body proper- 
ties of the gas, like the momentum distribution and the 
density profiles. An interesting feature exhibited by spin- 
polarized Fermi gases is the occurrence of anti-bunching 
effects in the two-body correlation function [see discus- 
sion on g-\-\(r) in Sec. IV B 2| . This suppression of the 
probability of finding two atoms within short distances 
is a direct consequence of the Pauli exclusion principle. 

In a recent experiment (Rom et al., 2006) the density 
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FIG. 33: Dipole oscillations of a Fermi gas of 40 K atoms at 
T = 0.3 Tf in the presence (solid symbols and solid line) and 
in the absence (open symbols and dashed line) of a lattice 
with height s = 3. The lines correspond to the theoreti- 
cal predictions and the symbols to the experimental results. 
The horizontal dot-dashed line represents the trap minimum. 
From Pezze et al. (2004). 
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FIG. 34: Density correlation function C(d) measured after 
expansion as a function of the distance d = \z — z'\ in a 
Fermi gas of 40 K atoms released from an optical lattice. The 
difference C(d) — 1 of the correlation function with respect 
to its uncorrelated value shows a clear evidence of the anti- 
bunching effect at d = £ = 2hKt/m. From Rom et al. (2006). 
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correlation function (n (z, t) n (V, t)) of a Fermi gas was 
measured after expansion from an optical lattice. The 
average is taken on different runs of the experiment and 
the distributions are integrated along the x,y directions. 
The density correlation function, measured at large ex- 
pansion times, is proportional to the momentum corre- 
lation function (n(p = mz/t)n(p' — mz'/t)). Due to the 
periodicity of the Bloch function u qzn {z) an atom with 
quasi-momentum p = hq carries momenta with values 
h(q + 2jK) for all integers j. As a consequence the same 
atom can be found at the points z = h(q + 2jK)t/m. 
Since the Pauli principle does not allow two fermions 
to occupy the same Bloch state, the correlation func- 
tion (n (z, f) n (z' , t)) should vanish for relative distances 
which are integer multiples of 2hKt/m. 

In the experiment of Rom et al. (2006) 40 K atoms 
were confined at T/Tp w 0.23 in a 3D optical trap. In 
the conditions of the experiment the atoms filled the first 
energy band. The key experimental results are presented 
in Fig. [MJ where the anti-bunching effect is visible at 
\z — z'\ = £ = 2hKt/m. Notice that at short relative 
distances \z — z'\ =0 the anti-bunching effect is masked 
by the positive contribution of the autocorrelation term 
to the density correlation function (|6"5)) . 



B. Interacting fermions in optical lattices 

The study of interacting Fermi gases in optical lattices 
is expected to become a growing field of research. The 
problem has been so far approached theoretically using 
two different perspectives. 

In a first approach one starts from the two-body Hamil- 
tonian where interaction effects are accounted for in 
terms of a single parameter, the s-wave scattering length 
a. More microscopic details of the interatomic poten- 
tial are unimportant provided that the lattice period is 
much larger than the effective range \R*\ of the interac- 
tion. The resulting many-body theories are in most cases 
well suited to treat the superfluid phase of the system, 
but have not so far extensively developed to investigate 
other states, like the Mott insulator or the antiferromag- 
netic phase. Basic applications of this approach concern 
the study of the dimer formation (see next Section) and 
the calculation of the BCS critical temperature (Orso and 
Shlyapnikov, 2005). 

A second approach is based on the development of 
more phenomenological models, like the Hubbard model, 
extensively employed in solid state physics. The Hubbard 
model is well suited to study the novel phases emerging in 
the presence of the optical lattice, like the Mott insulator 
and the antiferromagnetic phases (Werner et al., 2005). 
Key questions are the identification of the parameters 
of the model in terms of the microscopic ingredients of 
the problem (the s-wave scatttering length and the inten- 
sity of the periodic potential) and its applicability under 
extreme conditions, for example at unitarity, where the 
scattering length is much larger than the period of the 



lattice. 

A separate class of problems finally concerns the 
physics of low dimensional systems, in particular of ID 
systems, which can be experimentally produced using op- 
tical lattices and will be discussed in the next Section. 

From the experimental point of view first important re- 
sults on the role of interaction have concerned the dimer 
formation in 3D tight lattices (Stoferle et al., 2006) and 
the superfluid to Mott-insulator transition (Chin et al, 
2006). 



1. Dimer formation in periodic potentials 

In this Section we discuss how the formation of a dimer 
is perturbed by the presence of a periodic potential. The 
problem is non trival because, differently from free space 
or harmonic trapping, the two-body problem can not be 
simply solved by separating the relative and center of 
mass coordinates in the Schrodinger equation. In partic- 
ular, the center of mass motion affects the binding energy 
of the molecule. 

The problem of calculating the binding energy for ar- 
bitrary laser intensities of ID lattices was solved by Orso 
et al. (2005) (the limit of tight lattices was previously 
considered by Fedichev, Bijlsma and Zoller, 2004). Dif- 
ferently from free space, where dimers are created only for 
positive values of the scattering length, in the presence of 
the periodic potential bound dimers exist also for nega- 
tive values of the scattering length, starting from a criti- 
cal value a cr < 0. When the laser intensity becomes very 
large the dimer enters a quasi-2D regime. In this limit 
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FIG. 35: Binding energy of molecules measured with 40 K 
atoms in a 3D optical lattice. The data correspond to dif- 
ferent intensities of the optical lattice: s — 6Er (triangles), 
s — WEr (stars), s = 15Er (circles) and s = 22Er (squares). 
The solid line corresponds to Eq. (|138[) with no free parame- 
ters. At the position of the Feshbach resonance (a — > ±oo) the 
binding energy takes the value e;, = ~hu) op t- From Stoferle et 
al. (2006). 
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the two interacting atoms are localized at the bottom 
of the same optical well where, in first approximation, 
the potential is harmonic with frequency ui pt- Then the 
two-body problem can be solved analytically yielding, in 
particular, the value = —0.244 huj opt — ~QA%%^[sEr 
at unitarity (Petrov, Holzmann and Shlyapnikov, 2000; 
Idziaszek and Calarco, 2006). 

The formation of dimcrs affects dramatically also the 
tunneling of particles through the barriers produced by 
the optical lattice, particularly the effective mass M* 
which is defined through the dispersion law E(p z ) of a 
molecule as 1/M* = [d 2 E(p z )/dp 2 z ] Pz=0 . As a result the 
effective mass M* of the dimer is significantly larger than 
the value 2m* where m* is the effective mass of a single 
atom in the presence of the same lattice potential (see 
Sec. lX Al) . The difference is caused by the exponential de- 
pendence of the tunneling rate on the mass of a tunneling 
particle. Near the threshold for the molecular formation 
M* approaches the non-interacting value 2m* . 

Let us now discuss the case of a 3D lattice potential 
of the form (|131|) . We restrict the analysis to a lattice of 
high intensity s. In this case the atomic pair is confined 
near one of the minima of the lattice, where the potential 
can be considered harmonic and isotropic with frequency 
uj pt and the two-body problem can be solved analytically 
(Bush et al., 1998). A bound state is found for any value 
of the scattering length, the binding energy being given 
by the solution of the equation 

r (-e b /2huo opt - 1/2) a 

where T is the Gamma function and a op t = yJh/mLu opt . 
The resulting predictions are shown in Fig. [35] For small 
and positive scattering lengths (a <C a opt ) Eq. (|138|) 
yields the binding energy ej, = —ti 2 /ma 2 relative to free 
space, while at unitarity one finds the result = —ftw opt . 

The formation of molecules driven by the presence of 
the lattice was observed in the experiment of Stoferle 
et al. (2006), where the lattice was formed by three 
orthogonal standing waves and the binding energy was 
measured by radio-frequency spectroscopy. The radio- 
frequency pulse dissociates dimers and transfers atoms in 
a different hyperfine state which does not exhibit a Fes- 
hbach resonance. Therefore the fragments after dissocia- 
tion are essentially non-interacting. The results are pre- 
sented in Fig. for different values of s and show good 
agreement with theory. In particular, one can clearly 
see the existence of bound states for negative values of o 
which would be impossible in the absence of the lattice. 



2. Hubbard model 

The Hubbard model (Hubbard, 1963) provides a useful 
description for atoms in the lowest band of tight lattices. 



In the simplest version the Hamiltonian has the form: 

H = -t^2 (4i d rt + 4i d ii) + U J2 "*t"a ' ( 139 ) 

(id) '< 

where the indices i and j run over the N a sites of the 
cubic lattice and correspond to first neighbor sites. The 
operators c[ a (£i a ) are the usual creation (annihilation) 
operators of particles with spin a =|, J, on the site i, while 
hia = c\ a c.i a is the corresponding number operator. The 
term in the Hamiltonian (|139[) containing the parame- 
ter t (hopping term) describes the tunneling of atoms 
between sites and plays the role of the kinetic energy op- 
erator. If U = this term gives rise to the dispersion law 
(I132p with 6 — 2t. The parameter U describes instead 
the interaction between atoms and can have both posi- 
tive and negative sign. Physically this term corresponds 
to the energy shift produced by the interaction when two 
atoms of opposite spin are localized in one of the lat- 
tice sites. The perturbative calculation of the shift using 
the pseudopotential (fl"5"j) yields U = (4irh 2 a/m) J dr|V>o| 4 
where ipo is the ground-state wavefunction of an atom in 
the individual potential well. In the harmonic approxi- 
mation for the local optical potential, holding for large 

laser intensities, one finds U — \ -Eps 3 ^ 4 Ka. The use 

of perturbation theory is justified if the shift is small 
compared to the optical oscillator energy \U\ -C frhj opt , 
or equivalently if \a\ <C a op t- This condition is equivalent 
to the requirment that the energy U is small compared 
to the gap between first and second band. For higher 
values of U the applicability of the Hubbard model (|139p 
is questionable since in this case the Hamiltonian should 
account also for higher bands. 

At zero temperature the model is characterized by two 
parameters: the ratio u = U/t between the interaction 
and the hopping coefficients and the average occupancy 
p = N/N a , which in the first band picture should be 
smaller or equal to 2 due to Fermi statistics. The phase 
diagram predicted by the Hubbard model is very rich, 
including the superfluid, the Mott insulator and the an- 
tiferromagnctic phases. It is also worth noticing that in 
the strong coupling limit the Hubbard model is equiva- 
lent to the Heisenberg model (see, for example, Bloch, 
Dalibard and Zwerger, 2007). The detailed discussion of 
the various phases available in interacting Fermi gases 
trapped by a periodic potential lies outside the scope of 
this work and we refer to the recent reviews by Georges 
(2007) and Lewenstein et al. (2007). 

XI. ID FERMI GAS 

The physical properties of ID Fermi gases differ in 
many interesting aspects from the ones of their 3D coun- 
terparts. 

In practice, to create a one-dimensional gas, atoms 
must be confined in a highly elongated harmonic trap, 
where the anisotropy parameter A = lu z /lu± is so small 
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that the transverse motion is "frozen" to the zero point 
oscillation. At zero temperature this condition implies 
that the Fermi energy associated with the longitudinal 
motion of the atoms, Ep = NHlo z /2, is much smaller 
than the separation between the levels in the trans- 
verse direction, Ef <C Huj±. This condition requires 
A <C 1/N. Such a configuration can be realized using 
a two-dimensional optical lattice formed by two perpen- 
dicular standing- wave laser fields. If the intensity of the 
beams is large enough, the tunneling between the minima 
of the lattice is absent and the atoms, confined in differ- 
ent minima, form an array of independent tubes. In the 
experiment by Moritz et al. (2005), the typical number 
of particles per tube is less than 100, while A ~ 0.004, 
thereby ensuring a reasonably safe ID condition in each 
tube. 



A. Confinement Induced Resonance 

At low energy the scattering process between two 
fermions with opposite spin colliding in a tightly confined 
waveguide (uj z = 0) can be described by the effective ID 
interaction potential 



Vid(z) = gioS(z) 



(140) 



where the coupling constant giu is expressed in terms of 
the 3D scattering length a and of the transverse oscillator 
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FIG. 36: (color online). Molecular binding energy measured 
with 40 K in ID [solid (black) symbols] and in 3D [open (blue) 
symbols] configurations. The lower (black) line corresponds to 
the theoretical prediction of Bergeman, Moore and Olshanii 
(2003) without free parameters. The upper (blue) line corre- 
sponds to the law —h 2 /ma 2 + finite-range corrections. The 
vertical dashed (blue) line represents the position of the Fes- 
hbach resonance. From Moritz et al. (2005). 



length ax = \J h/mujj_ (Olshanii, 1998) 

_ 2h 2 a 1 
ma^ 1 — Ca/a± 



(141) 



Here, C = -<(l/ 2 )/V2 ~ 1.0326, with Q(x) denoting the 
Riemann zeta function. The condition of validity for the 
effective ID interaction (|140l) - (|14ip is provided by 



< 1 



(142) 



where k z is the longitudinal relative wavevector of the 
colliding atoms. The coupling constant g\u is reso- 
nantly enhanced for a — > a C i r = a±/C, corresponding 
to the so called confinement induced resonance (CIR), 
while it remains finite at the position of the 3D res- 
onance (a — > ±oo) where it takes the negative value 
giD = —2h 2 /(Cma±). Positive values of g\r>, cor- 
responding to an effective ID repulsive potential, are 



obtained only in the interval < a < 



Other- 



wise gin < 0, corresponding to an effective attraction. 
If \a\ <C aj_, the coupling constant takes the limiting 
form gi£> = 2h 2 a/(ma 2 L ) 1 which coincides with the re- 
sult of mean-field theory where the 3D coupling con- 
stant g = 4irh 2 a/m is averaged over the harmonic oscil- 
lator ground state in the transverse direction (see, e.g., 
Pitaevskii and Stringari, 2003, Chap. 17). 

In the region where g\o is negative, two atoms can 
form a bound state. The wavefunction of the relative mo- 
tion is obtained by solving the ID Schrodinger equation 
with the potential (|140p and is given by ip{ z ) = V^e - ^ 2 , 
where k = y '—me^/h. For the binding energy ef, one finds 
the result 



Ah 2 



9lD 



(143) 



yielding n = (m/2h 2 ) \gm\- Notice that is the energy 
of a dimer relative to the non-interacting ground-state 
energy huj±. The ID result (|143[) for the binding energy 
is valid under the condition na± <C 1, or equivalently 
|efc| hui±. The general problem of calculating e& in 
a tightly confined waveguide has been solved by Berge- 
man, Moore and Olshanii (2003) using the pseudopoten- 
tial (fT"5")) . A molecular bound state exists for any value of 
the scattering length a. Its energy approaches the free- 
space result —h 2 /ma 2 for a > and a a±, and the ID 
result (|143|) if a < and \a\ a±. At the Feshbach res- 
onance, 1/a = 0, these authors find the universal result 
tb — —O.Qhw^. 

These molecular bound states have been observed with 
40 K in the experiment by Moritz et al. (2005), where the 
binding energy ej, was measured using radio-frequency 
spectroscopy. In Fig. [3()]we show the experimental results 
obtained in highly elongated traps compared with the 
quasi-lD theoretical predictions of Bergeman, Moore and 
Olshanii (2003) (see also Dickerscheid and Stoof, 2005). 
The corresponding results for the molecular binding en- 
ergy in 3D configurations are also reported in the figure, 
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explicitly showing the existence of confinement induced 
molecules in the region of negative scattering lengths. 

The three- and four-body problems concerning the 
atom-dimer and the dimer-dimer scattering in quasi- ID 
configurations has been solved by Mora et al. (2004, 
2005) using techniques similar to the 3D calculation by 
Petrov, Salomon and Shlyapnikov (2004) discussed in 
Sec. MI CI In particular, one finds that the scatter- 
ing process between ID dimers with energy (|143j) can 
be described by the contact potential of the form (|140p 
with the same atom-atom coupling constant Q\d- Since 
giB < 0, the interaction between these dimers is attrac- 
tive. Notice, however, that the fermionic nature of the 
atoms prohibits the formation of bound states with more 
than two particles. 



B. Exact theory of the ID Fermi gas 



We consider a two-component Fermi gas with equal 
populations of the spin states (iVf = N± = N/2) con- 
fined in a tight waveguide of length L. At zero tempera- 
ture and in the absence of interactions, all single-particle 
states within the "Fermi line" —kp<k<kp are occu- 
pied. The Fermi wavevector 



kp — 2 UlD 



(144) 



is fixed by the linear density nio = N/L and the corre- 
sponding Fermi energy is given by Ep = (irkum) 2 /8m. 
The condition (|142p . allowing for the use of the effec- 
tive ID interaction (|140[) - (|141[> . implies the requirement 
niDd± <C 1. In this case the many-body problem is com- 
pletely determined by the Hamiltonian 



H 



ID 



2m 



N 



^dz? 



1=1 



N r N l 

i=l i'=l 



Zi') , (145) 



which contains only one dimcnsionless parameter 

mq\ n 

(146) 



7 



h 2 ni 



D 



Correspondingly the ground-state energy per atom can 
be written in the form 



L ID 



N 



2m 



e(7) 



(147) 



in terms of the dimensionless function e(7), and analo- 
gously for the chemical potential /i = d(nip>E / N) / dn\p). 
From Eq. (|146[) one notices that the weak-coupling 
regime (|7| <C 1) corresponds to high densities nip>, while 
the strong-coupling regime (I7I ^> 1) is achieved at low 
densities. This is a peculiar feature of ID configurations, 
resulting from the different density dependence of the 
ratio of kinetic to interaction energy in ID (cx nio) com- 
pared to 3D (cx n -1 / 3 ). 



The ground-state energy of the Hamiltonian (|145|) has 
been calculated exactly using Bethe's ansatz both for re- 
pulsive, gip, > (Yang, 1967), and attractive, g\u < 
(Gaudin, 1967 and 1983), interactions. Some limiting 
cases of the equation of state at T = are worth dis- 
cussing in detail. In the weak-coupling limit, I7I <C 1, 
one finds the perturbative expansion 
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+...), 



(148) 



where the first correction to the Fermi energy Ep car- 
ries the same sign of 7. The above expansion is the ID 
analogue of Eq. (|3l?l) holding in 3D. In the limit of strong 
repulsion, 7 ^> 1, one instead finds (Recati et al., 2003b) 



/i = 4E F 1 



16 In 2 

37 



(149) 



The lowest order term in the above expansion coin- 
cides with the Fermi energy of a single component non- 
interacting gas with a twice larger density (N a = N), 
consistently with the expectation that the strong atom- 
atom repulsion between atoms with different spins plays 
the role in ID of an effective Pauli principle. 

The regime of strong attraction, \j\ ^> 1 with 7 < 0, is 
particularly interesting. In this case one finds the follow- 
ing expansion for the chemical potential (Astrakharchik 
et al., 2004b) 



P 2 



e b , E 



4 

,1 

4 V 3 7 



(150) 



where eb is the binding energy (|143[) of a dimer. Ac- 
cording to the discussion at the end of Sec. IXI Al in 
this regime the system is a gas of strongly attractive 
bosonic dimers. A crucial point is that the formation of 
bound states of these composite bosons is inhibited by the 
fermionic nature of the constituent atoms. In Eq. (|150p 
the leading term Ep /4, beyond the molecular binding en- 
ergy eb, coincides with half the chemical potential of a gas 
of impenetrable bosons (Tonks- Girardeau gas) with den- 
sity n\pj/2 and mass 2m. This peculiar behavior is an ex- 
ample of the exact mapping between bosons and fermions 
exhibited by ID configurations (Girardeau, 1960). The 
consequences of the "Fermi-Bose duality" in ID, includ- 
ing the regimes of intermediate coupling and for both 
s-wave and p-wave scattering, have been investigated in 
a series of papers (Cheon and Shigehara, 1999; Granger 
and Blume, 2004; Girardeau, Nguyen and Olshanii, 2004; 
Girardeau and Olshanii, 2004; Mora et al., 2005). 

The strongly interacting regime, I7I 3> 1, can be 
achieved by tuning the effective coupling constant g\D- 
As in the experiment by Moritz et al. (2005) on ID 
molecules, one makes use of a Feshbach resonance to 
tune the scattering lenth in the region around the value 
a c i r corresponding to the confinement induced resonance 
where 7 = ±00 [see Eq. (|141[) ]. This resonance con- 
nects a BCS-like weakly attractive regime (7 negative 
and small), corresponding to weakly bound pairs with 
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energy (|143p and size kT 1 >> nTh, to a BEC-like regime 
(7 positive and small) of tightly bound bosonic molecules 
with energies of order fku± and size much smaller than 
the average distance between dimers. These dimers are 
expected to behave as a ID gas of bosons interacting with 
a repulsive contact potential (Mora et ai, 2005). How- 
ever, these tightly bound dimers can not be described by 
the Hamiltonian (|145|) which, if gip, > 0, only describes 
the repulsive atomic branch. The occurrence of such a 
BCS-BEC crossover in ID was first suggested by Tokatly 
(2004) and by Fuchs, Recati and Zwerger (2004). 

The properties at T = of the ID Fermi gas with 
short-range interactions are considerably different from 
the ones of usual 3D Fermi liquids. Actually, this gas is 
an example of a Luttinger liquid (Luttinger, 1963). The 
low-energy properties of this liquid are universal and do 
not depend on the details of the interaction, on the spe- 
cific Hamiltonian (lattice or continuum models) nor on 
the statistics of the atoms (Haldane, 1981). The only 
requirement is the existence of long-wavelength gapless 
excitations with linear dispersion. The Luttinger effec- 
tive Hamiltonian is expressed in terms of the compress- 
ibility and of the velocity of propagation of these gapless 
excitations. For a review on the properties of Luttinger 
liquids see, e.g., Voit (1995) and Giamarchi (2004). 

Hydrodynamic sound waves (phonons) are predicted 
by the Hamiltonian (|145|) . They propagate with the ve- 
locity c determined by the compressibility through the 
general relation mc 2 = niE>dn / 1 dn\D . In this case, the 
knowledge of the equation of state allows for an exact 
determination of the effective Luttinger parameters (Re- 
cati et ai, 2003a). If 7 > 0, c is larger than the Fermi 
velocity vp — Tthnipj/2m. For example, in the case of 
strong repulsion (7 > 1) the speed of sound takes the 
limiting value c = 21^(1 — 4 In 2/7). F° r small values 
of 7 the sound velocity tends to «f, while for 7 < it 
becomes smaller than the Fermi velocity. The inverse 
compressibility mc 2 , however, remains positive indicat- 
ing the stability of the gas even in the strongly attractive 
regime where one finds c = vp/2. This is in sharp con- 
trast with the behavior of a ID Bose gas with attractive 
contact interactions, where the ground state is a soliton- 
like many-body bound state (McGuire, 1964). 

The presence of phonons in a Luttinger liquid af- 
fects dramatically the long-range behavior of the cor- 
relation functions, fixed by the dimensionless parameter 
i] = 2hkp/mc (Luther and Peschel, 1974; Haldane, 1981). 
For example, the one-body density matrix behaves, for 
\z — z'\ > 1/niD, as 



(*t {z )y a (/))« 



ID \Z 



1 I J 

z'\)~ + ~ 



(151) 



For a non-interacting gas 77 = 2 and the one-body density 
matrix decays as sin (Tmiu\z — z'\) j\z — z \. This behav- 
ior reflects the presence of the jump from 1 to in the mo- 
mentum distribution at the Fermi surface k = ±kp. In 
the presence of interaction the correlation function (|151|) 



decreases faster. This implies that the jump at kp disap- 
pears and the momentum distribution close to the Fermi 
surface behaves as — n^ F cx sign(kp — k)\kp — fc|' 3 , 
with /3=i + 2- l>0. 

The Hamiltonian (|145[) supports also spin waves to- 
gether with sound waves. For 7 > these spin excitations 
also have a linear dispersion at small wavevectors. Their 
velocity of propagation c s tends to vf for small values of 
7, but for finite interaction strengths it is different from 
the speed of sound c. In the strongly repulsive regime 
one finds the small velocity c s = vpir 2 /7, to be compared 
with the corresponding sound velocity c = 2vp discussed 
above. This "spin-charge separation" is a peculiar fea- 
ture of Luttinger liquids. The possibility of observing 
this phenomenon in ultracold gases has been investigated 
by Recati et al. (2003a). In the case of attractive interac- 
tions, 7 < 0, the spin-wave spectrum exhibits a gap A gap 
(Luther and Emery, 1974). This "spin-gap" is defined ac- 
cording to Eq. (|30|) and is therefore analogous to the pair- 
ing gap in 3D Fermi superfluids. In the weak-coupling 
limit I7I -C 1, the gap is exponentially small being pro- 
portional to A/pyfexp (— 7r 2 /2|7|) (Bychkov, Gorkov and 
Dzyaloshinskii, 1966; Krivnov and Ovchinnikov, 1975). 

So far we have considered uniform systems. In experi- 
ments the gas is confined in the longitudinal z-direction 
by a harmonic potential. If the average distance between 
particles is much smaller than the longitudinal oscillator 
length, X/niD <C a z — y/ h/moj z (requiring small enough 
values of the trapping frequency lo z ), one can use the lo- 
cal density approximation to calculate the properties of 
the trapped system (Astrakharchik et al., 2004b). The 
density profile and the Thomas-Fermi radius can be ob- 
tained as a function of the interaction strength from the 
solution of Eq. (|T2"j) . One can also calculate the frequency 
of the lowest compression mode using the hydrodynamic 
theory of superfluids of Sec. IVII CI Since fi cx n\ D both 
in the weak- and strong-coupling limit, independent of 
the sign of interactions [see Eqs. (|149|) and (|150p ]. the 
frequency of the mode tends to 2lo z in these limits. The 
study of the breathing-mode frequency for intermediate 
couplings has been carried out by Astrakharchik et al. 
(2004b). 

Another interesting application of the local density ap- 
proximation to trapped ID configurations is provided by 
the study of spin polarized systems (Orso, 2007). The 
ground-state energy of the Hamiltonian (|145[) can be 
calculated exactly also with unequal spin populations 
A T ^ Ni (Yang, 1967; Gaudin, 1967). In the case of 
attractive interactions, g\o < 0, the T = phase dia- 
gram always consists of a single phase corresponding to: 
i) a fully paired state with a gap in the spin excitation 
spectrum if 2Vj = N± , ii) a non-interacting fully polarized 
state if N± — and iii) a gapless partially polarized state 
if iVf > Ni (Oelkers et ai, 2006; Guan et al, 2007). 
This latter state is expected to be a superfluid of the 
FFLO type (Yang, 2001). In harmonic traps spin un- 
balanced configurations result in a two-shell structure: a 
partially polarized phase in the central region of the trap 
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and either a fully paired or a fully polarized phase in the 
external region depending on the value of the polariza- 
tion (Orso, 2007). This structure is in sharp contrast 
with the behavior in 3D configurations, where the unpo- 
larized superfluid phase occupies the center of the trap 
and is surrounded by two shells of partially polarized and 
fully polarized normal phases (see Sec. IIX Cp . One can 
understand this behavior by noticing that the larger den- 
sities occurring in the center of the trap correspond in ID 
to a weak-coupling regime and, consequently, pairing ef- 
fects are smaller in the central than in the external region 
of the trap. The opposite situation takes place instead 
in 3D configurations. 



XII. CONCLUSIONS AND PERSPECTIVES 

In this review we have discussed some relevant fea- 
tures exhibited by atomic Fermi gases from a theoreti- 
cal perspective. The discussion has pointed out a gen- 
eral good agreement between theory and experiment, re- 
vealing that the basic physics underlying these quantum 
systems is now reasonably well understood. The most 
important message emerging from these studies is that, 
despite their diluteness, the role of interactions in these 
quantum gases is highly non trivial, revealing the dif- 
ferent facets of Fermi superfluidity in conditions that are 
now accessible and controllable experimentally. The pos- 
sibility of tuning the value and even the sign of the scat- 
tering length is actually the key novelty of these sys- 
tems with respect to other Fermi superfluids available in 
condensed matter physics. The long sought BCS-BEC 
crossover, bringing the system into a high T c superfluid 
regime where the critical temperature is of the order of 
the Fermi temperature, can be now systematically inves- 
tigated and many theoretical approaches are available to 
explore the different physical properties. The situation is 
particularly well understood at zero temperature where 
important properties of these Fermi gases in harmonic 
traps, like the equilibrium density profiles, the values 
of the release energy and of the collective frequencies, 
can be calculated in an accurate way using the equation 
of state of uniform matter available through Quantum 
Monte Carlo simulations and employing the local den- 
sity approximation. Other important features that can 
be now considered reasonably well understood theoret- 
ically and confirmed experimentally are the momentum 
distribution along the crossover, the collisional processes 
between pairs of fermions on the BEC side of the reso- 
nance and the basic properties of Fermi gases in optical 
lattices such as Bloch oscillations, the structure of the 
Fermi surface and the binding of molecules. In general 
the most interesting regime emerging from both the the- 
oretical and experimental investigations is the unitary 
limit of infinite scattering length, where the dilute gas 
becomes strongly correlated in conditions of remarkable 
stability. In this regime there are no length scales related 
to interactions entering the problem, which consequently 



exhibits a universal behavior of high interdisciplinary in- 
terest. 

Many important issues remain nevertheless to be ad- 
dressed or to be explored in a deeper and systematic way. 
A brief list is presented below: 

- Thermodynamics. More theoretical work needs to 
be done to determine the transition temperature along 
the crossover and the thermodynamic functions below 
and above T c . Very little is known about the tempera- 
ture dependence of the superfluid density (Taylor et al, 
2006; Akkineni, Ceperley and Trivedi, 2006; Fukushima 
et al., 2007) and its role on physically observable quan- 
tities like, for example, the propagation of second sound 
(Taylor and Griffin, 2005; Heiselberg, 2006; He et al, 
2007). Another open question remains the identification 
of a good thermometry in these ultracold systems, where 
the experimental value of the temperature is often sub- 
ject to large uncertainties. 

- Collective modes and expansion. The transition 
from the hydrodynamic to the collisionlcss regime on the 
BCS side of the resonance and the consequences of the 
superfluid transition on the frequencies and the damping 
of the collective oscillations, as well as on the behavior 
of the aspect ratio during the expansion, still requires a 
clear understanding. This problem raises the question of 
the inclusion of mesoscopic effects in the theoretical de- 
scription which become important when the pairing gap 
is of the order of the harmonic oscillator energy. Also 
the question of the transition at finite temperature be- 
tween the hydrodynamic and the collisionless regime in 
the normal phase near resonance requires more theoreti- 
cal and experimental work. Another important question 
is the temperature dependence of the viscosity of the gas, 
which is in principle measurable through the damping of 
the collective modes. The problem is particularly inter- 
esting at unitarity where a universal behavior is expected 
to occur (Son, 2007). 

- Diabatic transformation of the scattering 
length. The experimental information on the conden- 
sation of pairs (Regal, Greiner and Jin, 2004b; Zwier- 
lein et al., 2004) and on the structure of vortices (Zwier- 
lein et al., 2005b) along the crossover is so far based on 
the measurement of the density profiles following the fast 
ramping of the scattering length before expansion. The 
theoretical understanding of the corresponding process 
is only partial and requires more systematic investiga- 
tions through the implementation of a time-dependent 
description of the many-body problem. 

- Coherence in Fermi superfluids. Although 
the experimental measurement of quantized vortices has 
proven the superfluid nature of these ultracold gases, 
coherence effects have not yet been explored in an ex- 
haustive way. For example, the implications of the order 
parameter on the modulation of the diagonal one-body 
density in interference experiments remain to be stud- 
ied. Other relevant topics are the interference effects in 
the two-body correlation function (Carusotto and Castin, 
2005) and the Josephson currents (Spuntarelli, Pieri and 
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Strinati, 2007). 

- Rotational properties. The physics of interact- 
ing Fermi gases under rotation in harmonic traps is a 
relatively unexplored subject of research, both in the ab- 
sence and in the presence of quantized vortices. Ques- 
tions like the difference between the superfluid and the 
collisional hydrodynamic behavior in rotating configura- 
tions, the Tkachenko modes of the vortex lattice (Watan- 
abe, Cozzini and Stringari, 2007) and the nature of the 
phase diagram of the rotating gas at high angular veloc- 
ity are topics requiring both theoretical and experimental 
investigation. 

- Spin imbalanced Fermi gases. This subject of re- 
search has attracted a significant amount of work in the 
last few years. Many important questions still remain 
open, like the nature of the superfluid phases caused by 
the polarization at zero and at finite temperature in the 
different regimes along the crossover. Recent experiments 
have also raised the question of the role of surface ten- 
sion effects at the interface between the normal and the 
superfluid component (De Silva and Mueller, 2006b; Par- 
tridge et al, 2006b) and of correlation effects in the spin 
polarized normal phase (Schunck et al., 2007) 

- Fermi- Fermi and Fermi-Bose mixtures. A grow- 
ing activity, both on the experimental and on the theoret- 
ical side, is expected to characterize future studies of mix- 
tures of different atomic species. The quantum phases of 
Fermi-Bose mixtures in optical lattices (Albus, lllumi- 
nati and Eisert, 2003; Lewenstein et al, 2004; Gunter et 
al, 2006), the formation of dipolar gases (Ospelkaus et 
al., 2006c; Modugno, 2007) and the superfluid behavior 
of Fermi-Fermi mixtures of different atomic masses are 
important topics for future research. 

- RF transitions. RF transitions provide valuable 
information on the gap parameter and on pairing effects 
in interacting Fermi gases (Chin et al., 2004; Shin et al., 
2007) . The proper inclusion of final-state interactions in 
the calculation of the spectral response is a crucial ingre- 
dient to make the analysis of the experimental findings 
conclusive on a quantitive basis. 

- Quantum impurities. The investigation of the mo- 
tion of impurities added to a trapped quantum gas can 
open interesting possibilities for the determination of the 
viscosity coefficients and for the study of the Landau cri- 
terion of superfluidity. 

- Phase transitions in optical lattices. The recent 
experiments on the superfluid to Mott-insulator phase 
transition (Chin et al., 2006) have already stimulated first 
theoretical work involving the use of a multi-band Hub- 
bard model (Zhai and Ho, 2007). An important issue is 
the behavior of the transition at unitarity where the scat- 
tering length is much larger than the lattice spacing of 
the optical potential. Further important topics concern 



the study of the various magnetic phases, which can be 
implemented by playing with the geometry and the di- 
mensionality of the lattice, as well as the role of disorder 
(Lewenstein et al, 2007). 

- p-wave superfluidity. The recent experiment by 
Gacblcr et al. (2007) on the production and detection of 
molecules in a single-component 40 K gas using a p-wave 
Feshbach resonance and the measurement of the binding 
energy and lifetime of these p-wave molecules opens new 
exciting possibilities of realizing p-wave superfluids with 
ultracold gases. First theoretical investigations of the 
many-body properties of these systems as a function of 
the interaction strength predict a much richer phase dia- 
gram compared to s-wave superfluids, including quantum 
and topological phase transitions (Gurarie, Radzihovsky 
and Andreev, 2005; Cheng and Yip, 2005; fskin and Sa 
de Melo, 2006a; Gurarie and Radzihovsky, 2007). 

- Low-dimensional configurations. Two- and one- 
dimensional configurations of ultracold gases can be eas- 
ily produced in laboratories. The experiment by Moritz 
et al. (2005) on ID molecules provides the first exam- 
ple of a low-dimensional arrangement combined with the 
use of a Feshbach resonance. Many interesting features 
are expected to take place when the scattering length 
becomes larger than the characteristic length of the con- 
finement. The properties of the BCS-BEC crossover in 
2D and in ID as well as the physics of Luttinger liquids 
in ID can be addressed by further investigations. 
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